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ABSTRACT 


R.  L.  Sielken  Jr.,  H.  0.  Hartley,  R.  K.  Spoerl 


Statistical  PERT  la  a new  procedure  for  obtaining  Information 
about  the  distribution  of  a project's  completion  time  when  the  project 
Is  comprised  of  a large  number  of  activities  and  the  time  required  to 
complete  an  individual  activity  once  it  can  be  begun  Is  a random 
variable.  The  project  Is  Represented  as  an  acyclic  network  whose  arcs 
correspond  to  the  project  activities.  This  network  Is  simplified  by 
replacing  various  activity  configurations  by  single  equivalent 
activities  and  then  decomposed  into  several  subnetworks.  The 
distribution  and  moments  of  each  subnetwork's  completion  time  are 
bounded  and  approximated  on  the  basis  of  two  points  from  each 
activity's  completion  time  distribution  by  using  some  mathematical 
programming  techniques  and  a new  result  in  the  theory  of  networks. 

The  project's  completion  time  distribution  is  then  approximated  by 
combining  the  approximate  subnetwork  distributions. 

This  report  describes  several  refinements  in  the  subnetwork 


I 


analysis  procedure.  One  major  refinement  greatly  reduces  the 
computational  effor^^^obtainlng  bounds  on  the  project  completion 


time  moments  and  distribution.  A second  major  refinement  allows  the 
two-point  approximation  of  an  activity's  completion  time  distribution 


to  better  represent  skewed  distributions.  The  computer  programs 
required  to  implement  the  new  subnetwork  analysis  procedure  ere  listed 
and  documented. 


-c  ~ 


Statistical  PERT:  An  Improved  Subnetwork  Analysis  Procedure 


R.  L.  Sielken  Jr.,  H.  0.  Hartley,  R.  SC.  Spoerl 

The  well-known  Program  Evaluation  and  Review  Technique  (PERT)  Is 
concerned  with  a 'project*  comprised  of  a large  number  of  'activities'  which 
are  arranged  as  the  arcs  In  a complex  acyclic  network  (sea  s.g.  Figure  1). 
The  activities  at  any  network  node  'commence'  as  soon  as  all  activities 
'terminating'  at  that  node  are  completed.  The  time  required  to  complete 
an  activity  once  It  can  be  begun  Is  a random  variable,  and  hence  the  time 
needed  to  complete  the  entire  project  Is  also  a random  variable. 

In  Technical  Report  No.  48  "Statistical  Critical  Path  Analysis  In 
Acyclic  Networks:  Statistical  PERT"  a comprehensive  new  procedure  for 

obtaining  information  on  the  project  completion  time  and  Its  distribution 
was  described  and  Illustrated.  That  procedure  involved  the  following  five 
general  steps: 

Step  1:  Identification 

Represent  the  project  and  its  component  activities  In  tarms  of  an 
acyclic  network  with  one  source  and  one  sink.  Identify  each  activity's 
completion  time  distribution  or  at  least  two  points  on  each  activity's 
completion  time  distribution. 

Step  2:  Simplification 

Replace  various  activity  configurations  and  their  associated 
completion  time  distributions  by  a single  equivalent  activity  and 
completion  time  distribution. 


Step  3:  Decomposition 

Decompose  the  simplified  network  into  ssversl  subnetworks  by 
separating  parallel  subnetworks  and  than  separating  the  resulting  sub- 
networks at  each  cut  vertex.  A cut  vertex  is  any  node  such  that  every 
path  from  the  source  to  the  sink  passes  through  it. 

Step  4:  Analysis 

Each  subnetwork  arising  from  Step  3 is  analyzed  on  the  basis  of 
two  points  from  each  component  activity's  completion  time  distribution. 
The  result  of  this  analysis  is  an  approximation  of  aach  subnetwork's 
completion  time  distribution  and  the  moments  of  this  distribution. 

Step  5:  Synthesis 

Confclne  the  approximate  subnetwork  completion  time  distributions 
determined  in  Step  4.  Hie  result  la  an  approximate  completion  tima 
distribution  for  the  entire  project. 

The  purpose  of  this  report  is  to  document  several  refinements 
in  the  subnetwork  analysis  step.  Step  4. 

1.  Analysis  of  a Subnetwork 

The  analytical  procedure  described  in  this  section  yields  the 
following  Information  on  each  subnetwork  when  each  component  activity's 
cospletion  time  distribution  is  replaced  by  a discrete  two-point 
distribution: 

(a)  Upper  end  lower  bounds  on  the  mean  subnetwork  completion  time 
as  well  as  the  other  momenta  of  the  subnetwork  completion  time. 

(b)  Upper  and  lower  bounds  on  the  distribution  function  of  the 
subnetwork  completion  time. 
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(c)  An  approximate  distribution  function  of  the  subnetwork 
completion  time. 

Each  subnetwork  is  assumed  to  be  an  acyclic  network  with  one 
source,  one  sink,  and  no  cut  vertices. 

The  analysis  of  each  subnetwork  involves  essentially  two  parts: 

1.  The  formation  of  "clusters"  of  activities  whose  effect  on  the 
subnetwork  completion  time  seems  to  be  interrelated. 

2.  The  approximation  of  the  subnetwork  completion  time  moments 
and  distribution  on  the  basis  of  the  clusters. 

1.1  Formation  of  Clusters 

The  actual  completion  time  distribution  of  each  individual  activity, 

A,  in  the  subnetwork  is  replaced  by  a discrete  distribution  with  probability 

P at  the  lower  point,  l , and  Q • 1-P  at  the  upper  point,  u . 

A A 

Let  n be  the  number  of  activities  in  the  subnetwork.  Then  for 
each  of  the  2n  combinations  of  the  l^'a  and  u^'s  there  will  be  a 
subnetwork  completion  time  (a  critical  path  time).  The  r-th  moment  of 
these  2n  times  will  be  denoted  by  Tf , and  the  distribution  function  of 
these  times  trill  be  denoted  by  F.  The  approximation  of  the  T^'s 
(especially  T^,  the  mean)  and  F is  the  goal  of  the  subnetwork  analysis. 

Since  n will  usually  be  fairly  large,  the  collate  enumeration  of  the 
2n  critical  path  timeB  will  usually  be  unreasonable.  Hence  the 
activities  which  are  most  likely  to  be  on  the  critical  path  through 
the  subnetwork  are  identified  and  their  joint  behavior  investigated. 

The  mean  of  the  completion  time  distribution  for  activity  A is 
defined  to  be  m^  - Pi^  + Qu^ . The  standard  deviation  of  the  completion 
time  distribution  for 


activity  A is  defined  to  be  s. 


r2  2 2 
nl  + (* A-«A 


and  is  assumed  to  be  positive.  (The  assumption  that  sA  > 0 is  not  really 
a practical  restriction  since  the  difference  between  a fixed  activity 
completion  time  and  one  with  a very  small  dispersion  is  negligible  from 
a practical  viewpoint.)  The  ai&network's  critical  path  when  each 
activity's  completion  time  is  set  equal  to  its  mean  will  be  referred  to 
as  the  "original"  critical  path.  The  activities  on  this  critical  path 
will  be  referred  to  as  "critical  activities"  with  K equalling  the  number 
of  such  activities.  Some  non-critical  activities  might  become  critical 
if  some  of  the  completion  times  for  the  original  critical  activities 
were  decreased.  These  activities  are  identified  as  follows.  The  completion 
time  for  one  critical  activity,  say  A,  is  set  equal  to  max{mA  - XsA,  0} 
where  A is  a non-negative  algorithm  parameter  which  the  user  specifies. 

All  other  completion  times  are  set  equal  to  their  means.  Then  the  longest 
path  through  the  resulting  network  is  determined.  Any  activities  on  this 
path  which  were  not  on  the  original  critical  path  are  now  referred  to  as 
the  "associates"  of  A since  the  effect  of  these  "associates"  on  the  networks' 
completion  time  is  related  to  A's  completion  time.  This  procedure  is 
repeated  for  each  original  critical  activity. 

Each  critical  activity  and  its  associates  make  up  one  "cluster". 

These  K initial  clusters  are  now  "pooled"  by  combining  any  two  clusters 
with  at  least  one  activity  in  cosmon.  In  general  there  will  still  be 
more  than  one  cluster,  and  many  of  the  n - K non-critical  activities  will 
not  occur  in  any  cluster. 

The  associates  correspond  to  the  activities  which  become  critical 
when  the  completion  times  of  the  original  critical  activities  are 
lowered.  However,  some  of  the  originally  non-critical  activities  may 
also  become  critical  if  their  completion  times  exceed  their  means 


and  the  completion  times  of  the  original  critical  activities  are  at 
their  means.  These  activities  are  identified  next.  Each 
originally  non-crltical  activity  is  investigated  separately.  If  activity 
A is  being  investigated,  then  the  completion  time  for  A is  set  equal  to 
m^  + 6sa  where  0 is  a non-negative  algorithm  parameter  which  the  user 
specifies.  The  completion  times  for  all  other  activities  are  set  equal 
to  their  means,  and  the  corresponding  critical  path  determined.  This 
critical  path  will  either  be  the  original  critical  path  or  a new  path 
which  includes  A.  In  the  latter  case,  the  activities  on  the  original 
critical  path  which  are  not  on  a new  critical  path  containing  A are 
called  the  "ellmlnants"  of  A.  Thus,  the  effect  of  A's  ellminants  on 
the  networks  completion  time  is  related  to  the  completion  time  for  A. 
Hence,  A is  added  to  any  cluster  which  contains  at  least  one  of  A's 
ellminants.  After  this  procedure  has  been  repeated  for  each  originally 
non-crltical  activity,  the  resultant  clusters  are  "pooled"  again  by 
combining  any  two  clusters  with  at  least  one  activity  in  common. 

Although  the  number  of  clusters  is  reduced  when  the  pooling  on 
the  basis  of  the  associates  occurs  and  then  further  reduced  when  the 
pooling  on  the  basis  of  the  ellminants  occurs,  there  will  generally 
remain  more  than  one  cluster  and  several  activities  not  in  any  cluster. 

In  general  the  larger  the  values  of  X and  6 the  greater  the 
number  of  activities  in  the  clusters  and  the  smaller  the  number  of 
clusters.  In  particular  the  procedure  for  forming  the  clusters  has  the 
following  properties: 

Property  1:  If  X^  > X^,  then  any  activity  which  would  be  an  associate 

of  a critical  activity  A when  X ■ X^  would  also  be  an 
associate  of  A when  X ■ X^. 


. m-.  m ■ 


Property  2:  If  62  > 0^,  then  any  critical  activity  which  would  be  an 

eliminant  of  a non-critical  activity  A when  6 ■ 0^  would 
also  be  an  eliminant  of  A when  8 - 82* 

Property  3:  For  any  originally  non-critical  activity  A there  exists 

8^  such  that  A will  have  some  elloinants  for  any  8 >_  0^. 
Property  4:  For  any  fixed  value  of  X,  the  set  of  activities  in  the  union 
of  the  clusters  is  monotically  non-decreasing  as  8 ■+•  ®. 
Property  5:  There  exists  a finite  value  6*  such  that  if  6 0*,  then 

every  activity  would  be  in  some  cluster. 

Property  6:  The  number  of  clusters,  originally  K,  is  non-increasing  as 


Property  7:  There  exists  a finite  value  8*  such  that  if  0 > 8*,  then 

there  would  only  be  one  cluster. 

Most  of  the  properties  of  the  cluster  formation  procedure  are 
fairly  straightforward;  however,  Property  7 requires  some  special 
justification.  This  justification  is  based  on  the  following  definition 
and  theorem  which  is  proven  in  Appendix  A. 


Definition:  In  any  acyclic  network  a bridge  over  any  two  consecutive 

arcs  A^  and  A2  is  any  arc  A^  such  that  all  paths  from  the  source  to 
the  sink  passing  through  A^  do  not  pass  through  either  A^  or  A2. 

Theorem  1:  In  any  acyclic  network  with  no  cut  vertices  there  is  at  least 

one  bridge  for  any  pair  of  consecutive  arcs. 

Property  3 Implies  that  all  activities  will  belong  to  some  cluster  if 
0 > 8*  and 


8*  ■ max{0  : A originally  non-critical). 
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Now  consider  any  two  consecutive  activities  A^  and  on  the  original 
critical  path.  Theorem  1 implies  that  there  is  a bridge  over  and  A^, 
say  A^.  Since  the  original  critical  path  passes  through  A^  and  A^»  A^ 
cannot  be  on  the  original  critical  path.  Therefore,  if  0 >_  0»  ^ 0^  , A^ 
and  A^  will  be  eliminants  of  A^  and  hence  will  be  in  the  same  cluster  as 
A^.  Thus,  since  each  cluster  contains  at  least  one  original  critical 
activity  and  any  two  consecutive  original  critical  path  activities  belong 
to  the  same  cluster  when  8 >_  0*,  there  is  only  one  cluster  when  0^0* 
and  Property  7 is  established. 


1.2  Approximate  Subnetwork  Coimletlon  Time  Moments  and  Distribution 


1.2.1  A Lower  Bound  on  T^  and  an  Upper  Bound  on  F 


For  each  cluster  C let  n denote  the  number  of  activities  in  C,  and 

c 

n n 

let  v ■ 1,  ...,  2 c index  the  2 c configurations  of  activity  completion 
times  when 

(a)  the  completion  time  for  each  activity  A not  in  C is  equal  to 
its  lower  point,  and 

(b)  the  completion  times  for  the  activities  in  C are  at  each  of 

n 

the  2 c possible  combinations  of  their  upper  and  lower  points. 


tv  - critical  path  time  for  the  v-th  configuration 


Pv  - probability  of  the  v-th  configuration 


where 


■ A <pi(1  - !v,i>  + "i  Vi’ 


a 


Then 


and 


&v  i m 1 if  the  tine  for  the  i-th  activity  in  C is  in 
the  v-th  configuration 

- 0 if  the  tine  for  the  i-th  activity  in  C ia  in  the 
v-th  configuration. 


T (C)  = l p t 
r , v v 

v»l 


f (t;  c)  = z P it(ty) 

V”1 


where 

I_(t  ) - 1 if  t < t 
t v v - 

- 0 if  t > t 
v 

n 

are  the  r-th  moment  of  the  2 c critical  path  times  and  their  distribution 
function  respectively.  Let 


and 


T~(0,  X)  - max  T"(C) 
C 


F+  (t;  8,  X)  » min  F+  (t;  C) 
C 


which  depend  on  0 and  X since  the  composition  and  number  of  clusters 
depend  on  6 and  X. 

The  first  step  in  showing  that  1^(8,  X)  is  a lower  bound  for  is 
proving  the  following  theorem: 

't. 

Theorem  2:  For  any  cluster  C,  any  positive  integer  r , and  any  activity 

A not  in  C, 


fr”(C  U {A})  > T"(C) 


Proof : Consider  any  particular  critical  path  for  a particular  one  of 

n 

the  2 c combinations  of  upper  and  lower  points  involved  in  (C) . Consider 
the  following  two  cases: 

(i)  activity  A with  its  completion  time  equal  to  l A is  on  the 
critical  path,  and 

(ii)  activity  A with  its  completion  time  equal  to  is  not  on 
the  critical  path. 

Let  $A  ■ (uA  - &A)  and  the  particular  critical  path  time  be  t.  Then 
in  case  (i) 

(a)  if  the  completion  time  for  A is  set  equal  to  uA>  this  will 
Increase  the  r-th  moment  of  the  critical  path  time  by 
Q[(t  + «A)r  - trJ;  and 

(b)  if  the  completion  time  for  A is  set  equal  to  this  will  not 
alter  the  r-th  moment  of  the  critical  path  time. 

In  case  (il) 

(a)  if  the  completion  time  for  A is  set  equal  to  uA,  this  may 

increase  the  r-th  moment  of  the  critical  path  time  by 

Q [ ( t + 6.)r  - tr]  or  less,  and 
A 

(b)  if  the  completion  time  for  A is  set  equal  to  &A,  this  will  not 
alter  the  r-th  moment  of  the  critical  path  time. 

Therefore  in  either  case  the  contribution  to  T^CC  U {A})  - Tf(C)  for  this 
particular  critical  path  will  be  between  0 and  Q[(t  + 6A>r  - trJ.  Since 
the  contribution  is  non-negative  for  each  particular  combination, 

T'(Cl)  (Al)  >T'(C).  QED 

A straightforward  application  of  Theorem  2 yields  the  following 


theorem: 
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Theorem  3:  For  any  two  clusters  and  and  any  positive  integer  r, 

VC1  U C2)  - max^r(Cl)’  *r(C2)}* 

Property  2 of  the  cluster  formation  procedure  implies  that  if  e is 
Increased  the  clusters  expand  or  are  pooled.  Thus,  Theorems  2 and  3 
imply  that,  for  fixed  X,  T~(0,  x)  is  non-decreasing  as  0 increases. 
Furthermore,  Properties  5 and  7 together  imply  that  for  0 sufficiently 
large  there  is  only  chie  cluster  and  all  of  subnetwork  activities  are  in 

t * 

that  cluster.  Hence,  for  0 sufficiently  large  1^(0,  X)  ■ T^,  and  the 
following  theorem  is  true: 

Theorem  4: 

(a)  T~(0,  x)  is  a non-decreasing  function  of  0 for  any  fixed 
values  of  X and  r; 

(b)  there  exists  a finite  value  0*  such  that  0 :>  0*  implies 

T“(0,  X)  - Tr 

for  any  X,  and  r;  and 

(c)  for  any  0,  X,  and  r 

T~(0,  X)  < Tr. 

Similarly,  the  first  step  in  showing  that  F+(t;  0,  X)  is  an  upper 
bound  on  F(t)  is  the  following  theorem: 

Theorem  5:  For  any  cluster  C,  any  value  of  t,  and  any  activity  A not 

F+(t;  C U {A})  < F+(t;  C). 


in  C 


-vit* 


3 

'j  ' ^ ..  * ,.  ..  ...... 

i* 


I'! 

f 


Proof : Consider  any  particular  configuration  of  activity  times  used  to 

determine  F+ (t;  C)  before  C is  augmented  by  A.  When  A is  added  to  C, 
this  particular  configuration  will  appear  once  with  A at  its  upper 
percentile  and  once  with  A at  its  lower  percentile.  When  A is  at  its 
lower  percentile,  the  configuration's  critical  path  time  is  unchanged. 
However,  when  A is  at  its  upper  percentile,  the  configurations'  critical 
path  time  is  either  unchanged  or  possibly  increased  if  A were  on  the 
configuration's  critical  path.  Thus,  the  addition  of  A leaves  the 

( 

cumulative  probability  associated  with  critical  path  times  less  than  or 
equal  to  t either  unchanged  or  decreased.  QED 

A straightforward  extension  of  Theorem  5 ia  the  following 
theorem: 


i 


Theorem  6:  For  any  two  clusters  and  C 2 and  any  t, 

F+(t;  Cx  U C2)  < min{F+(t;  Cjh  F+(t;  C2> } . 

Since  Property  2 of  the  cluster  formation  procedure  implies  that 
the  clusters  expand  or  are  pooled  if  8 is  increased.  Theorem  5 and  Theorem 
6 together  imply  that,  for  all  t and  any  fixed  X,  F+(t;  6,  X)  is  non- 
increasing function  of  8.  Furthermore,  Properties  5 and  7 together  imply 
that  for  8 sufficiently  large  there  is  only  one  cluster  and  all  of  the 
subnetwork  activities  are  in  that  cluster.  Hence,  for  8 sufficiently 
large  F+(t;  8,  X)  ■ F(t)  and  the  following  theorem  is  true: 

Theorem  7: 

(a)  F+(t;  0,  X)  is  a non-increasing  function  of  8 for  every  t and 

any  X; 
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(b)  there  exists  a finite  value  9*  such  that  0 _>  0*  Implies 

F+(t;  9,  X)  - F(t) 

for  every  t and  X;  and 

(c)  for  any  0,  X,  and  t 

f+  (t;  e,  x)  > F(t). 

1.2.2  An  Upper  Bound  on  and  a Lower  Bound  on  F 

For  each  cluster  C let  n denote  the  number  of  activities  In  C.  and 

c 

n n 

let  v - 1 2 c Index  the  2 c configurations  of  activity  completion 

times  when 

(a)  the  conpletlon  time  for  each  activity  not  in  the  cluster  is 
equal  to  its  upper  point,  and 

(b)  the  completion  times  for  the  activities  in  the  duster  sre 

n 

at  each  of  the  2 c possible  combinations  of  their  upper  and 
lower  percentiles. 

Let  t , p , and  I (t  ) be  as  before  and  define 
v rv  t v 

n 

2 c 

T (C)  - E p tr, 
r v-1  v v 

T*(9,  X)  - min  T+(C) 
r C r 

n 

2 c 

F (t;  C)  - E plt<tv>' 
v-1 

and 


F~(t;  9,  X)  - max  F~(t;  C) 
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Then  an  argument  completely  analogous  to  that  used  to  prove  Theorem  4 
and  Theorem  7 leads  to  the  following  theorems: 

Theorem  8: 

(a)  T*(9,  X)  is  a non-increasing  function  of  6 for  any  fixed  values 
of  X and  r; 

(b)  there  exists  a finite  value  8*  such  that  8 £ 8*  implies 

T+(e,  X)  - Tr 

for  any  X and  r;  and 

(c)  for  any  8,  X,  and  r 

Tr  £ T+(8,  X). 

Theorem  9: 

(a)  There  exists  a finite  value  8*  such  that  8 £ 8*  implies 

F~(t;  8,  X)  - F(t) 

for  every  t and  any  X;  and 

(b)  for  any  8.  X,  and  t 

F~(t;  8,  X)  < F(t). 

1.2.3  Summary 

Theorems  2-9  together  imply  that  for  any  value  of  8 and  X chosen 
by  the  algorithm  user: 

(a)  T~(8,  X)  < T £ T*(8,  X),  for  any  positive  integer  r;  and 

(b)  F~(t;  8,  X)  £ F(t)  £ F+  (t;  8,  X)  for  any  t. 


#*!'■ 
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They  also  imply  Chat,  for  any  r,  X,  and  t. 


and 


Tr(0,  X)  - Tr(e,  X) 


F (t;  0,  X)  - F-(t;  0,  X) 


would  decrease  monotonically  to  zero  if  6 were  increased.  In  fact. 
Theorems  2-9  imply  that  there  exists  a value  8*  which  doesn't  depend 
on  r,  X,  or  t such  that  0^0*  implies  that 


and 


Tr(0,  X)  «=  Tr  - Tr(0,  X) 


F"(t;  0,  X)  - F(t)  - F (t;  0,  X). 


A reasonable  approximation  for  F(t)  is 


F(t;  8,  X)  - >j[F"(t;  0,  X)  + F (t;  0,  X)]. 


2.  Comparison  with  the  Original  Subnetwork  Analysis  Procedure 


The  original  subnetwork  analysis  procedure  documented  In  Technical 

Report  No.  48  formed  the  clusters  in  essentially  the  same  way  as  the  new 

subnetwork  analysis  procedure  described  in  this  report  except  that  in  the 

original  procedure  P always  equalled  ‘i  and  0 and  X multiplied  the  point 

__________ 

difference,  u^  - t^,  instead  of  the  standard  deviation,  /Pl^  + Qu^  - m^. 

In  the  original  procedure  n,j  ■ Z nc  denotes  the  number  of  activities 

in  the  union  of  the  clusters,  and  T+(0 , X)  is  defined  to  be  the  average 

nU 

of  the  r-th  power  of  the  2 critical  path  times  when 

(a)  the  completion  time  for  each  activity  not  in  the  union  of  the 
clusters  is  equal  to  its  upper  point,  and 


(b)  the  completion  times  for  the  activities  in  the  union  of  the 

nu 

clusters  are  at  each  of  the  2 possible  combinations  of  their 
upper  and  lower  points . 

Correspondingly,  F (t;  0,  A)  was  defined  to  be  the  proportion  of  these 

nU 

2 critical  path  times  that  were  less  than  or  equal  to  t.  Analogously, 

+ nU 

F (t;  0,  A)  was  the  proportion  of  the  2 critical  path  times  less  than 
or  equal  to  t when 

(a)  the  completion  time  for  each  activity  not  in  the  union  of  the 
clusters  is  equal  to  its  lower  point,  and 

(b)  the  completion  times  for  the  activities  in  the  union  of  the 

nl) 

clusters  are  at  each  of  the  2 possible  combinations  of 
their  upper  and  lower  points. 

nu 

The  only  problem  with  this  procedure  is  that  2 may  be  quite  large  even 

for  relatively  small  values  of  (0,  A).  For  example,  if  the  original 

critical  path  contains  10  activities  and  each  critical  activity  has  one 
nU  20 

associate,  then  2 *2  - 1,048,576,  and  the  determination  of 

T*(0,  A),  F (t;  0,  A),  and  F+(t;  0,  A)  requires  the  evaluation  of  over 
2 million  critical  path  times.  On  the  other  hand,  in  the  new  procedure  the 
determination  of  T~(0,  A),  T+(0,  A),  F (t;  0,  A),  F+(t;  0,  A)  requires 

nc 

the  evaluation  of  only  2 £ 2 critical  path  times,  80  in  the  example. 

c 

Thus  the  new  procedure  greatly  reduces  the  computational  effort  required 

to  bound  the  project  completion  time  moments  and  distribution. 

nU 

A practical  alternative  to  evaluating  sll  2 critical  path  times 

nU 

called  for  in  the  original  procedure  is  to  randomly  sample  the  2 
critical  path  times  when  ny  is  large  and  base  the  bounds  on  the  sample 
critical  path  times.  A coogmter  program  for  the  original  procedure  with 
a sampling  option  is  documented  in  Appendix  C.  (This  program  supercedes 


! 

\ 
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Che  subnetwork  analysis  program  given  in  Technical  Report  No.  48.)  It 
should  be  noted  that  the  original  subnetwork  analysis  procedure  with 
the  sampling  option  in  effect  is  still  superior  to  a simple  Monte  Carlo 
simulation  of  the  subnetwork  since  the  subnetwork  analysis  procedure 
(i)  provides  the  information  in  terms  of  associates  and 

eliminants  about  which  activities  play  an  important  role 
in  determining  the  subnetwork's  completion  time  and  the 
interactions  among  activities,  and 
(ii)  samples  only  those  activities  which  most  Influence  the 
subnetwork  completion  time. 

Furthermore,  the  loss  In  accuracy  due  to  sampling  in  the  subnetwork 
analysis  procedure  seems  to  be  quite  minimal  even  for  relatively  small 
sample  sizes  - see  for  example  Table  2 and  Table  3. 

All  specific  computational  results  documented  in  this  report  are 
for  the  project  network  given  in  Figure  1.  (This  network  is  the 
simplified  network  from  a large  naval  PERT  problem.)  A listing  of  each 
activity's  two-point  approximation  except  for  its  P and  Q is  given  in 
Table  1. 

Sampling  can  also  be  used  in  the  new  subnetwork  analysis  procedure 

n 

c 

in  the  somewhat  unlikely  event  that  6 and  X are  chosen  so  large  that  2 

for  some  cluster  C is  too  large.  In  this  case  the  probability  that  the 

n 

v-th  configuration  of  activity  completion  times  v - 1,  ....  2 c is 
selected  on  a sampling  trial  is 

n 

p - nc  [P,(l  - 6 .)  + Q.6  ,] 

rv  i v,i  Hi  v,i 

where 


''r  ^ P1**  ^ -^W.  Wf  mrimaf*r*+i*  •***?• 


Table  1.  The  Upper  and  Lower  Points  in  the  Two-Point 
Approximations  to  the  Activity  Completion 
Time  Distributions  for  the  Project  Network 
in  Figure  1 


Activity 

Origin 

Terminal 

Lower 

Upper 

Number 

Node 

Node 

Point 

Point 

1 

1 

2 

0.0 

0.0 

2 

8 

12 

336.27 

429.47 

3 

2 

3 

57.47 

89.96 

4 

2 

4 

57.47 

89.96 

5 

2 

5 

57.47 

89.96 

6 

2 

6 

57.47 

89.96 

7 

2 

7 

57.47 

89.96 

8 

3 

8 

68.96 

107.95 

9 

4 

8 

68.96 

107.95 

10 

5 

8 

68.96 

107.95 

11 

6 

8 

68.96 

107.95 

12 

7 

8 

68.96 

107.95 

13 

2 

8 

150.36 

193.49 

14 

6 

10 

333.96 

403.85 

15 

3 

9 

333.96 

403.85 

16 

7 

11 

355.10 

409.85 

17 

11 

18 

141.75 

221.90 

18 

10 

13 

672.36 

783.11 

19 

9 

14 

560.89 

660.00 

20 

9 

15 

560.89 

660.00 

21 

9 

16 

560.89 

660.00 

22 

11 

17 

542.80 

638.71 

23 

12 

18 

111.10 

173.92 

24 

18 

19 

256.03 

346.98 

25 

12 

19 

302.80 

400.67 

26 

12 

20 

311.95 

410.71 

27 

11 

21 

423.58 

530.74 

28 

12 

22 

315.35 

415.71 

29 

19 

20 

7.66 

11.99 

30 

20 

21 

16.77 

22.55 

31 

21 

22 

11.49 

17.99 

32 

22 

23 

39.54 

48.87 

33 

12 

26 

301.91 

400.86 

34 

5 

26 

767.09 

892.74 

35 

12 

27 

350.31 

460.06 

36 

12 

24 

382.32 

464.98 

37 

12 

25 

385.28 

461.54 

38 

25 

24 

11.49 

17.99 

39 

24 

23 

16.28 

23.00 

40 

26 

27 

7.66 

11.99 

41 

27 

25 

20.86 

28.29 

42 

4 

28 

810.17 

976.10 

43 

23 

29 

15.32 

23.99 

j 

» 

i 
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I 


f 

I 


Activity 

Origin 

Terminal 

Lower 

Upper 

Number 

Node 

Node 

Point 

Point 

44 

28 

30 

15.32 

23.99 

45 

14 

31 

57.47 

89.96 

46 

16 

31 

49.81 

77.96 

47 

15 

31 

53.64 

83.96 

48 

17 

31 

88.12 

137.94 

49 

31 

32 

3.83 

6.00 

50 

32 

13 

49.81 

77.96 

51 

13 

33 

109.01 

152.63 

52 

11 

32 

745.50 

811.32 

53 

10 

32 

714.74 

799.83 

54 

30 

14 

0.0 

0.0 

55 

30 

16 

0.0 

0.0 

56 

30 

15 

0.0 

0.0 

57 

29 

17 

0.0 

0.0 

58 

29 

28 

0.0 

0.0 

20 


Table  2.  The  Effect  of  Sampling  in  the  Original  Subnetwork  Analysis 

Procedure:  Percentiles  of  the  Project  Completion  Time  Distribution 


Percentiles 

1 Sample  Sices 

(e.i) 

- (0,0) 

(9 

.*)  “ 

(.25,0) 

2"U-26 

25  2* 

2”“  - 2“ 

1000 

500 

200 

100 

.05 

1360 

1363  1377 

1381 

1381 

1381 

1384 

1384 

.10 

1383 

1390  1413 

1409 

1406 

1411 

1406 

1414 

.15 

1403 

1413  1423 

1424 

1421 

1429 

1424 

1423 

.20 

1416 

1423  1436 

1441 

1436 

1443 

1441 

1446 

.25 

1430 

1430  1446 

1451 

1448 

1451 

1451 

1453 

.30 

1446 

1453  1466 

1458 

1456 

1461 

1458 

1463 

.35 

1459 

1479  1472 

1468 

1466 

1471 

1468 

1471 

.40 

1476 

1486  1476 

1478 

1473 

1478 

1476 

1478 

.45 

1499 

1499  1512 

1483 

1483 

1483 

1483 

1483 

.50 

1542 

1542  1516 

1488 

1486 

1491 

1491 

1488 

.55 

1578 

1578  1575 

1496 

1496 

1498 

1498 

1498 

.60 

1578 

1578  1575 

1501 

1501 

1505 

1501 

1503 

.65 

1602 

1602  1579 

1513 

1510 

1513 

1513 

1513 

.70 

1605 

1605  1602 

1515 

1515 

1515 

1523 

1515 

.75 

1605 

1618  1605 

1525 

1525 

1525 

1528 

1525 

.80 

1621 

1618  1605 

1533 

1533 

1530 

1538 

1530 

.85 

1622 

1622  1618 

1543 

1543 

1543 

1545 

1543 

.90 

1645 

1622  1622 

1558 

1558 

1558 

1570 

1577 

.95 

1648 

1645  1648 

1585 

1585 

1582 

1587 

1597 

.975 

1648 

1648  1648 

1607 

1607 

1605 

1605 

1622 

.99 

1648 

1648  1648 

1625 

1625 

1625 

1622 

1647 

1.00 

1648 

1648  1648 

1649 

1649 

1649 

1649 

1649 

Table  3.  The  Effect  of  Skewed  Activity  Completion  Time 

Distributions  on  the  Subnetwork  Analysis  Procedures. 


Estimated  Percentiles  of  the  Project  Completion 
Time  Distribution 

Original  New 

Subnetwork  Subnetwork 

Analysis  Analysis 

Procedure  Procedure 


Monte 

(Sample  Size 

(Sample  Size 

Percentile 

Carlo 

■ 1000/cluster) 

- 1000/cluster) 

1st  Run 

2nd  Run 

1st  Run 

2nd  Run 

.05 

1405 

1402 

1396 

1370 

1374 

.10 

1423 

1416 

1416 

1390 

1398 

.15 

1435 

1430 

1433 

1402 

1410 

.20 

1447 

1438 

1443 

1425 

1437 

.25 

14j7 

1448 

1449 

1437 

1449 

.30 

1467 

1453 

1456 

1450 

1453 

.35 

1477 

1462 

1466 

1461 

1461 

.40 

1486 

1469 

1473 

1484 

1500 

.45 

1496 

1476 

1479 

1504 

1524 

.50 

1507 

1481 

1483 

1528 

1528 

.55 

1516 

1487 

1489 

1528 

1528 

.60 

1527 

1495 

1496 

1528 

1528 

.65 

1538 

1500 

1499 

1563 

1560 

.70 

1549 

1510 

1512 

1567 

1567 

.75 

1560 

1517 

1516 

1579 

1579 

.80 

1580 

1526 

1529 

1587 

1587 

.85 

1603 

1541 

1542 

1598 

1602 

.90 

1627 

1550 

1552 

1610 

1614 

.95 

1675 

1575 

1582 

1622 

1673 

1.00 

2019 

1646 

1648 

1720 

1720 

The  "1st  run"  and  "2nd  run"  correspond  to  two  different  samples  with 
different  initializations  of  the  random  number  generator. 
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6 . - 1 if  the  time  for  the  i-th  activity  in  C is  u.  in 

v,i  ,,  i 

the  v-th  configuration 

- 0 if  the  time  for  the  i-th  activity  in  C is  4^  in 
the  v-th  configuration. 

Then  with  I^t^)  as  before  the  estimated  bounds  from  C are 


T“(C) 


c 2 


E w p t 
v»l  ^VV 


and 


n 

n 2 c 

F+(t;  C)  - ^g-  Z wvPvIt(tv> 

V»1 


where  is  the  number  of  times  the  v-th  configuration  appears  in  the 
sanple,  N is  the  sample  size,  and  t^  is  the  critical  path  time  correspond- 
ing to  the  v-th  configuration  of  activity  completion  times  with  the 
completion  time  for  each  activity  A not  in  C equal  to  4^.  Similar 
modifications  are  made  for  T*(C)  and  F (t;  C) . The  corresponding 
estimators  Tr(e,  X),  T*(e,  X),  F (t;  0,  X),  and  F+(t;  9,  X)  are  all 
unbiased. 

A computer  implementation  of  the  new  subnetwork  analysis  procedure 

including  the  sampling  option  is  documented  in  Appendix  B.  In  addition 

to  reducing  the  conputational  effort,  the  new  subnetwork  analysis  procedure 

allows  the  user  to  specify  any  probabilities  (F,  Q)  for  (£,  u)  Instead 
* 

of  requiring  (1/2,  1/2).  This  would  not  be  of  great  significance  if  all 
activity  completion  time  distributions  were  symaetric.  However,  since 
many  completion  time  distributions  are  skewed,  the  ability  to  specify 
(P,  Q)  can  be  a real  advantage.  To  exemplify  this  advantage,  the 
completion  time  for  each  activity  in  Figure  1 was  taken  to  be  a linearly 
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transformed  chi-square  random  variable,  (y^  - ci^c2*  wll«re  80(1  c2 
were  determined  so  that  the  points  in  Table  1 corresponded  to  the  15-th 
and  85-th  percentiles  respectively.  This  made  the  activity's  completion 
time  distribution  highly  skewed.  Then  the  corresponding  project 
conq>letion  time  distribution  was  approximated  using 

(i)  a Monte  Carlo  simulation  of  size  1000, 

(ii)  the  new  subnetwork  analysis  procedure  with  (P,  Q)  * (1/2,  1/2) 
for  all  activities,  and 

(iii)  the  new  procedure  with  each  activity's  (P,  Q)  chosen  so  that 
the  mode  and  first  two  moments  of  its  two-point  approximation 
equalled  the  mode  and  first  two  moments  of  its  transformed 
chi-square  distribution. 

The  results  are  given  In  Tabla  3.  If  the  Monte  Carlo  approximation  is 
used  as  a basis  for  comparison,  the  value  of  being  able  to  specify  (P,  Q) 
is  obvious. 

The  Monte  Carlo  PERT  simulation  program  used  in  the  above  experiment 
is  documented  in  Appendix  D and  supercedes  the  Monte  Carlo  PERT  program 
given  in  Technical  Report  No.  48. 

In  the  special  case  where  every  activity's  (P,  Q)  is  (1/2,  1/2), 
the  T^(C)  and  F (t;  C)  can  be  computed  on  the  basis  of  the  2 c critical 
path  times  corresponding  to 

(a)  the  completion  time  for  each  activity  A not  in  C being  equal 
to  the  mean,  m^,  and 

(b)  the  completion  times  for  the  activities  in  C being  at  each  of 

n 

the  2 c possible  combinations  of  their  upper  and  lower  points. 
This  is  a change  from  the  usual  computation  in  that  the  activities  not 
in  C are  at  their  means  here  instead  of  their  lower  points.  This  change 
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will  tend  to  improve  the  estimators  Tr(0,  X)  and  F+(t;  6,  X).  However, 
this  change  is  only  guaranteed  not  to  invalidate  Theorems  4 and  7 when 
all  (P,  Q)  are  (1/2,  1/2).  The  computer  implementation  of  the  new 
subnetwork  analysis  procedure  includes  the  option  to  make  this  change. 

3.  Conclusion 

The  new  subnetwork  analysis  procedure 

(i)  forms  associates,  eliminants,  and  clusters  in  essentially 
the  same  way  as  the  original  procedure, 

(ii)  bounds  the  project's  completion  time  moments  and  distribution 
primarily  on  the  basis  of  the  individual  clusters  Instead 
of  on  a pooled  cluster,  and 

(iii)  approximates  an  activity's  completion  time  distribution  by 

a two-point  distribution  with  possibly  unequal  probabilities 
for  the  two  points  instead  of  always  equal  probabilities. 

The  advantage  of  (li)  is  that  much  much  fewer  critical  path  times  need 
to  be  evaluated  in  determining  the  bounds  on  the  project  completion 
time.  The  advantage  to  (iii)  is  the  ability  to  better  approximate 
skewed  activity  completion  time  distributions. 

Computer  Implementations  of  the  new  subnetwork  analysis  procedure, 
the  original  subnetwork  analysis  procedure,  and  a Monte  Carlo  simulation 
algorithm  are  documented  in  Appendices  B,  C,  and  D respectively.  Both 
subnetwork  analysis  procedures  include  options  to  use  sampling  for 
large  clusters. 
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APPENDIX  A 
Proof  of  Theorem  1 

The  principle  objective  of  this  appendix  is  to  prove  the  follow- 
ing theorem: 

Theorem  1:  In  any  acyclic  network  with  no  cut  vertices  there  is  at 

least  one  bridge  for  any  pair  of  consecutive  arcs. 

The  networks  considered  in  this  appendix  are  assumed  to  be 

t 

acyclic,  have  no  cut  vertices,  and  have  one  source  and  one  sink.  Also 
the  two  arcs  and  A^  are  any  two  adjacent  (consecutive)  arcs  with  A 
preceding  Ay 

Definition  1:  A bridge  over  A^  and  A^  is  any  arc,  say  Ay  such  that 

all  paths  from  the  source  to  the  sink  passing  through  A^  do  not  pass 
through  either  A^  or  Ay 

Definition  2:  An  origin  violator  of  A^  and  Aj  is  an  arc,  say  Ay 

such  that  there  exists  a path  from  the  terminal  node  of  A^  to  the  sink 
which  passes  through  Ay 

I 

Definition  3.  A terminal  violator  of  A^  and  A ^ is  an  arc,  say  Ay 
such  that  there  exists  a path  from  the  source  to  the  terminal  node  of 

IA^  which  passes  through  Ay 

An  intuitive  feeling  for  these  definitions  can  be  obtained  by  con- 
sidering any  path  P*  from  the  source  to  the  sink  which  passes  through 

| 

A,  and  A2.  If  an  arc  A^  is  an  origin  violator,  then  there  is  a path 
from  the  source  to  the  sink  which  follows  along  P*  through  the  terminal 
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node  of  Ax  and  then  goes  through  Ar  This  path  "originates"  from  P* 
too  late  for  A3  to  be  a bridge  over  Ax  and  A 2*  Similarly,  if  A3  is 
a terminal  violator,  then  there  is  a path  from  the  source  to  the  sink 
which  passes  through  A3  and  then  joins  into  P*  before  P*  passes  through 
the  terminal  node  of  A^  This  path  "terminates"  into  P*  too  early  for 
A3  to  be  a bridge. 


BRIDGES  over  A^  and  A2:  H,  I,  N,  0,  P,  and  Q 

ORIGIN  VIOLATORS  of  A3  and  A2:  A2,  F,  G,  J,  K,  R,  and  S 

TERMINAL  VIOLATORS  of  A^  and  A2:  A^,  D,  E,  L,  and  M 

The  following  three  lemmas  are  straightforward  consequences  of  the 
definitions  of  a bridge,  an  origin  violator,  and  a terminal  violator. 


Lemma  1:  Every  branch  in  Che  network  is  either  a bridge  over  A^  and 

an  origin  violator  of  Aj  and  A2>  or  a terminal  violator  of  A^ 
and  A^. 

Proof  of  Lemma  1:  Suppose  that  A^  is  not  a bridge.  Then  there  exists 

a path  P from  the  source  to  the  sink  which  contains  A^  and  either  A^ 
or  A2> 

Suppose  that  P contains  A^.  If  A^  - A^  or  A^  precedes  A^  on  P, 
then  P contains  a path  from  the  source  to  the  terminal  node  of  A^  which 
passes  through  A^,  and  A^  would  be  a terminal  violator.  On  the  other 
hand,  if  A^  follows  A^  on  P,  then  P contains  a path  from  the  terminal 
node  of  A^  to  the  sink  which  passes  through  A^,  and  A^  would  be  an 
origin  violator. 

Suppose  that  P contains  If  A^  “ A2  or  A^  comes  after  A2  on 

P,  then  P contains  a path  from  the  terminal  node  of  A^  (the  origin 
node  of  A 2)  to  the  sink  which  passes  through  Aj,  and  A^  would  be  an 
origin  violator.  If  A^  comes  before  A2  on  P,  then  P contains  a path 
from  the  source  to  the  terminal  node  of  A^  (the  origin  node  of  A2> 
which  passes  through  A^,  and  A^  would  be  a terminal  violator.  This 
completes  the  proof  of  Leoma  1. 


Lemma  2:  A^  is  a terminal  violator  of  A^  and  Aj,  and  A2  is  an  origin 
violator  of  A^  and  Aj. 

Lemma  3:  Any  arc  A^  cannot  be  both  an  origin  violator  of  A^  and  A2 


and  a terminal  violator  of  A^  and  Aj. 
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Proof  of  Lemma  3:  Suppose  that  an  arc  is  both  an  origin  violator 

and  a terminal  violator.  Since  A^  is  an  origin  violator,  there  exists 
a path  from  the  terminal  node  of  A^  to  the  origin  node  of  A^.  Since 
A^  is  a terminal  violator,  there  exists  a path  from  the  terminal  node 
of  A^  to  the  terminal  node  of  A^.  The  existence  of  these  two  paths, 
however,  Implies  the  existence  of  a circuit  which  contradicts  the  given 
acyclic  structure  of  the  network.  This  completes  the  proof  of  Lemma  3. 

Proof  of  Theorem  1:  Since  the  terminal  node  of  A^  cannot  be  a cut 

vertex,  there  exists  a path  P from  the  source  to  the  sink  which  does 
not  pass  through  the  terminal  node  of  A^.  Denote  the  arcs  on  P by 
C^,  C2 0^  with  Ci_1  preceding  on  P. 

Suppose  that  none  of  C^,  C^,  . ...  Cp  are  bridges  over  A^  and  Aj. 
Then  Lemma  1 and  Lemma  3 together  imply  that  each  of  , C2 , . . . , Cp 
is  either  an  origin  violator  of  A^  and  A2  or  a terminal  violator  of 
Aj  and  A2  but  not  both.  Since  the  origin  node  of  is  the  source, 
cannot  be  an  origin  violator  and  must  be  a terminal  violator. 
Similarly,  since  the  terminal  node  of  is  the  sink,  Cp  cannot  be  a 
terminal  violator  and  must  be  an  origin  violator.  Hence,  there  exists 
j _>  1 such  that  C^,  C2,  ...,  Cj  are  all  terminal  violators  and  is 

an  origin  violator. 

Since  Cj  is  a terminal  violator,  there  exists  a path  from  the 
terminal  node  of  to  the  terminal  node  of  A^.  Furthermore,  since 

is  an  origin  violator,  there  is  a path  from  the  terminal  node 
of  A^  to  the  origin  node  of  C (the  terminal  node  of  Cj ) . These 
two  paths  imply  the  existence  of  a circuit  from  the  terminal  node  of 
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to  the  terminal  node  of  and  then  back  to  the  terminal  node  of 
Cj . (The  definition  of  P implies  that  the  terminal  node  of  is  not 
the  terminal  node  of  A^.)  This  contradicts  the  given  acyclic  structure 


of  the  network  and  completes  the  proof  of  Theorem  1. 


APPENDIX  B 


New  Subnetwork  Analysis  Program 


The  New  Subnetwork  Analysis  Program  is  an  implementation  of  the 
analytical  procedure  described  in  Section  1 of  this  report.  The  basic 
required  input  is 

(a)  an  acyclic  network  with  one  source  and  one  sink, 

(b)  two  points  from  each  component  activity’s  completion 
time  distribution 

(c)  probability  P to  be  associated  with  the  lower  point, 
and 

(d)  specified  values  for  the  algorithm  parameters  0 and  A. 

The  output  is  mainly 

(a)  upper  and  lower  bounds  on  the  moments  of  the  network 

completion  time,  T*(0,  A)  and  T^O,  A)  r » 1,  2 10; 

(b)  upper  and  lower  bounds  on  the  distribution  function  of 
the  network  completion  time,  F+(»;  0,  A)  and  F (•;  0,  A); 
and 

(c)  an  approximate  network  completion  time  distribution, 

F(-J  0,  A)  - |[F+(-;  0,  A)  +f"(-;  0,  A)]. 

The  basic  computational  technique  for  determining  critical  path 
times  is  the  Simplex  Algorithm.  This  algorithm  is  applied  to  the  dual 
problem.  The  Simplex  Algorithm  is  used  instead  of  the  standard  network 
analysis  techniques  because  the  Simplex  Algorithm  is  ideally  suited  for 
the  type  of  parametric  programing  required  to  evaluate  several  critical 
path  times  when  only  the  activity  times  vary  from  one  problem  to  the  next. 

A listing  of  the  Subnetwork  Analysis  Program  and  a program  flowchart 
are  given  at  the  end  of  this  appendix. 


Specific  Input  Instructions: 

Card  1.  Col.  1-3:  The  number  of  activities  in  the  network,  Format  (13), 


Col.  4-6:  The  number  of  nodes  in  the  network,  Format  (13). 

For  each  activity  one  card  with: 

Col.  11-15:  The  origin  node  of  the  activity,  Format  (15). 

Col.  21-25:  The  terminal  node  of  the  activity,  Format  (15). 

Col.  31-40:  The  lower  point  on  the  activity's  completion  time 

distribution,  Format  (F10.0). 

Col.  41-50:  The  upper  point  on  the  activity's  completion  time 

distribution,  Format  (F10.0). 

Col.  51-60:  The  probability  P to  be  associated  with  the 

lower  point  (1-P  will  be  associated  with  the 
upper  point),  Format  (F10.5). 

Next  Card.  Col.  1:  OPTON1.  0PT0N1-1  Implies  that  the  program  will 

terminate  after  the  clusters  have  been  formed  on 
the  basis  of  associates  and  ellminants.  OPTON^l 
implies  that  the  program  will  follow  the  normal 
procedure . 

Col.  2:  0PTON2.  0PT0N2-1  implies  that  the  lower  bounds 

on  the  moments  of  the  project  completibn  time  and 
the  upper  bound  on  its  distribution  will  be 
determined  by  using  all  activity  times  outside  the 
cluster  at  their  means  instead  of  their  lower 
points.  This  is  only  guaranteed  to  be  a valid 
procedure  when  all  (P,  Q)  - (1/2,  1/2).  0PT0N2»*1 

Implies  that  the  program  will  follow  the  normal 
procedure . 
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1 -r, . 


Next  Card.  Col.  1-3:  IEDF,  The  program  computes  an  absolute  upper  and 

lower  bound  for  the  network  completion  time.  This 
range  is  subdivided  into  IEDF  equal  parts  and  the 
approximate  distribution  function  (F+,  F , F) 
values  are  printed  at  each  of  these  dividing  points. 
IEDF  would  usually  be  between  10  and  100.  IEDF, 
Format  (13). 

Next  Card.  Col.  1-5:  0,  Format  (F5.2). 

Col. 6-10:  X,  Format  (F5.2). 

Next  Card.  Col. 1-10:  SAMSIZ.  The  number  of  activity  time 

configurations  to  be  randomly  selected  for 

explicit  consideration  in  each  cluster  analysis. 

n 

If  SAMSIZ  < 0 or  SAMSIZ  > 2 c,  all  activity  time 
configurations  will  be  explicitly  considered  - no 
random  sampling  will  be  done.  Format  (110). 

The  nodes  should  be  numbered  1,  2,  . ...  n with  the  source  being  number 
1,  the  sink  being  number  n,  and  the  other  node  numbers  being  arbitrary. 

The  activities  should  be  numbered  1,  2,  ...  in  any  order  desired. 

Current  Dimension  Restrictions : 

Currently  the  program  is  dimensioned  for  a maximum  of 
60  Activities 
40  Nodes 
25  Clusters 

25  Activities/Cluster  and  IEDF  <_  100. 

Example: 

The  Program's  input  and  output  are  illustrated  in  terms  of  the 


network  in  Figure  B-l. 
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NF*  SUBNETWORK  ANALYSIS  PROGRAM 
IMPLICIT  REAL*8  (A-H.O-Z) 

FOR  THE  SAKE  OF  IDENTIFYING  THE  APPPOPRIATE  DlMFNSIONSt  LET 
M = THE  NUMBER  OF  ACTIVITIES  IN  THF  NETWORK 
nmm  = number  of  nodes  in  the  network 

NMMPl  = NMM  ♦ I 
N - M ♦ NMM 

L = THE  LENGTH  OF  THF  CRITICAL  PATH 
C.  = THE  MAXIMUM  NUMBER  OF  BRANCHES  IN  A CLUSTER 
I E OF  = THE  NUMBER  OF  DIVISIONS  IN  THE  FMPI R IC AL 
DISTRIBUTION  FUNCTION 


INTFGFR  TAIL!  M ) , HF ADI  Ml.ASSGpPJ  L.L  ) . CL  I NCL ( L .L ) , EGRP ( L > 
DIMENSION  nincli  c>,inclus(  l.  O.NCLINCC  L) 

DIMENSION  FDf IFDF) ,NLE*  D(t  . IEDF),NSAVE{  IFDFl.NIRlL) 

DIMENSION  AVGI  L).THAT(L) 

DIMENSION  LEFtIM  J.LEFTOIM  ),NONCP(M) 

DIMENSION  INhASF(NMM) ,XNODE<NMM) 

DI  MFNS I ON  XB1  (NMMPl  ) , y ] (NMMPl  > , REDCOS(N) . I ST  AT( n I 
DI MENS  I ON  ICRITO(LI,NiNAg(M).iCRITN(LF1).CTIMF(N»,  COT  IN ) 

D I MENS  ION  K 1 ( L ) < I OH  ( L I < FLO( M > , FH I ( M ) , S I GMA ( M > , B 1 I NV ( NMMPl  . NMMP 1) 
REAL  MOMFNT(L.IO) 

DIMFNSICN  PP(M),PO(M)  i 

f 

OF  COUP  36  THESE  DIMENSIONS  ArE  MERELY  UPPEP  BOUNDS 


COMMON  B 1 I NV.PFDCaS.CT  ImE.XBI  . INBASF.HE  Ad  » TA  I L » NMMP I * NMM  « N • I ST AT 
COMMON  M.MPI 

INTEGER  tAIL(6C)  ,HEAO((>C  ),ASSGRP( 25,25)  .CLINcLC  25. 25). EGRP (25) 
INTEGER  SAMSl Z.PANSAM 

01 MFNS  I ON  NTNCLC  25)  , INCLUSI 2 5,25)  , NCLlNCI  2b) 

DIMFNSICN  FQ( 1"0  ).NIB(?5  ) 

R E AL* B NLEFD(2s.l90 ) .NS AVE< 100  ) 

DIMFNSICN  AVG( 25 ) ,THAT(  25) 

DIMFNSICN  INBASE(AO) 

0 IMF  NS  I ON  XNODE ( 40 ) 

D I Mr NS  I ON  XB1 ( 4 1 ) , Y I ( 41  ) .RFDCOSI 1C0>.ISTAT(10C> 

DI MENS ICN  B1INV( A l , A 1 ) , KB<  2b) . IBBI 25) ,FLO( 60) .PHI (60  » .SI GMA ( 60) 

D I Mr  NS  ION  iCRi TP (25) .NI  NAGI 5C  ) . ICRITN(26> .CDTC 100>«CTIME(100) 

D I MFNS ICN  LEFT(60),LEFtO(  60  ) .NONCPI60) 

REAL*B  LAMBDA. MCMENT(25, 10 ) 

DIMFNSICN  PP ( 60 ) . PQ ( 6C ) 

INTFGFR  OPTONl . OPTON2 

M = THE  NUMBER  OF  ACTIVITIES  IN  THE  NETWORK 
NMM  = THF  number  oF  NODES  IN  THF  PpPT  NETWORK 
REA0(5.1C0)  M.NMM 
FORMAT (211) 

N=NMM+M 
MP1 =M+ 1 
NMMP 1 = N MM* l 

THE  ACTIVITIES  ARE  DESCRIBED  IN  TERMS  OF  THEIR  NODES 
!I=ThF  TAIL  NODE,  THE  ORIGIN  NODE 
JJ=TEF  HEAD  NODE,  THE  TFRMlNAL  NODE 
Fl.O  = The  l OWER  PUINT 
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FHI  = THE  UPPER  POINT 

SIGMA  = (FHI  - FL0)*DSORT(  pp*(l-pp))  = STD.  DEVIATION 
PP  = THE  PROBABILITY  OF  THF  LOWER  POINT 


OO  61 r I =1 , M 

READ  C 5 « 2EP 1 ) II.JJ.FLOII  I.FHIU  >,PP(I) 

PQ(  I )=1  .OC-PP( I ) 

SIGmA(  I ) = ( FHI  ( n-FLO<  I ) ) *DSOR T { PP<  I loPQC  I > ) 
01  FORMAT ( 10X. I5.5X.I5.  5x  , E 1 0 . 0 . F 1 0 . 0 , F 1 C . 5 > 
COT  ( I l=PP(  I )*FLO(  t)«-PQ<  I )*FHI(  I ) 

CTI ME ( I > = COT ( I ) 

TAIL! I > = TI 
610  HEAD!  I ) = JJ 

COT  = THF  ORIGINAL  RIGHT-HAND  SIDES.  I.E. 
OF  FLO  AND  FHI 

CtIME  = THE  CURRENT  RIGHT-HAND  SIDES 
DO  556 1 0 I = MP ] . N 


THE  MEANS 


DO  556 l 0 

5561 0 CTI  MF ( I 1 


optoni  =1  implies  that  the  program  will  terminate  after 

THE  CLUSTERS  HAVE  BEEN  FORMED.  NO  BOUNDS  ON 
THE  PROJECT  COMPLETION  TIME  MCMENTS  OR 
DISTRIBUTION  WILL  BE  determined. 

CPTON1  NTT  a 1 IMPLIES  THAT  THF  NORMAL  PROCEDURE  WILL  BF 
followed . 

OPTONs  =1  IMPLIES  THAT  THE  LOWER  BOUNDS  ON  THE  MOMENTS 
ANO  THE  UPPpP  ROUND  ON  THE  DISTRIBUTION  WILL 
BE  L'FTERMINFO  USING  ALL  ACTIVITY  TIMES  OUTSIDE 
THF  Cl  USTFR  AT  THEIR  MFAN.  THIS  PROCEDURE  IS 
ONLY  GUARANTFED  TO  BE  VALID  WHEN  ALL 
(P.O)  = < .5,  .51  . 

OPTON2  NOT  = 1 IMPLIES  THAT  THE  NORMAL  PROCEDURE  WILL  BE 
FOLLOWED . 


QPTON2 


OPTON2 


RFAOI5,  77651  ) OPTON 1 . OPTCN2 
77551  FORMAT ( 1 C 1 1 ) 

I=OPTONl + QPT  ON 2 
IF(I.GE.l)  WrITF(6, 77553J 
77553  FORMAT ( 1H1 ) 

IFIOPTCNl  .50.1)  WRTTElb. 77552) 

7755?  FOPMATl  1H0 « 1 0 Y » • OPT  I ONJ  = i AND  THE  PROGRAM  WILL  TERMINATE  AFTER  THF 

* CLUSTERS  HAVE  BEEN  FOPMFD . • . / * II  X . * NO  BOUNDS  ON  THE  PROJECT  COMPL 
*E  T I ON  T I Mp  MCMENTS  OP  DISTRIBUTION  WILL  BE  DETERMINED.*) 

IFIOPTCN2  . EO . 1 ) WRITE!  6, 7755A) 

7755A  FOFMAT ( I HO  . 1 0 X . • OPT  I ON2 - 1 AND  THE  LOWER  BOUNDS  ON  THE  PROJECT  COMP 
WETION  TIME  MOMENTS  AND  THF  UPPER  ROUND  ON  THE  PROJECT  COMPLETION  T 

* I ME  m STR  IBUTI ON • ./ . I IX.  • WILL  Bf  DETERMINED  USING  ALL  ACTIVITY  TIM 

*FS  OUTSIDE  THE  CLUSTER  AT  THEIR  ME AN  . • , / . 1 1 X , • TH I S PROCEDURE  IS  ON 
*LY  GUARANTEED  Tq  PE  VALID  WHEN  ALL  (P.O)  = (.5.  .5)  .*) 

C I EOF  = THF  NUMBER  OF  DIVISIONS  IN  THE  EMPIRICAL 

C DISTRIBUTION  FUNCTION 

REAn(F.lOO)  I EOF 
WRITF(6,P7-'P  ) 

270 r*  FOFMAT(  IHl  ,15X. ‘INITIAL  INPUT*) 

WP  T TF ( 6 . 276 l ) 

2701  FORMAT ( IHO. IPX, ‘ACTIVITY  ORIGIN  TERMINAL  LOWEP  POINT  UPPER  P 
* 01  NT  MEAN  STANDARO  DEVIATION  PPOB.  LOWER  *!,') 


UPPER 


no  ? 7C  A 1 = 1. M 

WRI  TF(  6,270?)  !,TAJL(  I)  ,HCAD<  I I . FLUI  I I.FHK  I)  , CT  I ME  ( II  .SIGMA  ( I'l  «PP 


iai  afete --f  ft 


n n n n n n 


2 70? 


C 

C 

c 

r 

c 


c 

c 

c 

r 

c 

r 

c 

c 


c 

c 

r 


r 

r 


FORMAT ( lH  , 1 3X . T3.5X,  13. 7X,  13. 5X.F1C  .4.3X  .F10.4,3X,F1C.4,4X,F1 0.4, 
*!4K,Ff .4) 


60  1 C 
I 0 4 
200  1 
?002 


12 

10 


3C 


THE  FOLLOWIN'".  INDICATORS  APfc  USED! 


IPARM  = 1 


IPAuM  = 2 


IPARM  = 3 


IMPLIES  THF  CRITICAL  PATH  TIME  WHEN  ALL  ACTIVITY 
COMPLETION  times  are  set  equal  to  thfir  means  is 
OPING  DETERMINED 

ROUND 
DF  ING 
BOUND 
OF  INC, 


(PA PM  > i 


IMPLIES  THAT  THE  LOWER 
FOR  THE  SUBNETWORK  IS 
IMPLIES  that  the  UPPER 
FOR  THE  SUBNETWORK  IS 
WHEN  INOCXL=C  IMPLIES 
BEING  DETERMINED 


TIME 


ON  THE  COMPLETION 
DETEPMINFD 

ON  THE  COMPLFT I ON  TIME 
DETERMINFO 


THAT  THE  ASSOCIATFS  ARF 


INnEXL=' 
I NLi  E XL  = 1 


INpt£XL=3 


implies  that  initial  clusters  are  still  being  formed 
implies  that  the  lletovers.  their  eliminants,  and 
pooled  clusters  are  being  determined 
implies  THAT  THF  2**NINCL(  I RUNS  FOP  EACH  CLUSTER 
APE  BEING  MADE 


ICOCP  = 0 


I CROP  = 1 


IMPLIES  THAT  THE  PROCEDURE  FOR  DETERMINING  UPPER 
POUNDS  ON  THE  MOMENTS  OF  THE  NETWORK  COMPLETION  TIME 
AND  LOWER  BOUNDS  ON  THE  DISTRIBUTION  OF  COMPLETION 
TIMES  HAS  NOT  BEEN  BEGUN 

IMPlILS  THAT  THF  PROCEDURE  FOR  DFTERMINING  UPPER 
unuNrjs  on  the  moments  of  the  nftwork  completion  time 

ANO  LOWER  ROUNDS  ON  THE  DISTRIBUTION  OF  COMPLETION 
TIMES  IS  BEING  INITIALIZED 


CONTI  N(JF 
ICRCP 


IMPLltS  THAT  THE  PROCEDURE  FOR  DETERMINING  UPPFR 
BOUNDS  ON  THE  MOMFNTS  OF  THE  NFTWORK  COMPLETION  TIME 
ANO  LOWER  BOUNDS  UN  THE  DISTRIBUTION  OF  COMPLETION 
TIMFS  IS  BCING  CARRIED  CUT 


t PAPM= I 
I NDFXL=0 
I CBCP=  0 
CONTINUE 
DO  104  1=1 ,NMM 

TNPAsF ( I >=M4 l 
DO  2C0I  J=1.M 
1ST  AT  < J1=C  . 

DC  200 2 J = MP 1 . N 
I ST  ATI J 1 = 1 
DO  K 11=1 , NMMP1 
DO  1 2 L-l i KMMP1 
B 1 INVIl  .11)  = 0 . 
B1  I NV(  I I . I I ) = 1 

DO  30  1=1. NMM 

xni ( i ) = o . 

XR| (NMMP1 1 = | . 

TOLR 1 = 1 • CD- I C 


120 

121 
1 22 

123 

124 

125 

126 

127 

128 
129 
13C 

131 

132 

133 
1 34 
1 35 
I 36 
1 37 

138 

139 

140 
1 41 

142 

143 

144 

145 

146 

147 

148 

149 
1 50 
151 
I 52 
I 53 

154 

155 

I 56  | 

157 

158 
1 59 
160 
1 6 1 
162 

163 

164 

165 

166 

167 

168 

169 

170 
1 71 

172 

173 

174 
1 75 

176 

177 

178 
179, 


4v*fcjr.n 


C START  THE  SIMPLEX  ALGORITHM 

C SOLVE  THE  DUAL  PROBLEM 

C Th£  NUMBFR  OF  VARIABLES  IS  M REAL  ♦ NMM  SLACKS 

' for  A TOTAL  of  N VARIABLES 

c 

3 SO  CONTINUF 

28Qr  DO  2 3 J=1,N 
RATS  = 0. 

IF  (ISTATI J) .EQ.l)  GO  TO  52800 
IF  (J.GT.M)  GO  TO  22 

RATS  =-Bl  INVl  1 .HEAP! J ) + l >+01 INV<  I .TA IL( J> FI > ♦ CTlME(J> 

GO  TC1  S2300 

22  RATS  =-Rl  INV(  1 , J-M+j  ) 

S?800  RFOCOS I J ) = RATS 

23  CONTINUE 
22800  CONTINUE 

I PMAX= 1 

RMAX=EpOCOS< 1 ) 

00  24  J=P.N 

IE(PEOCnS(JI  .LF.  OMAX)  GO  TO  24 
RMAX  = PEOCCS<  J I 
IRMA  X=  J 

7 4 CONTINUE 

IFIRMAX  .LF.  T OL PI)  GO  TO  401 
22824  CONTINUE 

00  26  L = 1 < NMMP 1 
IF  (IPMAX.GT.MI  GO  TO  50026 

Y1(l)  =-Bl INV(L,TaIL<  IRMAX )+l)+B!  !NV<L.HEADURMAX)M) 

GO  TO  26 

50026  V 1 < L ) = HI INV  (L.  IRMAX-MFl) 

26  CONTINUE 

Yl(l)  = Y 1 ( 1 ) - CTIME(IRMAX) 

NUMBER  = 0 

00  27  L -2  »NMMP  I 

27  IFIYi(L)  .LE.  T0LR1I  NUMBFR=NUMBERF I 
I F ( NUMBER  .FO.  NMM)  GO  TO  403 
RMIN=.RRD  ?4 

1 R M I N=  0 . 

OO  32  11=2. NMMP I 

IF(V1(  II  ).LF.  TOLRl)  GO  TO  32 

RATS  = XB 1 ( I I ) / Y 1 (II) 

RR  = RAT  S-RM IN 

IF(RP  .GE.  0.00)  GO  TO  32 
RMI N=PATS 
I RM  t N= I I 

32  CONTINUE 

DO  33  J = 2»  NMMP 1 

WW  = H1 I NV<  IRMI N ,J)/YI(IRMIN  I 
DO  37  L = I . NMMP I 

37  B 1 I N V ( l . J) =31  I NV(L, J ) -W W«Y  t (L) 

33  BIINV(IRMTN  ,J)=MW 
C 

C UPDATE  THE  BASIC  VARIABLES:  INBASE  AND  XBI 

C 

ISTATI INBASEt IPMIN-1))=0 
ISTATI IRMAX) =1 
INBASF<  I PM  IN— 1 1= IRMAX 
W=XB1(IrMIN  )/YI(|RMIN  ) 

OO  38  I=J,NMM^I 
XBI  ( I ) sXBI ( I >-Yl( I l*W 


38 


n r> 


t 

s 1 


XBMIRMIN  ) = »' 

GO  TO  .150 
WRITFI6.530 ) 

FORMAT! 1h0.5X.‘ND  FEASIBLE  SOLUTION  EXISTS.  CHECK  yOuR  INPUT  DATA 
*.  • ) 

WRITE! 6.e5C ) 

FORMAT ( IH1  ) 

GO  TO  999 

END  OF  THE  SIMPLEX  ALGORITHM 
CONtINuF 

IF<  ICBCP.EO.  1 ) G(1  TO  6008 
IF(  INDExt .E0.2  ) GO  TO  3204 


KKK  = THE  NUMBER  OF  NODES  ON  THE  CRITICAL  PATH 

KB(L)=  THE  L-TH  node  IN  THE  CRITICAL  PATH.  COUNTING  BACKWAPDS 
FROM  THE  TERMINAL  NODE 

KKB=  THF  NUMRER  OF  ACTIVITIES  ON  THE  CRITICAL  PATH 
I BB(  L )■=  THE  L-Th  ACTIVITY  ON  THE  CRITICAL  PATH,  COUNTING 
BACKWAPDs  FROM  THE  TERMINAL  NTDE 

CONTINUE 


INBASE  IS  A SET  OF  M INTEGER  VARIABLES  WHICH  INDICATE  THE 
COMPOSITION  OF  THE  CURRENT  BASIS.  FOR  EXAMPLE, 

INRASE(K)  = 7 IMPLTFS  THAT  THE  K-TH  COLUMN  tN  THE  BASIS  B 
CORRESPONDS  TO  THE  7-TH  VARIABLE 


ISTAT  INDICATES  THE  BASIC  STATUS  of  EACH  VARIABLE 

ISTAT(K)  = 1 IMPLIES  THAT  THE  K-TH  VARIABLE  IS  IN  THE 
DUAL  BASIS 

ISTAT(K)  = 0 IMPLIES  THAT  THE  K-TH  VARIABLE  IS  NOT  IN  THE 
DUAL  BASIS 


the  FOLLOWING  STATEMENTS  DETERMINE  THE  NODES  AND  ACTIVITIES  ON 
THE  CRITICAL  PATH 


THE  DUAL  SOLUTION  IMPLIES  THE  FOLLOWING  OPTIMAL  SOLUTION  TO  THE 
PRIMAL  PERT  PROBLEM.  HOWEVER  SOME  OF  THE  NODE  Tl MESI OTHER  THAN 
THE  LAST  ONE)  MAY  BE  HIGHER  THAN  NECESSARY.  THUS  IN 
DETERMINING  THE  CRITICAL  PATH  AN  ALTERNATIVE  OPTIMAL  SOLUTION 
MAY  HA  VF  TO  BE  IDENTIFIFD. 

R1INV  IS  NOT  CHANGED. 

DO  R3C02  1 = 1 « N MM 

XNCDE ( II=»IINV(1.I*1) 

KKK=  1 

KOI  1 > = NMM 

IK=KR(KKK) 


DETERMINE  WHETHER  THF  TIME  TO  REACH  NODE  I K 

AS  LARGE  AS  INDICATED  FROM  THE  DUAL  SOLUTION 


IS  NECESSARILY 


SMI N*9 99999 . 

I SM!N=9 

DO  B1C00  1=1 .M 


loco;  l i:  v v 
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IFIHFAUI 1 ) ,NF . IKl  GO  TO  83000 

81!  AC.«=XNnr>E<HF  AO  ( I ) )-XNODElTAII  ; t ) )-CT  I ME  ( It 
IFISLACK.GE .SMIN)  GO  TO  83000 
SMI N=SL AC< 

ISMI N= I 

83000  CONTINUE 

IF< SMI N .LT.O .00 0 1 ) GO  TO  83003 

THF  TIME  FOR  NODE  IK  WAS  UNNECCESSAR1LY  large 

XNPOF ( lKl  = XNOOF(  IKl-SMIN 
KKK=KKK-1 
GO  TO  33001 
3003  IBRIKKK »=ISMIN 
KKK=KKK+1 

KP(KKK) =TAlL( ISMIN) 

I F < TAIL ( ISMIN ) .GT. I > GO  TO  83001 
KKR^KKK-1 

if ( INDEXL.cO. 1 ) GO  TO  3121 
IPAPM= I PAPM+  1 

IF  ( IPAPM.GT.4  ) GO  TO  20  1 C 
IF<  IPAR v.EQ. 3 1 GO  TO  6A00 
IF ( IPaRM.FO.4 ) GO  TO  6401 

ICP!TP(L>=  the  L-TH  ACTIVITY  ON  THF  ORIGINAL  CRITICAL  PATH 
KCOR-  THE  NUMBER  OF  ACTIVITIES  ON  THF  ORIGINAL  CRITICAL  PATH 

TOTAL  = R I I N V ( 1 , KMMP1  I 
KfP  6 = K K H 
ICPITNI 1 1 = NMM 

on  2802  I=1,KCPR 
ICPITNI I+l ) = KE ( I ♦ 1 » 

2802  ICPITPI  I ) = IRH<  I ) 

X = T OT  AL 

WRI TFI 6,851 ) X 

851  FCpMAT(  IHO  ,5X  , • Th^  CRjTICAL  PATH  TIME  WHEN  EACH  ACTIVITY*‘S  COMPLE 
*TlON  TIME  IS  SET  EQUAL  TO  ITS  MF AN  IS  = ‘,015,5) 

WPl TF( 6 ,7606 1 KKK 

7606  FORMAT!  IHO . 1 OX ,« THP  ‘,l3,»  NODFS  ON  THF  CRITICAL  PATH  ARE  AS  FOLLO 
*w ? RFGINNING  WITH  THE  TERMINAL  NODE:') 

WRITEI6 ,7707  1 < KRI I ) , 1=  I ,KKK  ) 

7707  FORMAT ( 15X, 201  I 3,‘ ,•  )) 

WRI TF 16,7710)  KKO 

77l 0 FORMAT! IHq.Iox , ‘THE  ‘,!3,»  CRITICAL  ACTIVITIES  ARE  AS  FOLLOWS  REGf 
INNING  w I TH  THE  TFRMJNAL  ACTIVITY!*) 

WPI TFI 6. 7707  ) I IBB  I I ) , I = 1 , KKR  ) 

RFADI5,2<520  ) THFTA. LAMBDA 
2920  FORMAT  I 2 FE  • 2 ) 

WPITE(6,3C71  ) THFTA, LAMBDA 

3071  FPFMAT  I IHO ,1  OX , • THFTA  = ‘,E|5.5.‘  LAMBDA  a »,EI5,5) 

C 

C 4 AMS  I Z = THE  NUMUEp  OF  ACTIVITY  time  CONFIGURATIONS  il  BE 

r RANDOMLY  SELFCTEO  for  CONSIDERATION  IN  EACH  CLUSTER 

r 

r note:  since  this  is  a random  sample  , SOME  PERCENTILF 

c COMBINATIONS  MAY  BE  CONSIDERED  MORF  THAN  ONCE, 

C 

READ  15,3239)  SAMSIZ 
3209  FORMAT  1110) 


300 

1 

30  1 

I 

302 

1 

303 

■ 

304 

■ 

305 

306 

1 

307 

308 

1 

309 

I 

310 

1 

311 

I, 

112 

313 

314 

315 

316 

■s 

t*j 

31  7 

11 

318 

j 

319 

320 

321 

322 

32  3 

) J 

324 

; :| 

325 

1 1 

326 

1 

327 

■ 1 

328 

329 

:1 

330 

i . 

331 

' i 

332 

il 

333 

ill 

334 

m 

335 

!i| 

336 

. 1 1 

337 

1 1 

338 

:jl 

339 

' I 

340 

1 

341 

342 

^ | ■ 

343 

344 

'll 

345 

! 1'S 

346 

If 

347 

[1 

348 

1 i 

349 

i ■? 

350 

' >; 

351 

352 

■; 

353 

354 

il 

355 

356 

| 

357 

1 l 

350 

■ 

w 

359 

C 


r>  n <j>  r>  n ,1  n o .1  n n 
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C THF  CCMPLFTION  TIME  FOR  ALL  ACTIVITIES  IS  SET  TO  THEIR  LOWER 

C PERCENTILE.  THF  RESULTING  CRITICAL  PATH  TIMF  IS  A LOWER 

C BOUND  ON  THF  FX=FCTFD  CRITICAL  PATH  T I ME  • 

C 

DO  6402  1 = 1, M 
640?  CTIMFI  I ) = FLO<  I ) 

CALL  B I NV  A ( &2800 ) 

6400  CPl  B=  B 1 I NV( 1 .NMMP1) 

WP1TFI6.64C5)  CPLB 

6405  FORMAT ( 1H0 ,5X  ,» A LOWER  BOUND  ON  THE  EXPECTED  CRITICAL  PATH  TIME  IS 

* = '.FI5.5I 

WRITEI6 ,7606)  KKK 
WPITFI6.7707)  ( KB ( I ) » I = 1 , KKK ) 

WRITFI 6,7710)  KKB 
WP!TF(6,7707)  ( IBB{  I ),!  = ), KKB) 

C 

C.  THE  CCMPLFTION  TImF  FOR  ALL  ACTIVITIES  IS  SET  TO  THEIR  UPPFR 

c PFPCFNTTLF.  THF  RFSULTING  CRITICAL  PATH  TIMF  IS  A UPPER 

C BOUND  ON  THF  FXPECTI  D CRITICAL  PATH  TIME. 

C 

DO  6 4C 6 I =1  . M 
64C6  CTIMFI II  = FHI ( I ) 

CALL  flINVA(&2BCO) 

6401  CPUB=  B1INVI1, NMMP1 ) 

WPl TF ( 6 , 64  09 ) CP UR 

6409  FORMAT ( 1HC ,5x t 'A  UPPER  BOUND  ON  THE  EXPECTED  CRITICAL  P*TH  TIMf  IS 

* = • , E 1 5 , 5 ) 

WRITFI 6,7606)  KKK 

WRITEI6.7707)  ( K R ( I ) • I - 1 , KKK ) 

WP I TP( 6 , 77 10)  KKR 

WRITFI  6.7707)  I IRBI  I ) , ’ = 1 , KKB ) 

C 

C FDI I ) = THE  LOWFP  ROUND  ON  THE  EXPECTED  CRITICAL  PATH  TIME 

C PLUS  1/IFDF  OF  THF  DISTANCE  TO  THE  UPPFR  BOUND 

C NLFFDIIR.I)  = THE  SUM  OF  ITHE  CRITICAL  PATH  TIME  FOR  A 

C CONE  1 1.  UR  ATI  ON  * THE  PROBABILITY  OF  THF 

CONFIGURATION  WHFN  THE  CRITICAL  PATH  TIME  IS 

<=  FDI I)  ) * (THE  NUMBER  OF  POSSI 9LF 
CONFIGURATIONS)  / ITHE  SAMPLE  SIZE) 

FOP  THE  IR-TH  CLUSTER 

FD  AND  NLFFO  ARE  USED  TO  OuILD  AN  •EMPIRICAL*  DISTRIBUTION  OF 
THE  CRITICAL  PATH  TIMES 

C=( CPUB-CPL3)/IEDF 
DO  6412  K=) , KCPR 
OB  64 1 2 1=1 , IEOF 

FDI  I ) = CPL  R ♦ I *C 
412  NLEFDI  K, I 1 = 0 .00 

THE  ASSTCIATF  GROUPS  ARF  NOW  FORMED 
C 

WPITEI 6,3165) 

3 1 f 5 FOPMAT ( 1H1 »5X,  • T HE  ASSOCIATES  ARE  NOW  IDENTIFIED:*) 

I I T I ! = l 

DO  2825  I =1  ,M 
2825  CT IMF!  I )=COT{ I ) 

!WWWO= ICRITPI 1 ) 

CHANG=  L AMRDA*  S I GMA ( I WWWO ) 

TFX=COT( IWWtiQ )-CHANG 
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, IF(TFX  .LT.r.C  ) CHANG=COT  < IWNWO) 

CT  IMF(  INWWQ1=C0T!  I W WWQ ) — CHANG 
TALL  8INVA1&23C0) 

23C1  CONTINUE 

IF ( 1ST AH  I toWWQ ) .ED.  I ) CALL  S INV 1 ( G22 825 , COT!  t WWWQ) »CTl ME( INN WO  I . 
*IWWWQ) 

PFDCOS! I NNNO)  = REDCOS ( I NWNQ > +COT ( t N NNO ) -CT IM E C I WWWQ ) 

22825  CTIMF(  IWWWQ)=COT(  IWNNO) 

I NNNQ= I CPITPI I 1 1 I I I 

CHANG=  L AM30 A* S IGm A { IWNNO) 

tfx=cot< i nwno) -Chang 

IF (TEX .LT . o. 0 ) CHANG=COT ( InnNOI 

C T l ME ( I NWNO )=COT(  fW  WWQ ) -CHANG 

IF  ( ISTATI  tWWWO)  .EQ.l  ) CALL  8 I N V 1 < & 2 2 BO 0 , C T I ME<  I NNN Q > .COT ( I N NNO » . 
At  NWNQ) 


C 

C. 

r 

291  C 


2913 

i 2912 
2oll 

t 

29]  5 

i 

i 

l 

i 

» 

i 

2p  16 
| ? 8 1 0 


t C 


RFDCOS ( IWWWQJ^PEDCOS!  l N N NO ) -COT ( I * NNO ) ACT  I ME ( IWWWQ) 

C.n  TC  ->2800 

DETERMINE  ASSOCIATE  GPOUP 

NINAGI  I I I I I 1=0 
DO  291 1 K= 1 , KKB 
KKt  1 

IF( IBDIK) .EO. ICRITP(KK)  ) GO  TO  2911 
IF! KK.CE.KCP3)  GO  Tn  2912 
KK=KK ♦ 1 
GO  TO  2913 

NINAGI  I III  I ) = N I NAG(  III1I1+1 
ASSCRPI  III  TI  , N I N AG!  II  II  I ) )= I BR ( K ) 

CONTINUE 

WRITE! S .2915)  I I II  I , ICR  1 TP< I I I I I ) ,NINAG(  I » I I I ) 

FORMAT! lHf .l^X, 'THE  NUMHFR  OF  ASSoCIATFS  ASSOClATFD  NITH  ThE  '.J3, 
* • — TH  CRITICAL  PATH  ACTIVITY,  I.E.  ACTIVITY  '.13,',  IS  = *.13) 

I DUCK=  N l NAG!  IIIII) 

IF!  I0UCK.EQ.0  ) GO  To  28  I C 

NPJ  TE  ! 8 ,2916  ) I ASS  GPP  ( I 1 I I I , I ) , 1 = 1 , I DUCK  ) 

FORMAT ( 1H0.15X, 'THE  ACTtVJTJFS  IN  THE  aSSOcIATf  GROUP  ARE  AS  FOLLO 
*WS'  ,/, 1 SX ,50!  I 3, • , • ) ) 

III  I 1=1  I I T M-l 

IF!  I I I I I .LF.KCP3)  GO  TO  280  1 
DETERMINE  THE  CLUSTERS 

THE  CLUSTFRS  ARE  POOLED  TOWARD  JHE  TERMINAL  NODE 
NCL'JS  = THE  NUMDEp  OF  NON-EMPTY  CLUSTERS 

NINCL(I)  = THF  NUMRFR  OF  ACTIVITIES  IN  THE  I-TH  CLUSTER 
I NCLu  S ! I , J ) = THE  J-TH  ACTIVITY  IN  THE  I-Tm  CLUSTER 
NCL  t NC ! I ) = THE  NUMBER  OF  CLUSTERS  COMPRISING  THE  I-TH 
CLU8TFP  AFTF9  POOLING 

CLINCLII.J)  = THE  J-TH  CLUSTER  WHICH  HAS  BEEN  POOLED  INTO 
THF  I-TH  CLUSTER 

NCL I NC  AND  CLINCL  HELP  KFFP  TRACK  OF  WHICH  CLUSTER  THE 
CRITICAL  PATH  ACTIVITIES  APF  IN 


8FLOW  FORMS  CLUSTEPS  BY  PUTTING  EACH  CRITICAL  PATH  ACTIVITY  IN 

separate  cluster  and  thfn  aoding  each  critical  path  activity's 

ASSOCIATES  to  ITS  CLUSTER 


4 20 

421 

422 

423 
4 24 

425 

426 

427 

428 

429 
4 30 

431 

432 

433 

434 
4 35 

436 

437 

438 

439 

440 

441 

442 

443 

444 

445 

446 

447 

448 

449 

450 

451 

452 

453 

454 

455 

456 

457 

458 

459 

460 

461 

462 

463 

464 

465 

466 

467 

468 

469 
4 70 

471 

472 

473 

474 

475 

476 

477 

478 

479 


i 


30  25 
3027 


3028 

I ! 


303  0 
r 
c 
c 


3«35 


NClUS=KCpB 
no  3020  T=1,KCPF) 

NCL I NC ( I ) = 1 

CLINCI  ( 1,1  > = I 

N I NCL ( I ) =NINAG<  I ) ♦ 1 

I NCL  US ( 1,1  ) = t CP  I TO  ( I ) 

tr (NJNAG< I ) .PO.O ) GO  TO  3020 

TOuCK=NtNCL<  I > 

on  30  2 1 J=2 , I DUCK 

JJ=J-1 

I NCL UP ( I.J)=ASSGRP( I.JJl 
CONT INUE 

RFl  Ow  POOL*;  CLUSTERS  FORMED  FROM  ASSOCIATES 

T A=0 
I A;  IAA 1 

IF(  tA.GF.KrPO)  GO  TO  3C  3C 
if(  nci  us.ro.  1 ) GO  TO  3C  30 
Inf  A=N  I NCL  ( IA) 

IF( ICTA.FO.nl  GO  TO  3031 
I AA=I A* 1 

00  3023  I I = I A A « K CPU 
IDT I=NI NCL (II) 

IF(  101  I .FQ.O)  GO  Tn  3C23 
00  3025  1=1  • I 0 I A 
00  3i-25  J = 1 , IDII 

IF  ( INCLUS(  I I . J) .FQ. INCLUSt  IA,  I ) ) GO  TO  30  27 

CONTINUE 

GO  TO  3323 

NCL  U 5=  NCL  US— 1 

OO  30?  A J=  1 , 1 13  I I 

DO  3020  1=1 .IOTA 

IF(  INCLUS(  IT ,J).FO. I NCLUS!  !A» I ) ) GO  TO  3028 
CONTINUE 

N I NCI  ( I A ) = NT  NCL ( I A)  + 1 

tNCLUS ( I A, N I NCL < TA>  ) = I NCLUS ( I I . J ) 

CONT INUE 
NT  NCL ( I I ) = 0 

NCL  I NC.  ( I A ) = NCL  t NC  ( I A ) ♦ 1 
CLINCI  ( T A,  NCL  INCH  A)  ) = II 
NCL I NC ( I I ) =0 

continue 

GO  TO  3031 
CONTI NUE 
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BFl  0*  DESCRIBES  CLUSTERS  AFTFR  POOLING  BASED  ON  THF  ASSOCIATES 
WRITE! f, 3033 > NCLUS 

FORMAT < 1 HI  ,10X, ‘THERE  A 'F  *,I3,*  NONEMPTY  CLUSTERS  AFTER  POOLING  O 
AN  THE  F'ASIS  OF  ASSOCIATES  ONLY.*) 

I 1=0 

no  3<-3A  I = 1,KCPB 
I F ( N j NCL (I  ) • FO • C ) GO  TO  3034 
T !=  1 1 ♦ 1 

inuCJ=NlNCL(  I ) 

WP1TE!6.3C35>  I , ( I NCLUS 1 1 . J » * J=1 . IOUCJ  ) 

FORMAT( 1H0.ICX. ‘THE  ACTIVITIES  IN  THE  CLUSTER  ARE  AS  FOL 

ALOWS :»,/,15X,50!l3»*.,>) 

CONTINUF 


526 

527  , 

528 

529 

530 

531 
53? 

533 

534 

535 

536 

537 

538 
839  j, 


c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 


p 3110 
3 If.  2 


.110  7 
3101 

3103 


C 

c 

r 

c 

c 

c 

r 

r 

r 


3122 

3123 


3323 


OESCRIflFS  WHERF  each  activity  is  before  eliminants  are 

CONSIDERED 


examine  EACH  ACTIVITY  AND  DETERMINE  WHICH  CLUSTER. IF  ANY.  IT  IS 
IN. 

LEPT(I)  = 0 IMPLIES  THAT  THE  I-TH  ACTIVITY  IS  NOT  !N  ANY 

CLUSTER 

LEET(I)  = j IMPLIES  THE  I-TH  ACTIVITY  IS  IN  THE  J-TH 
CL  USTER 


WRItE<6,3104 1 

3104  EORMaT(  1HC . I0X , • ThE  CLUSTER  TO  WHICH  EACH  ACTIVITY  BELONGS!  15 X 

* . * I ZERO  IMPLIES  that  THE  ACTIVITY  IS  NOT  In  ANY  CLUSTER)*  I 
DO  3101  I=1,M 

LE^TI I )=f 
DO  3102  J=1,KCPB 
IE ( NINCLI J ) .EO.O ) GO  TO  3102 

i puck  - n i ncl ( j ) 

DO  3110  Ksl.IDUCK 

IF(  I .FO.INCLUSI J.K) ) GO  TO  3107 

continue 
continue 

GO  TO  3101 
LEFT! 1 )=J 

WPtTE(6,3133)  I .LEFT  ( I ) 

FORMAT ( 1 H » 1 SX  * * THE  • . I 3 . • -TH  ACTIVITY  IS  IN  THF  *,T3.*-TH  CLUSTER 
**  ) 

iNDEXL-t 


tc 

Mc 

c 

c 


LEFTOVERS  ARE  ACTIVITIES  NOT  IN  CLUSTERS  AfTEP  ASSOCIATES  HAVE 
HEF  N CONSIDERED  BUT  lEFOWE  ELIMINANTS  HAVE  BFEN  CONSIDERED 


DETERMINE  THE  NUMBER  qF  LEFTOVERS.  NLEFT 

LEFTO(L)  = J IMPLIES  THAT  THE  L-TH  LEFTOVER  IS  THF  J-TH 
At  T I V I T Y 


Nl  EFTO 
DO  3122  J=1,M 

I F( LEFT ( J ) ,NF .0 1 GO  TO  3122 
nleft=  nleftm 

LEFTCt  NLFFTl-J 

continue 

WPI TE( fi ,31 23 » NLFFT 

format  i iho. iox. *there  are  *,13.*  activities  not  in  any  cluster  ye 

*T  . • ) 

WRI TE( e , 3323  > 

FORMAT!  1HI ,SX, * THF  ELIMINANTS  OF  EACH  NON-CRI  T I C AL- PATH  ACTIVITY  A 
*RE  NOW  determined:*) 


ELIMINANTS  FDR  EACH  MON-C R I T I C AL- PATH  ACTIVITY  ARE  NOW 
DETERMINED 

NNNCP  = THE  NUMF'FP  OF  ACTlVlTlFS  NOT  ON  THE  CRITICAL  PATH 
NONCPILE)  = THE  LE-TH  ACTIVITY  NOT  ON  THE  CRITICAL  PATH 


KNNCF=M-KCP3 

LF=r 

DO  5000  1=1, M 


540 

541 
54? 

543 

544 


>45 


546 

547 


546 


549 

550 


551 

552 

553 

554 

555 

556 

557 

558 

559 

560  ' 

561 

562 

563 

564 

565 

566 

567 

568 

569 

570 

571 

572 

573 

574 

575 

576 

577 

578 

579 

580 

581 

582 

583 

584 

585 

586 

587 

588 

589 

590 

591 

592 

593 

594 

595 

596 

597 

598 

599 


it&L:  1:6 J,  - 


1 


J=1 

5001  IF(  I .FQ.ICPITPl J))  GO  TO  5000 
J = J+  1 

IFI J.LE.KCPB)  GO  TD  5301 

5002  LF=LF+1 
NONCPI LE>=! 

5oik  CONTINUE 

WPTTFI6, 50051  NNNCP 

5005  FORMAT  l 1H0  o5x  .*  THERE  APF  * , I 3 , • ACTIVITIES  NOT  ON  THE  CRITICAL  PA 
♦TH.  THFY  *RF  AS  FOLLOWS:*) 

IFINNNCP.FG.O ) GO  TO  3124 
DO  5006  1=1, LF 

5006  WRITE! 6.5C07)  I, NONCPI II 

500  7 FOPMATI1H  , 1 5 X , I 3 » * • • • I 3 ) 

IFINNNCP.FO.O ) GO  TO  3124 
LF=0 

3126  LF-rLFM 

IF  ( ISTATI  IWWWO) .F3.I  ) CALL  B I NV 1 ( £2 3 1? 7, COT I I WWWQ 1 . CT I ME< IWWWO ) • 
♦IWWWO) 

PFOCOS  IIWWWO)  = RF  DC  OS  ( IWWWQ)“CTI  ME  ( IWWWO)  4-COT  ( IWWWO ) 

23127  CONTINUE 

CTTMFI  I W W W 0 ) = COT  I I W WW  Q ) 

CT I me t NCNcPILE ) ) = COT  ( NONCP ( L E I ) * THFT A *S I GMA ( NON CP t LE  > ) 

IF  ( ISTATI NONCP(LEl)  ,EO . 1 ) CALL  B I NV 1 I S7756 . CT T ME ( NONCPI LE ) ) • 

* COT  ( NCNCPf  LF1)  ,N  DNCPILF  ) ) 

RFOC Of ( NCNCPILE ) ) =RFDCOs I NONCPI  LE  ) )-COT  INONCPIlF)  ) 4-CT  I ME  I NONCPI  LE  ) 
♦> 

7756  IWWWO  = NCNCPILE) 

WRITE  I 6, 3 15?)  NONCPI  LF)  , CT I ME ( NONCP ( LE ) ) 


FORMAT! 1HC.///.  5X , • THE  COMPLFTION  T I MF  FOP  THE 
* HAS  HFEN  CHANGED  TO  *,E15.5) 

GO  TO  228CC 
CONTINUE 


•TH  ACTIVITY 


OFTERMINE  THF  FlImInAnTS  OF  THF  LE-TH  ACTIVITY  NOT  ON  THE 
CRITICAL  PATH 

NE  = THE  NUMBER  OF  FLIMINANTS  FOR  THE  LE-TH 

ACTIVITY  NOT  ON  THF  CRITICAL  PATH 
FGRRIJ)  = r HF  J-TH  FLIMlNANT  FOR  THF  LE-TH  ACTIVITY 

NOT  ON  THE  CRITICAL  PATH 


NF=r 

DO  3133  K=1 i K CPO 
DO  3131  1=1, KKH 

IF!  TRBI  I 1 .EG.  ICRITPIX  ) 1 GO  TO  3130 
3131  CONTINUF 


SBBRSS 


NE=NE+ 1 

FGPP(NE)=ICR1TPI<) 

CONTINUE 

WPITFI 6,31 33)  NE , NONCPI  LE ) 

FORMAT  I IHO ,10X , 'THERE  ARE  *.I3,*  FLIMINANTS  CORRESPONDING  TO  ACTI  V 
* I TY  • , I 3) 

IFINE.EO.O)  GO  TO  3171 
DO  3r>5  K=1,NE 

WRITE! 6 ,31 36)  K, NONCP  It  E ) , tGpp ( K ) 

rORMATIlH  , 1 4 X , * THE  *.13.*-TH  FLIMlNANT  CORRESPONDING  TO  ACTIVITY 
* • , I 3 , * IS  ACTIVITY  *,131 

OFTERMINE  WHETHER  NONCP(LE)  IS  AN  ASSOCIATE 
JA  = 1 IF  NONCP ( LF ) IS  AN  ASSOCIATE 


56 


3 182 


I 

I 

I 

I 

I 1103 

f ' 

j 


3 1 8 A 
3172 

3171 


n 


M71 


■ t 

til1*74 


JA  = 2 IP  NONC  P ( L F ) IS  NOT  AN  ASSOCIATE 

K=NCNCP(LE) 

JA=1 

IF(LFFT(K  ) .EO.Ol  JA  = 2 
IPIJA.FO.?)  GO  TO  S010 
I T = LFPT (K  ) 

THF  IT-TH  CLUSTER  IS  FXPANOFO  TO  INCLUDE  EL  I M I NA  NT  S 

GO  TP  S011 
CONTINUE 

ITTTT  IS  THF  ACTIVITY  NUMB  F P OF  THE  FIRST  FLTMINANT 

IT  IS  THF  CLUSTER  Tn  WHICH  The  FIRST  ELlMtNANT  CURRENTLY  BELONG 

ITTT=rGRP< 1 1 
lT=l.EF  T( ITTT) 

LFFT (NPNCP(LE) ) = I T 

THE  IT-TH  CLUSTER  IS  FXpANOEO  TO  INCLUDE  FL I M I NANTS 

N I NCI.  ( I T ) = N I NCL  < I T ) ♦ 1 
INCLUSI  IT.MNCK  IT)  )=NCNCP(LF) 

IFINE.ED.I ) GO  TO  3171 
DO  3172  J=JA,NE 

IU  IS  THE  ACTIVITY  NUMBFR  OF  THE  NEXT  ELIMINANT 

if  iu  is  in  cluster  <,  thfn  cluster  k is  pooled  into  cluster  it 

IU=FGRP( J) 

K=LFFT{ IU) 

IF( IT.FO.K)  GO  TO  3172 

NCLUS=NCLUS-1 

I W=NCL I NC(<) 

on  3103  I A- 1 « I w 

LFFT( ICR  I TP (CL  I NCL  I K,  I A 1 ) l = IT 

NCI  I Nc  I IT) =NCL I NC ( I T 1 ♦ 1 

CLINCH  IT.NCLINC(Ir))=CLlNCL(K,  I A > 

NCL INC ( K ) = 0 
IW  = MNCL(K) 

Nt  NCL ( K Is" 

DO  3104  1A=1,TW 

LF  F T ( I NCLUSIK , I A) * = I T 

N I NCL ( IT  ) = N I NCL  C I T I ♦ 1 

I NCL  US  ( IT,  NINCLI  IT  1 ) = INC.LUS(K  , I A 1 

CONTINUE 

CONTINUE 

IFILF.LT.NNNCR)  GO  TO  3126 

FND  OF  POOLING  BASED  ON  ELIMINANTS  FXCEPT  FOR  THE  FOLLOWING 
DESCRIPTION 

WR I TF( 6 ,31 73)  NCLUS 

FOPMATI1H1 ,05X,*THERF  ARF  *,I3,'  CLUSTERS.*! 

DO  3176  I*1,KCPB 

IFININCH  I l.EO.Ol  GO  TO  3176 

1 DO=NI NCL ( I 1 

WRI  TF  < 6, 31  74)  NlNCl.  ( I ) • I • ( I NCLUS  < I . J)  • J*I  , |OD  ) 

FORMAT!  |HO,|OX,*THfcHF  API:  *,13,*  ACTIVITIES  IN  THF  *#|1»*-TM  CLU5T 


660 

661 

662 

663 

664 

665 

666 
66  7 
660 

669 

670 

671 

672 

673 

674 

675 

676 

677 
670 
679 
6B0 
681 
682 

683 

684 

685 

686 

687 

688 

689 

690 

691 

692 

693 

694 

695 

696 

697 

698 

699 

700 

701 

702 

703 

704 

705 

706 


707 

700 

709 

710 

711 

712 

713 

714 

715 

716 

717 

718 

719 


*EP 


57 

. THEY  ARE  AS  FOLLOWS:* ,/,2CX.50< 13. *.*> ) 

I OUC K = NCI  I NC  ( I ) 

WRITF(6 ,31 76)  NCL I NC ( I ) , (CLINCl  ( I , J>  , J=1 , I DUCK  > 

317b  FORMAT ( 1H0 , 1 5X ,! 3,  • CLUSTERS  HAVE  RFEN  POOLFD  TO  MAKE  THIS  CLUSTER 
*.  THcr  *ERE  AS  FOL-OWS ! • ./ .20X , 50< 13. * . * > > 

3176  CONTINUE 

IF  (SAMSIZ.LE.O)  WRITE! 6.6045  ) 

IF  (SAMSIZ.GT.C)  WRI  TF(  6.6056)  SA MSI Z , SAMS!  Z 
6056  FOTMATI  5x  . * THE  FOLLOWING  TABLES  WERE  DETERMINED  CONSIDERING  AT  M 
*OST  * , I 8 » • ACTIVITY  CONFIGURATIONS  PER  CLUSTER. *,/.  1 1 X .* tF  THERE  A 
*RF  NO  MOPE  THAN  ',13,'  ACTIVITY  CONFIGURATIONS  IN  A CLUSTFR. * . / . I I 
*X,<  THEN  all  activity  CONFIGURATIONS  arf  explicitly  considered,  and 

* NO  sampling  is  done.** 

6046  FORMAT  ( 5 X « * THE  FOLLOWING  TABLES  WERE  DETERMINED  CONSIDERING  ALL 

* ACTIVITY  CONFIGURATIONS.*) 

C 

C EARNING:  THE  SYSTEM  SUBROUTINE.  CLOCK  MAY  NOT  BF  A PART  OF 

C -.LI  SYSTEMS.  XR AN  NEEDS  TO  BF  A RANDOM  SEED. 

C 

CALI  CLOCK(XRAN) 

IYUTS  = XR  AN 

WP I TF  (6.3238)  IYUTS 

3?38  format  < i ht  .sx . * the  initialization  parameter  for  any  sampling  is  i 

*Y  = • , I 1C  ) 

r 

C STATFMFNT  number  3124  MARKS  THE  END  OF  POOLING  CLUSTERS  BASED 

C ON  LEFTOVERS  and  eliminants 

r. 

3124  CONTINUE 

I F ( CPTCNl  .EQ.I1  GO  TO  RR9 

c 

c 

C the  final  CLUSTFRS  HAVE  NOW  BEEN  DETERMINED 

C ThF  ?*4NINCL<!)  RUNS  APE  NOW  MADF  FOR  ALL  I WITH  NINCL(I)>C. 

C 

C 

6008  IF < ICBCP.EO. I ) !C8CP=? 

I NDFXL=2 

I F ( CPTON2  .EO.l)  GO  TO  877C 1 
IF( ICBC°.E0.?)  GO  TO  87701 
DO  17701  1=1. M 

) 17701  COT(I)=FLO(I) 

87701  CONTINUF 
IR  = 0 


IR= IRF l 

IF< IR.GT.KCOR)  GO  TO  32C8 

!F(MnCL(  IP)  .EO.O)  GO  TO  3230 

DO  6031  1=1.10 

MOMFNT ( IP  , I ) = C .03 

DO  R1C1  1 = 1.1  EOF 

NLFFDI  IR,  I ) = «'  .00 

IP  = C 

N!  R ( T R) =C 


NUMRFR  RF  ACTIVITY  CONFIGURATIONS  IN  THE  SAMPLE 
VECTOR  CONTAINING  THE  UPPER  BOUNDS  ON  THE  NETWORK 
COMPLETION  TIME  DISTRIBUTION  TO  BE  averaged  WITH  THE 
LOWER  BOUNDS  ON  THE  NETWOPK  TO  YIELD  THE  AVERAGE  NETWORK 
CCMPIFTION  TIME  DISTRIBUTION. 


58 


c 


r 

c 

c 

c 

r 


3?2? 


C 

C 

r 

r 

3?0* 


603? 

9900 


6 4?  C 


1000  1 
20000 


r 

c 

c 

r 


20500 


621  0 
621  1 


20501 

r 

c 

c 

r 

i 

5 

i 

) 


c 

c 

*r 


IC=NINCL ( !R  > 

1B=2**NINCL < IR ) 

I DOALL  = 1 MEANS  ALL  ArTIVITV  C C NF I GURA T I ONS  ARF  EXPLICITLY 
CONS IDCREO. 

f DC ALL  = 0 MEANS  TO  SAMPLE. 


I DOALL  = 0 

IF  (SAMSI7.LF.C.OP.SAMS  IZ.GE. tfl)  IOOALL=l 
DO  3222  1=1, M 

CT IMF(  I > = COT ( I 1 

C.O  TO  ?ocoo 

STATEMENT  32C4  IS  I HE  RE-ENTRY  POINT  FROM  THE  SIMPLEX 
AIG3RITHM  WHEN  CLUSTER  BASED  BOUNDS  ARE  BEING  COMPuTFD 

continue 

DO  6032  I =1  . 1 C 

MOMENT ( I R , I 1 = MOMtNT(IP.I)  ♦! B 1 I NV<  1 , NMMP 1 ) ** I 1 *SPROB 
X=  B 1 I NV ( 1 » NMMP 1 1 
X=  X-l . 00-10 
1=0 
1 = 1+1 

IF( x .GT  .FO(  I > 1 GO  TO  (,42c 
NLEFOI  IR.  1 )=NLEEr><  IR,  I 1 +SPPUB 
CONTINUE 
IP=IP+ 1 

N I B ( IP)=NIB< IR )+| 

GENERATE  NFXT  ACTIVITY  CoNF I GUHAT I CN  TO  BF  EXPLICITLY 
CONSlOFPEO. 

IF  ( ID3ALL.E0.01  GO  TO  205C0 
IF  (IP.GT.I13)  GO  TO  320  7 
PANSAM  = IP 
GO  TO  20501 

IF ( N In ( IP  1 .GT. SAMSIZ ) GO  TO  32C 7 

IYUTS  = I YUTS*655 39 

IF  (IYUTS)  6210,6211.6211 

IYUTS  = I YUTS  + 21 A74R364 7 + 1 

XRAN  = IYUTS 

XR  AN  = XRAN*. 466661 3F  — 9 

PANSAM  = XRAM*DFLOAT(  IU-1)  ♦ 1 

CONTINUE 

CONVERT  THE  PANDOM  NUMBER,  RANSaM,  to  a binary  number  to  define  an 
ACTIVITY  CONFIGUPAT ION, 

KRAN  = RANSAM 
SPPOB= 1 .DO 
OO  B5C  5 1 = 1,  IC 

l HALF  = K R AN/ 2 
17  = KRAN  - IHALF*2 

L = INCLUS(IR.I) 

CTIMEILl  = tZ*FHt(L)-IZ*FLO(L>  ♦ FLO(L» 

SPROR=  SPROB* ( 1 7*PQ(L  > ♦(  1 -IZ ) *PP(L  1 1*2,00 

SPROH=  '5*  *N  I NCL  ( I R I * THE  PROBABILITY  OF  THIS  CONFIGURATION 


797 

788 

7H-» 

790 

791 

792 

79  3 

794 

795 

796 

797 

798 

799 

800 

80  1 
80? 

803 

804 

805 

806 

807 

808 

809 

810 
81  1 
812 

813 

814 

81  5 
816 

817 

818 

819 

820 

82  1 
822 

823 

824 

825 

826 

827 

828 

829 

830 

831 
8 32 

833 

834 
935 
8 36 

837 

838 

839 


I 


1 


l OOP'1 

r. 

c 

c 


KRAN  = IHALF 
CALL  BINVA(f,2«00> 

NIP( IP)=NIO( IR)-1 

on  fiD’i  i = t , 1 o 

MOMENT ( ! P « I ) = MOMENT ( I R. I >/NI B{ IR ) 

GO  Tr  32CC 
WRITE  (6.6362) 

FORMAT  ( 1 M 1 ) 

IT  (TCOCP.FQ.2)  GO  TO  6011 

on  irr  oo  j = i , l o 

I T M A X = 0 

TM AX  = 0 . 

00  602 1 1=1 ,KCPB 

IF (NINCL( I ) .EQ.O)  GO  TO  5021 

IF  ( MOMENT < I . J > .LF.TMAX  ) GO  TO  EC2I 

tmax  = mcmf:nt(  i , j> 

ITMAX= I 
CONT  TNUE 

WRITE  (6,6023)  J . J , MO MF NT ( I T M AX , J ) 

FORMAT  (1H0,6X,*A  LOWFR  BOUND.  T- ( • , I ? , • 5 THFT A , L AMBDA ) . ON  THE  • 
*.  I?  » * -T  H MOMENT  OF  THF  NFTWORK  COMPLETION  TIME  = • .E15.5) 

conttnuf 


BFGIN  THF  PROCEDURE  FOR  DETERMINING  UPPER  ROUNDS  ON  THE 
NETWORK  COMPLETION  TIME  DISTRIBUTION 


0011  CONTINUE 

00  <3007  I R — 1 , K CPR 
!F(N INCL  ( IP  ) .FO.O ) GO  TO  9CC7 
00  oqon  1=2. I EOF 
11=1-1 

9990  NLFFD(  IR,  I ) = NLEFO(  Ir.  I ) tNLEFDC  IP.  II  ) 

WNNJ  R = OFLF AT (NIP ( IR ) ) 

DO  ToROO  1 = 1 , IFDF 

39990  NL  FFO ( IR.  I )=NLEFD(  IP,  I )/WNNIR 
9007  CONTINUE 

DO  101 1 1 IW=1 , KCPB 
IF ( NINCL ( IR  ) .GT ,P ) GO  TO  IP112 
10111  CONTINUE 
1011?  IPR=  I P 
C 

r i rp  = non-fmpty  cluster  with  the  smallest  inoex 


DO  101  19  19  = 1 « K CPR 

IF(NINCL( IR1.EO.O)  GO  To  10119 

DO  1 01  IP  1 = 1 . IFDF 

I F ( NLF F 0 ( IR,  I ) ,LT.NLEFD(  IRR, I ) ) NLEFD(  IRR.I ) = NLFFD(  IR.I  I 
10110  CONTINUE 
10119  CONTINUE 

WRITF  (6.626A) 

WPITF( 6.9A23) 

9423  FORMAT ( 1H0 .5X . 'AN  UPPER  BOUND  ON  THE  NETWORK  COMPLETION  T I MF  DISTP 

♦ irution:  f*-(  . ;thfta  jlamhoa ) • ) 

DO  942  1 1=1  . I EOF 

NSAVEI  I )=NLEFD( IRR,  I ) 

X = N1  EFO(  IRR.  I ) 

9421  WR1  TF(  6.94??)  rr)(  I ) . THr  T A • L AMBDA  . X 

942?  FORMAT ( 1 7X«*F4(  • ,F1 h. 6, .E 15.5. .FI5.S. • I ■ SElS.SI 


60 


c 

r 

c 

c 

r 

r 

c 

c. 

C 

C 

C 

c 

c 

c 

c 

r. 


600  i 

6011 

6264 


7 101 
1 000  9 
601? 


6900 


49990 

66900 

F 

! 

j 

?0  11  0 
?0l  15 


6 42.1 

$ 

I 


6421 

6422 

. r 

i 


BEGIN  THE  PROCEDURE  FOR  DETERMINING  UPPER  BOUNDS  ON  THE  MOMFNTS 
CF  THF  NETWORK  COMPLETION  TIME  AND  LOWER  BOUNDS  ON  THE 
DISTRIBUTION  OF  THE  completion  TIMES. 


THE  UPPER  BOUNDS,  T* ( R . THFT A .LAMBDA > . ARE  NOW  DETERMINED. 


FOP  THF  SAKE  OF  NUMERICAL  ACCURACY  THE  P£RT  PROBLEM 

WITH  new  activity  times  is  initially  solved  from  scratch 

INSTEAD  OF  UPDATING  AN  OLD  SOLUTION.  AFTER  THIS 
PFINITI ALI NATION,  THE  RFMAINING  CPITICAL  PATH  TIMES  ARF 
DETERMINED  flY  UPDATING  THIS  SOLUTION. 


ICBCP=  1 
DO  6001  1 = 1. M 

CT IMF!  I > = FHI ( I ) 

COT ( I ) =FHI ( I ) 

GO  Tn  69 1C 
WRITF  (6.6?64> 

FORMAT ( /// ) 

DO  1 0 c C 9 J=l,  10 
ITMIN= I RP 

T M I N=MC  ME  NT(IPP.J) 

DO  71  Cl  1 = 1 .KC.PB 
IFIMNCLII  l.EO.O)  GD  TO  7101 
IF(MOMENT(  I , J)  .GT.TMIN1  C.D  TO  7101 
I T M I N=  I 

T Ml N=MC  MENT I I , J) 

CONTINUE 

WPITF  (6,6012)  J , J , MO  ME  NT(  ITMlN.J) 

FORMaT  (1HO  .5X  ,' AN  UPPER  BOUND,  T •*•(•.  12  ,»  {THETA, LAMBDA!  , ON  THE  • 
*»l2,*-TH  MOMFNT  of  the  NETWORK  COMPLETION  TIME  = *.E15.S) 

00  66900  IR=l,KCPO 
IF <NINCL< IR) .EO.o)  GO  TO  669CP 
DO  6RC0  1 = 2,1 FDF 
11=1-1 

Nl  FFD(  I R,  I )=NLFFD{ Ip, I ) ♦NLEKDI  1R. I I ) 

W NN I B=  DFLOAT(NIBIIW)  ) 

DO  499  »0  I =1  , I EDF 

Nl FFD ( f R, I ) = NLEFD(  IP, I >/WNNT  B 

CONTINUE 

DO  201  I 5 I R= 1 ,KCPB 

IF(NINCL(  IR)  .EO.O)  GO  TO  20115. 

DO  20110  1=1 . I EDP 

TFINLFFCI  IR,  I ) .GT.NLFFDl  I PR . I ) ) NLEFD( IRR, ! I=NLFFD<  IR,  I I 

CONTINUE 

CONT INUE 

WRITr  (6,636?) 

WR1TFI  6.6  423  ) 

FORMAT  ( 1 HC  ,5X  , • A LOWER  BOUND  ON  THE  NETWORK  COMPLETION  TIME  DISTR 

• ibuticn:  f-(  . ;theta;lam  ida)»  ) 

nn  6421  1= l , I FDF 

X = MFFC  I IRR,  I 1 

WO  I TE  I 6,6.4  2?)  FD(  I ) .THFTA.LAMBDA.X 

format ( I 7X , «F- ( • ,F!5.8* • ; • ,E 15. 5.*  I • .E15.S, • 1 * »,FI5.5) 


n 


•i 


■ -ij.  -V-vAV  u-  iRfrwyl 


r»  r»  n n r>  o n r>  .on 


C 

c 


THr  APPROXIMATE  NETWORK  COMPLETION  TIME  DISTRIBUTION 


61 


WRITE  (6.6 16?  ) 

WR  I TE(  S , S4  72  ) 

0472  FORMAT ( 1 HO ,SX .• AN  APPROXIMATE  NETWORK  COMPLETION  TIME  DISTRIBUTION 
*:•.//.  1SX.  *F( .ITHFtA.LAMBOA)  = .5  * < F ♦<.; THET A .LAMBDA ) ♦ F-C.JTH 

♦ FTA.LAMODa  > )*,//) 

00  9471  1 = 1 , I F OF 

X-  .SDC  ANSAvP  ( I ) ♦ . 500  *NLEFD  < I RR  • I ) 

9 471  wr  1TP(6.R473)  FO ( I 1 , THE T A . L AMBD A . X 

7473  FORMAT ( 1 7x  . • F ( 4 ,E  1 5 . 5,  • J 4 , t 1 5 . 5 . 4 5 4 , E 1 5.  5 . 4 I = 4.F15.5> 

999  Wf?I  TF(  6 ,8S~  ) 

STOP 
F NO 

SUBROUTINE  4tNVA<*> 

IMPLICIT  RF4L*8  (A-H.O-Z) 

COMMON  41  I NV . 9F DCOS , C T I ME , X B 1 . I NB AS F . IHFA D , I T A I L . NM MP 1 .NMM.N. I STAT 
COMMON  M.MPt 

DIMENSION  1ST  AT < IOC ) , IHFADI  60  > , ITAI  L (60 ) .XB1I 41) 

01  MENS  I ON  41  I NV ( 4 1 ,41  ) , I NJASEI 4 C ) .C  T I MF ( 1 00 > ,REDCOS(  100  I 
C 

C UpO^TE  THE  FIRST  ROW  OF  R1TNV  AFTER  CHANGING  CTIME 

r 

DO  1 t=2,NMMPl 
BI INV( 1 . I > -0.0 

DO  1 J=2,NMMP1 

1 BllNV(l.I)  = BllNV(l.I)  ♦ Ftl  INVI  J.  I)  *CT  ImEC  InBASEI  J-l  ) > 

U°D  A T T VAl  UE  OF  THE  OBJECTIVE  FUNCTION 

XB1 ( 1)  = 01 INVI 1 ,NMMP! ) 

RFTURNl 
FNP 

SUBPUUTINF  4INV1  I 4 , T MNTW , T MOL D , l D > 

IMPLICIT  RPAl*B  ( \-H, 0- 7 ) 
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COMPUTE  THF  RFDCOS  CORRESPONDING  TO  ONE  CHANGE  IN  CTIME 

DO  2 1=1 ,NMM 

? IF ( IN44SFI I ) • FQ . ID)  11*141 
DIFF  = TMNFW-TMOLD 

TMNEW  IS  THE  NFw  TIME  AND  T MOL  0 IS  THE  OLD  TIME  CORRESPONDING 
td  thf  single  change  in  ctime 

DO  1 K = I « M 

IF( 1ST AT(K) .EO. I ) GO  TO  l 

RFDCOS ( K ) = RFOCOSI K ) -D I FF*(  01  INV(  l I • IHEADI  K>M>  - 

* 41 INVI I t . I TAIL!  K ) 41  ) I 
1 CONTINUE 

OP  4 K = MP 1 , N 

IF  I ISTAT(K) .FO. 1 ) GD  TO  3 

RFDCOS(K)  = RFDCOS! K ) - 0 IFF4D I I NV < I I • K-M4 I > 

3 CONTINUE 

UPOATF  THE  FIRST  ROW  OF  BIINV  AFTFR  CHANGING  CTIME 
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New  Subnetwork  Analysis  Program:  Flowchart 


Read  subnetwork  description 


Set  each  activity's  completion  time  to  its  mean 


Determine  the  activities  on  this  original  critical  path 


Read  the  algorithm  parameters  6,  X 


Set  each  activity' e completion  time  to  its  lower  point 


Determine  the  activities  on  this  critical  path.  (This  critical  path 
completion  time  is  an  lower  bound  on  the  subnetwork  completion  time.) 


Set  each  activity's  completion  time  to  its  upper  point 


Determine  the  activities  on  this  critical  path.  (This  critical  path 
completion  time  is  a upper  bound  on  the  subnetwork  completion  time.) 


Determine  the  associates  of  each  original  critical 
path  activity  and  the  initial  clusters 


Pool  clusters 


Determine  the  eliminants  of  each  non-critlcal  path 
activity  and  pool  the  new  clusters 


Describe  the  final  clusters 


Determine  T~(«;  8,  X)  and  F (•;  8,  X) 
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APPENDIX  C 

Original  Subnetwork  Analysis  Program 

The  Original  Subnetwork  Analysis  Program  is  an  implementation  and 
extension  of  the  analytical  procedure  described  in  Section  3 of  Technical 
Report  No,  48.  The  basic  required  input  is 

(a)  an  acyclic  network  with  one  source  and  one  sink, 

(b)  two  points  from  each  component  activity's  completion 
time  distribution,  and 

(c)  specified  values  for  the  algorithm  parameters  6 and  X. 

The  output  is  mainly 

(a)  upper  and  lower  bounds  on  the  moments  of  the  network 
completion  time,  T*(0,  X)  and  Tr(8,  X)  r ■ 1,  2,  ...,  10; 

(b)  upper  and  lower  bounds  on  the  distribution  function  of  the 
network  completion  time,  F+(*;  9.  X)  and  F (•;  9,  X);  and 

(c)  an  approximate  network  completion  time  distribution, 

F(.;  0,  X)  - l/2[F+(s  0,  X)  + F~(-;  6,  X)]. 

The  main  extension  of  this  program  is  the  inclusion  of  an  option 

n 

to  consider  only  a random  sample  of  the  2 c activity  time  configurations 

for  a cluster  C instead  of  explicitly  evaluating  the  critical  path  time 
n 

for  all  of  the  2 c activity  time  configurations. 

The  basic  computational  technique  for  determining  critical  path 
times  is  the  Simplex  Algorithm.  This  algorithm  is  applied  to  the  dual 
problem.  The  Simplex  Algorithm  is  used  instead  of  the  standard  network 
analysis  techniques  because  the  Simplex  Algorithm  is  ideally  suited  for 
the  type  of  parametric  programming  required  to  evaluate  several  critical 
path  times  when  only  the  activity  times  vary  from  one  problem  to  the  next. 


A listing  of  the  Original  Subnetwork  Analysis  Program  and  a program 
flowchart  are  given  at  the  end  of  this  appendix. 

Specific  Input  Instructions: 

Card  1.  Col.  1-3  : The  number  of  activities  in  the  network,  Format  (13) 
Col.  4-6  : The  number  of  nodes  in  the  network,  Format  (13). 


For  each  activity  one  card  with: 

Col.  11-15:  The  origin  node  of  the  activity,  Format  (15). 

Col.  21-25:  The  terminal  node  of  the  activity,  Format  (15). 

Col.  31-40:  The  lower  point  on  the  activity's  completion  time 

distribution.  Format  (F10.0) 

Col.  41-50:  The  upper  point  on  the  activity's  completion  time 

distribution,  Format  (F10.0) 

Next  Card.  Col.  1:  0FT0N1.  0PT0N1-1  implies  that  the  program  will 

terminate  after  the  clusters  have  been  formed 
on  the  basis  of  associates  and  eliminants. 

OPTONi^l  implies  that  the  program  will  follow 
the  normal  procedure. 

Next  Card.  Col.  1-3:  IEDF.  The  program  computes  an  absolute  upper 

and  lower  bound  for  the  network  completion  time. 

This  range  is  subdivided  into  IEDF  equal  parts 
and  the  approximate  distribution  function  (F+,  F , t) 
velues  ere  printed  et  eech  of  these  dividing 
points.  IEDF  would  usually  be  between  10  end 
100.  IEDF,  Format  (13). 
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Next  Card.  Col.  1-5  : 8,  Format  (F5.2). 

Col.  6-10:  X,  Format  (F5.2) 

Next  Card.  Col.  1-10:  SAMSIZ.  The  number  of  activity  time 

configurations  to  be  randomly  selected  for 
explicit  consideration  in  each  cluster 

n 

analysis.  If  SAMSIZ  < 0 or  SAMSIZ  > 2 c, 
all  activity  time  configurations  will  be 
explicitly  considered  - no  random  sampling 
will  be  done.  Format  (110). 

The  nodes  should  be  numbered  1,  2,  ...»  n with  the  source  being 
number  1,  the  sink  being  nunfcer  n,  and  the  other  node  numbers 
being  arbitrary.  The  activities  should  be  numbered  1,  2,  ...  in 
any  order  desired. 


Current  Dimension  Restrictions: 


Currently  the  program  is  dimensioned  for  a maximum  of 
60  Activities 
40  Nodes 
25  Clusters 

25  Activities/Cluster  and  IEDF  < 500. 


Example: 

The  program's  input  and  output  are  illustrated  in  terms  of 


the  network  in  Figure  C-l. 
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ORIGINAL  SUBNETWORK  ANALYSIS  PROGRAM 
IMPLICIT  REAL*fl  (A-H.O-Z) 

FOR  THE  SAKE  OF  IDENTIFYING  THF  APPROPRIATE  DIMENSIONS.  LET 
M = THE  NUMBFR  OF  ACTIVITIES  IN  THE  NETWORK 
NMM  = NUMBER  OF  NODES  IN  THE  NETWORK 
NMMPl  = NMM  «■  1 
N = M ♦ NMM 

L = THF  LENGTH  OF  THE  CRITICAL  PATH 
C = THF  MAXIMUM  NUMBER  OF  BRANCHES  IN  A CLUSTER 
IEDF  * THE  NUMBER  OF  DIVISIONS  IN  THF  EMPIRICAL 

distribution  function 


INTFGFR  TAIL!  M),HEAD(  m».ASSGrP(  L.L  > . C L l NCL ( L ,L > , EGRP ( L > 
DIMENSION  NINCLI  O.INCLUSI  L.  C).  NCL  INC!  L) 

DIMENSION  NLEFOI  IEDF) . FDI  IEDF  ) .NSAVC  I I EOF  I 
DIMENSION  iZZ(  C).AVG(  d.thatil) 

OIMFNSICN  I NB  AsE ( NMM ) , X NODE I NMM  ) 

DIMENSION  LFFT(M  ),LEFTO<M  ).NONCP(M) 

DIME NS  1 CN  X B 1 (NMMP1  ) , Y1  ( NMMP1  ) , REDCOS(N) , I STATI N ) 

DI  MFNSI CN  TCRI Tp  (L)  ,Nt  NAG( M) , ICRITNI  L + I » . CT IME( N),COT(N) 
DIMENSION  KBIL.  ),IBB(L),F25(M),F7S(M),sIGMAIM),B1  INV  ( NMMP1  , NMMPl  ) 
REAL  MOMENT (_.  1 0 > 

OF  COURSE  THESE  DIMENSIONS  ARE  MERELY  UPPER  BOUNDS 


COMMON  e 1 I NV.BFOCOS.CT  I ME .XBl  . INBASF .HF Ad. TA  IL. NMMPl  . NMM , N . I ST AT 
COMMON  M.MPI 

INTEGER  TAIL! ft D)  ,HEAD(60  J . A 5SGRP f 25 . 2 5 > .CLINCLI  25. 25) .EGRPI 25) 

integer  s amsi z . ransam 

DI  MF  NS  I CN  Nl  NCL (25) , l NCL  uS  ( 2 5 , 2 5 ) , NCL INC! 25) 

DIMENSION  FOI500). NLFFQ ( 50 C ) ,NS AVEl 500 ) 

DIMENSION  AVG ( ?S)«THAT(  25) 

DIMENSION  I NB ASF ( 40 ) 

DIMENSION  XNOOE(AO) 

DIME NS  ION  XBl  < 4 1 ) ,Y  1 ( 41  ) ,REOCOS<  ICO ) . I STATI 100) 

DI MENS  I CN  51  INV (4  1 ,41 ) .KB<25> . IBB (25  > , FL O ( 60  ) ,EH I ( 60 ) , S I GMA  ( 60) 
DIMENSION  ICRI TP (25) , N I N AG ( 50 ), ICRITN(26)  . COT ( 1 00 ) , CT| ME I 100) 
DIMENSION  left ( 60 ) , LEFTO ( 60 >.NONCP(60> 

RFAL*B  LAMBDA, MOMENT(25, 10 ) 

M = THE  number  of  ACTIVITIES  IN  THE  NETWORK 
NMM  = THE  NUMBER  OF  NODES  IN  THF  PERT  NETWORK 
READ ( 5 » 100)  M.NMM 
FORM  AT (213) 

n=nmm+m 

MPI=M+ I 
NMMR 1=NMM+ 1 


the  ACTIVITIES  ARE  DESCRIBED  IN  TERMS 
11=  THE  TAIL  NODE,  THE  ORIGIN  NODE 
JJ=THE  HEAD  NODE,  THE  TERMINAL  NODE 
FLO  * THF  LOWFR  PERCENT  ILF 
EH I = THF  UPPER  PERCFNTILF 
SIGMA  » FLO  - FHI 


OF  THEIR  NOOES 
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DO  FID  tsl  , M 
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5 ’ '-wr  N 


2 SO  | 


610 


5561  0 

C 

C 

C 

C 

C 

c 

c 

r 

77551 

7756? 


?700 
? 70 1 


2704 
2 7C  2 
C 

c 

c 

c 

c 

-c 

r 

c 

.c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c. 

o 

.c 

c 

c 

f. 


READ  ( 5 » 250  1 1 I I ♦ JJ»FLU(  I I.FHKI  1 
SIGMa(I)  = FH1 < I )-FLO! 1) 

FORMAT ! IOX. I5.5X . 15.  5X,F1C.0,F1C.0> 

COT  ( 1 ) = (FLO(  DFFHIC  I)  )/2. 

CT1MF! I ) = COT ( I ) 

TAIL! I 1 = II 
HEAD ( I )=JJ 

COTIM)  = THE  ORIGINAL  RIGHT  HAND  SIDES.  I .F.  THE  AVERAGE 
OF  FLO  AND  FH I 
DO  666 1 C I =MP 1 . N 
CTTMFI I ) = 0. 

OPTON1  =1  IMPLIES  THAT  THE  PROGRAM  HILL  TERMINATE  AFTER 
THE  CLUSTERS  HAVE  BEEN  FORMED.  NO  BOUNDS  ON 
THE  PROJECT  COMPLETION  T I MF  MOMENTS  OR 
DISTRIBUTION  MILL  BE  DETERMINED. 

OPTONl  NflT = 1 IMPLIES  THAT  THE  NORMAL  PROCEDURE  WILL  BE 
FOLLOWED. 

READ! 5, 776  51  ) OPTONl 
FOPMATI 1011) 

IF ( OPTONl  .EQ.l)  WR ITEI 6.77552) 

FORMAT ( 1H1  . 10X . 'OPT ION1  = 1 AND  THE  PROGRAM  WILL  TERMINATE  AFTER  THE 
* CLUSTERS  HAVE  BEEN  F OR MFD . • , / , 1 1 X , • NO  BOUNDS  ON  THE  PROJECT  COMPL 
*E T ION  TIME  MOMENTS  OR  DISTRIBUTION  WILL  BE  DETERMINED.*) 

RE  AO ( 5 . 100)  IEDF 
WRITE! 6,?7f 0) 

FORMAT! 1HI . 15X, • INI TIAL  INPUT*) 

WRITFI6.27P1  ) 

FORMAT!  1H0.1CX,  ‘ACTIVITY  ORIGIN  TERMINAL  LOWER  PERC.  UPPER  P 
*ERC.  AVERAGE  PERCENTILE  DIFFERENCE* I 
DO  ?7C4  1=1. M 

WP I TF! 6,2702)  I . TAIL!  I ) .HEAD! I ) ,FLO(  I) . FH I I I ) , CT I ME { I) .SIGMA!  I ) 
FORMAT  ! 1H  ,13X,I3,6X.T3,7X,I3.5X.F10.4.3X.F10.4»3X,FI0.6f4X.F|0.4) 

THE  FOLLOWING  INDICATORS  ARF  USED! 


IPARM  = 1 


IPARM  = 2 
IPARM  = 3 
IPARM  > 3 


IMPLIES  THE  CRITICAL  PATH  TIME  WHEN  ALL  ACTIVITY 
COMPLETION  TIMES  ARE  SET  EOUAL  TO  THEIR  AVERAGES  IS 
BFING  DETERMINED 

IMPLIES  THAT  THE  LOWER  BOUND  ON  THE  COMPLETION  TIME 
FOR  THE  SUBNETWORK  IS  BEING  DETERMINED 
IMPLIES  THAT  ThE  UPPFR  BOUND  ON  THE  COMPLETION  TIME 
FOR  THF  SUBNETWORK  IS  BEING  DETERMINED 
WHEN  I NDE  XL=0  IMPLIES  THAT  THE  ASSOCIATES  ARF 
BEING  DETERMINED 


I NDE XL  = D IMPLIES  THAT  INITIAL  CLUSTERS  ARE  STILL  BEING  FORMED 
j ND  F X L = 1 IMPLIES  THAT  THE  L FFTOVFHS . THEIR  ELIMINANTS.  AND 

poolfo  Clusters  are  being  determined 

I NOE XL  = 2 IMPLIES  THAT  THE  2**NINCL(  ) RUNS  FOR  EACH  CLUSTER 
ARF  BEING  MADE  AND  AVER AGFD 


ICBCP  = 0 IMPLIFS  THAT  THE  PROCFDURE  FOR  DETERMINING  UPPER 

BOUNDS  ON  THE  MOMFNTS  OF  THE  NETWORK  COMPLETION  TIME 
AND  LOWFR  BOUNDS  ON  THE  DISTRIBUTION  OF  COMPLETION 
TIMFS  HAS  NOT  BEEN  BEGUN 


60 

61 


62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

1 10 
1 1 1 
112 

113 

114 
116 
116 
117 

lie 

119 


IC.bcp 


ICHCo  = 


IMPLIES  THAT  THE  PROCEDURE  FOR  DETERMINING  UPPER 
ROUNDS  ON  THF  MOMENTS  OF  THE  NETWORK  COMPLETION  TIME 
AND  LOWER  BOUNDS  ON  THE  DISTRIBUTION  OF  COMPLFt ION 
TIMFs  Is  BEING  INITIALIZED 

IMPLIES  THAT  THE  PROCEDURF  FOR  DETERMINING  UPPER 
ROUNDS  ON  THE  MOMENTS  OF  THF  NETWORK  COMPLETION  TIME 

and  lower  bounds  on  thf  distribution  of  completion 

TIMES  IS  EEING  CARRIED  OUT 


= 0 IMPLIES  that  THE  PROCEDURE  for  determining  a upper 
BCUND  ON  thf  NETWORK  COMPLETION  TIME  DISTRIBUTION  HAS  NOT 
BEGUN 

= i implies  that  the  uppf.r  bound  on  the  network 

COMPLETION  TIME  DISTRIBUTION  IS  BEING  DETERMINED 


IPARM= I 
INDFXl.O 
ICBCP=n 
I DL  B = 0 
CONTINUE 
DO  104  1=1 ,NMM 

INrASF ( I )=M*J 
DO  2 CO  1 J=1 . M 
ISTATI J»=R. 

DO  2002  J = MP I ,N 
ISTATI J > =1 
DO  IC  IT=1,NMMP! 
DO  i?  L=1 • SMMP1 
B1 INVI L . I I I = C . 
B I INVI 11,11)  = 1 

DO  TO  1=1, NMM 

xri  1 1 > = r . 

XBHNMMPI)  = t. 

T OL  R I = 1 ,CC-I0 


-C 

C 

35c 

. 2800 


start  thf  simplex  algorithm 
SOLVE  THE  OUAL  PHORLEM 
THE  NuMBEP  OF  VARIABLES  IS 
FOR  A TOTAL  OF  N VARIABLES 


M RFAL 


NMM  SLACKS 


35C  CONTINUE 
. 2800  DO  23  J - 1 , N 

RATS  = 0 • 

IF  IISTATI J) ,E0. t ) GO  To  52800 
IF  (J.C.T.M)  GO  TO  22 

RATS  - -R1  INV(  I .HEAD! j)Tl  ) FB  l In V ( 1 , T A IL ( J I ♦ 1) 
GO  TO  5? SO 

22  RATS  =-Bl  INVI  1 . J-MM  ) 

52800  RFDCOSI J )=  PATS 

23  CONTINUF 
2280C  CONTINUE 

l RMAX= l 

RMAX=RF OCOSI  1 ) 

DO  24  J =2  , N 

IFIRFDCOSIJ)  ,LF,  RMAX)  GO  TO  ? 4 

, RMAX=REDCnS(  J > 

I RMAX= J 

24  CONTINUE 

IF  I RMAX  .LF,  TOLRII  Go  to  AO  1 


CTIMfI  J) 


fOUOU  *"  * IP  « CL1 
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324  CONTINUF 

DO  26  L = 1 « NMMP I 
IF  1IRMAX.GT.M)  GO  T 3 50026 

V 1 < L ) =-Bl INV !L .TAIL!  I R MAX  )M)+B1I  NV  !L»HEAO{  IRMA  X ) ♦ I I 
GO  TO  26 

326  Yl(L)  = 01  INV!L , IRMAX-M+1  ) 

26  CONTINUF 

Yin)  = Yl(l)  - CTIMF(IRMAX) 

NUMREP=0 

on  27  l =2,NMMP1 

IF(Y1(L>  *LF.  TOLR1)  NUMBFR  = NUMBEP  «•  1 

I F( NUMBER  ,F0.  NMM)  GO  TO  403 
RMIN=.99D  20 
I R M I N=  0 • 

OO  32  II  =2 « NMMP 1 
IF(V1(  I l ).l  F.  T0LR1  ) GO  TO  32 
RATS  =XB1 ( I I J/Yl ( I l ) 
rr=rats-fmi n 

IF(PP  .GE.  0.00)  GO  TO  32 

RMIN=P4TS 

I RMIN=I I 

CONTINUF 

DO  33  J = 2 , NMMP 1 

WW=B1 INVIIRMIN  ( Jl/Yl  (1  RM I N > 

DO  37  L-1.NMMP1 

81  INV!L.J)=81  INV!L,J)-WW*Y1<L) 

51  I NV ( IRMIN  . J)  = WW 


UPDATE  THE  BASIC  VARIAOLFSJ  INBAS6  AND  XB1 


I ST  ATI  INFUSE!  IRMIN-1 > )=P 
I ST  ATI  IRMAX)=1 
I NBASEI I PM  IN- I )=IRMAX 
W=XBI(IRMIN  )/YI(IRMIN  3 
DO  34  1=1, NMMP1 


XBl  (I)  = XB1 (I  ) -Y1 ( I )*W 
XB11IRMIN  ) = W 
GO  TO  350 
•RITE! 6,570) 

FORMAT! 1H0.5X, *NO  FEASIOLF  SOLUTION  EXISTS* 

*.  • ) 

WRITE! 6.R50) 

FORMAT ( 1H1 ) 

GO  TO  9*P 


COPY  AVAILABLE  TO  DOC G DOES f NOT 
PERMIT  FULLY  LEGIBLE  PRODUCTION 


CHECK  YOUR  INPUT  DATA 


FNO  OF  THE  SIMPLEX  ALGORITHM 


CONTINUE 


KKK  = THE  NUMBER  OF  NOOES  ON  THE  CRITICAL  PATH 

KR!L>=  THF  L-TH  NODE  IN  THE  CRITICAL  PATH.  COUNTING  BACKWARDS 
FROM  THE  TERMINAL  NODE 

KKB=  THE  NUMBER  OF  ACTIVITIES  ON  THE  CRITICAL  PATH 
I BB ! L ) = THE  L-TH  ACTIVITY  ON  THE  CRITICAL  PATH,  COUNTING 
BACKWARDS  FROM  THF  TERMINAL  NODE 


IFI  ICBCP.EO.l ) GO  TO  6008 
IF!  I NO F XL • FQ  « 2 I GO  TO  3204 


1 80 
181 
182 

183 

184 

185 

186 
187 
1 88 

189 

190 
1 91 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 
20  7 
208 

209 

210 
211 


IN8ASF  IS  A SFT  OE  M INTEGER  VARIABLFS  WHICH  INDICATE  THF 


COMPOSITION  OF  T HE  CURRENT  BASIS.  FOR  EXAMPLE. 

INBASE(K)  = 7 IMPLIES  THAT  THE  K-TH  COLUMN  IN  THE  BASIS  B 
CORRESPONDS  TO  THE  7-TH  VARIABLE 


ISTAT  INDICATES  THE  BASIC  STATUS  OF  EACH  VARIABLE 

ISTAT(K)  = I IMPLIES  THAT  THE  K-TH  VARIABLE  IS  IN  THE 
DUAL  BASIS 

ISTAT(K)  = 0 IMPLIES  THAT  THE  K-TH  VARIABLE  IS  NOT  IN  THE 
DUAL  BASIS 


thf  following  statements  DETERMINE  THE  NODES  AND  ACTIVITIES  ON 
The  CRITICAL  PATH 


THF  DUAL  SOLUTION  IMPLIES  THE  FOLLOWING  OPTIMAL  SOLUTION  TO  THE 
PPlMAL  PERT  PROBLEM.  HOWEVER  SOME  OF  THE  NODE  TIMES! OTHER  THAN 
THF  LAST  ONF)  MAY  BF  HIGHFR  THAN  NECESSARY.  THUS  IN 
DETERMINING  THF  CRITICAL  PATH  AN  ALTERNATIVE  OPTIMAL  SOLUTION 
MAY  HAVE  TO  BF  IDENTIFIED. 
biinv  is  not  changed. 


DO  P 30  0 2 1=1. NMM 
X NOOF ( I 1=01  I NV ( 1 . IFI  ) 
KKK  = I 
KP( 1 )=NMM 
IK=KB(KKK) 


DETERMINE  whether  THF  time  TO  peach  NODE  I X 
AS  LARGE  as  INDICATED  From  THE  DUAL  SOLUTION 


is  necfssarily 


SMI N=9999P9 , 

I SMI N=0 

DO  8 3000  T = 1 . M 

IFIHEADI  I ) .NF.ItO  GO  TO  83000 

SLACK=X NOOF (HEAD! I ) ) - XN ODE < T A l L < I H-CTlMEC  I) 
IF( SLACK  .GE.sMI N » GO  TO  83000 
SMI N=SLACK 
I SMI N= I 
CONT INUF 

TFISMIN.LT.O.OOO! > GO  TO  83C03 


the  time  FOR  NODE  IK  WAS  UNNFCCE S S AR I L Y LARgF 


XNCDE<  IK)  = XNODE ( IKl-SMIN 

KKK=KKK-1 

GO  TO  8 3C0  1 

IBB ( KKK  ) = ! SMI N 

KKK=KKKM 

KBlKKK  >=TA!L<  ISMIN) 

|F| TAIL (ISMIN) .GT.l  > GO  TO  83001 
KKB=KKK-t 

IF< INOEXL.EO.I ) GO  TO  3121 
I PAP  M=  I PAP MF  1 
I F ( TPARM.GT .4 ) GO  TO  ?9 I 0 
IF( IPARM.E0.3 ) GO  TO  6*00 
IF < I PARM.FQ.A)  GO  TO  6*01 


ICRITP(L)=  THE  L -TH  ACTIVITY  ON  THF  OPIGINAL  CRITICAL  PATH 


c 

c 


84 


C *CPB=  THE  NUMBER  OF  ACTIVITIES  ON  THE  ORIGINAL  CRITICAL  PATH 

C 

TOTAL  = Ft  I I NV  ( 1 , NMMP 1 ) 

<CPe=K.<B 
I CPITNI 1 ) = NMM 

on  280?  i=i,kcpr 

I CR  I TN  (14-11  = KB!  IM  ) 

280?  I CR I TP ( [ )=ln8(  I ) 

x=total 

WRITE! 6.851 > X 

851  format  ( iH0  .sx  .•  the  critical  path  time  whfn  each  activity's  comple 

*T I ON  TIME  IS  SET  EOUAL  TO  ITS  AVERAGE  IS  = *.E15.5> 

WRITE! 6.7606)  KKK 

7606  FORMAT!  1H0.10X ,*  THE  «,t3,*  NODES  ON  THE  CRITICAL  PATH  ARE  AS  FOLLO 
*WS  BEGINNING  WITH  THE  TERMINAL  NODES*) 

WRITFI6.77C7)  (KB! I ) . I = I .KKK  I 
7707  FORMAT ( I5X , 20!  I 3.*  , • ) > 

WRITE! 6 .7710)  K KB 

7710  FORMAT  ( 1H0. 1 CX  .*  THE  M3,'  CRITICAL  ACTIVITIES  ARE  AS  FOLLOWS  BEGI 
♦ NNINC,  WITH  THE  TERMINAL  ACTIVITY!*) 

WRITF(6.7707)  ( IBI»(  I ) . I = 1 ,KKB  ) 

RE AD ( 5 . 2R?  0 ) THETA, LAMBDA 
2920  FORMAT ( 2F5. 2 ) 

WRI TE(6.3071 ) THETA.LAMBDA 

3071  FORMAT ( 1HC .1  OX ,« THETA  = *,El5.5>*  LAMBDA  = *,E15.5) 

C 

C SAMSI7  = THF  NUMBER  OF  ACTIVITY  TIME  CONFIGURATIONS  TO  BF 

c RANDOMLY  SELECTED  FOR  CONSIDERATION  IN  EACH  CLUSTER 

C 

c notf:  since  this  is  a random  samplf  . some  percentile 

c COMBINATIONS  may  be  considered  more  than  once. 

c 

BF AO  (5,32091  SAMSIZ 
3209  FORMAT  (110) 


c 

C the  COMPLETION  time  FOR  ALL  ACTIVITIES  IS  SFT  TO  THEIR  LOWER 

r PFPCENTILF.  THF  RESULTING  CRITICAL  PATH  T IMF  IS  A LOWER 


C BOUND  ON  THF  EXPFCTED  CRI 

c 

DO  6402  1=1 ,M 
640?  CTIMF!  I 1 - FLO!  I ) 

CALL  0 I NVA ( t?8C  0 ) 

64  CO  CPLO=  BIINV!  1 .NMMP1  ) 


ICAL  PATH  TIME. 

COPY  AVAILABLE  TO  DDC  DOES  NOT 
PEMUT  FULLY  LEGIBLE  PBODUCTION 


WRITF! 6.6405)  CPLh 

6405  FORMAT ( IHC  »5X. * A LOWER  BOUND  ON  THF  EXPECTED  CRITICAL  PATH  TIME  IS 
* * • ,F15.f>) 

WRITE! 6,7606)  KKK 

WRITE! 6 .7707 ) ( Kfl( I > . I = 1 .KKK  ) 

WRITF! 6.7710)  KKB 

WRITF ( 6,7707)  ( IBB!  I ) . I *! ,KKB  ) 


O 

c THE  COMPLETION  TIME  FOR  ALL  ACTIVITIES  IS  SET  TO  THEIR  UPPFR 

C PFPCENTILE.  THE  RESULTING  CRITICAL  PATH  TIME  IS  A UPPER 

-C  BOUND  LN  THE  EXPECTED  CRITICAL  PATH  TIME. 

C 

DO  6406  1 = 1, M 
6406  CTIME! I ) = PHI ( I ) 

CALL  B I NVA!  028C0  ) 

6401  CPUH=  31 INVI1  ,NMM»I  ) 

WRI  TF  ( 6 , 64 09  ) C.PUB 


300 

301 

302 

303 

304 

305 

306 

307 

30  <2 

309 

310 

31  1 

312 

313 

314 

315 

316 

317 
31  8 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 

339 

340 

341 

342 

343 

344 
34b 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 
35R 
359 


n n 
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6409 


FORMAT { 1HC . SX  , • A 
* = • ,F1 5.5) 
WRITF! 6.7616) 

WRI TF! 6,7707) 

WRI TFI6.7710) 

WRI TF ( 6 • 7707 ) 


UPPFR  BOUND  ON  THF  EXPECTED  CRITICAL  PATH  TIME  IS 


KKK 

(KB( I >.I  = 1 .KKK) 
KKR 

< IRB< I ) . 1=1 , KKR ) 


C 

C 

c 

c 


path  time 
BOUND 
TIMFS 


C 

C 

C 


FD(  I ) = THE  LOwFR  BOUND  ON  THE  EXPECTED  CRITICAL 
PLUS  1/IEDF  OF  THF  DISTANCE  TO  THE  UPPFR 
NLFFD(I)  = THF  OBSERVED  NUMBER  OF  CRITICAL  PATH 
THAT  ARE  < OR=  FD  t I ) 

FD  AND  NLFFD  ARF  USED  TO  BUILD  AN  'EMPIRICAL*  DISTRIBUTION  OF 
THE  CRITICAL  PATH  TIMES 


6412 


C- ( CPUS-CPL8) / IEDF 
DO  6412  K= 1 , KCPB 
DO  641?  1=1.1 EDF 
FD ( I ) = CPL13  + I*C 
NL  EFD<  I » = 0 


THE  ASSOCIATE  GROUPS  ARE  NOW  FORMED 


C 

3165 


2825 


2 801 


22825 


WRI TF< 6 , 31 65) 

FORMAT!  INI .6X.  * The  ASSOCIATES  ARE  NOW  I DE NT  I F I E D S • ) 
i r t 1 1=  i 
no  2825  1 = 1. M 
CTIMFI  I )=COT(  I ) 

IWWWQ=ICRITP( 1 ) 

CHANG=  LAMBDAWS IGMA! IWWWQ) 

TFx=COT ( TWWWO ) -CHANG 
IF(TEX.LT. 0.0 ) CHANG=COT ( IWWWO) 

CTIMFI I WWWQ) =COT| l WWWQ) -CHANG 
CALL  R I NVA ( 02  800 ) 

CONTINUE 

IF(  I ST ATI  I WWWQ ) .EO.  1 ) CALL  B I Nv 1 (622  825,COT(  IWWWQ) .CTIME( IWWWQ)  . 
♦IWWWO) 

REPCOSI I WWWQ)  = REDCQS! I WWWQ ) TCOT ( I WWWO ) -CT IME ( IwWWQ) 

Ct IMF! I WWWQ )=COT ( IWWWQ) 

I WWWO- I CRITP<  I I I I I ) 

CHANG=  LAMBOATSIGMA! IWWWQ  ) 

TFX  = COT  ( IWWWQ ) -ChANG 
IFITEX.LT. 0.0)  CHANG=COT( IWWWQ) 

C T I ME ( IWWWQ)=C0T< I WWWQ) -CHANG 

IF  ( ISTATI  IWWWQ)  .EQ.l  ) CALL  B I Nv 1 ( 62 2800  , CTI ME < I WWWQ ) .COT ( I W WWQ ) . 
♦I WWWQ) 

RFDCOSI  IWWWQ )=REDCOS(  I W W WQ ) -COT ( I WWWQ ) + CT I ME ( IWWWQ) 

GO  TO  22800 


C 

c 

c 

2910 


DFTERVINE  ASSDCIATF  GROUP 


291  3 


29 1 2 
29  5 1 


NI  NAG( I I I I I )=0 
DO  2911  K=l . KKR 
KK=  1 

IF < IBR (K  ) . EQ. ICRITP(KK)  ) GO  TO  2911 
IF! KK  .GE.KCPB ) GO  TO  29  12 
KK  =KK ♦ 1 
GO  TP  29 1 T 

NINAC.I  I I III  >*N|NAC.(  I III  I )♦  I 

ASS CMP I 1 I I 1 I . NINAG!  I I II  I ) )= IMB( K > 

CONTINUE 


360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 
3 75 

376 

377 

378 

379 

380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 

392 

393 
39  4 

39  5 
396 
3 97 

398 

399 

400 

401 

402 

403 

404 

40  5 

406 

407 
1 j 


’i 
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WRItF<  6,29 161  I III  I , ICR  I Tp<  I II  I I ) ,NI  NAG(  I II  II) 

2915  FORMAT  I 1H0, 10X ,• THE  NUMBER  OF  ASSOCIATES  ASSOCIATFO  WITH  THE  *,I3. 
* • — TH  CRITICAL  PATH  ACTIVITY,  I.E.  ACTIVITY  • « I 3 « • • IS  = *.I3) 

IDUCK=N  IN AG ( I ! 1 1 I ) 

IFI IDUCK.FO.O 1 GO  TO  2810 

WRITEI6.29161  ( ASSG«P< 11111.11.1  = 1 . I DUCK  1 

2916  FORMAT ( 1H0 . 15x .* THE  ACTIVITIES  IN  THF  ASSOCIATF  GROUP  ARF  AS  FOLLO 
*WS*  15X.50I 13.* .• 1 1 

281 C 11111=111114-1 

TFI  I I I t I .LF.KCPO)  GO  TO  2801 

c 

C DETERMINE  THE  CLUSTFRS 

C 
C 

c 
c 
c 
c 
c 
c 
c 
c 

C NCLINC  and  ClInCL  HELP  KEEP  TRACK  OF  WHICH  CLUSTER  THE 

C CRITICAL  PATH  ACTIVITIES  ARE  IN 

C 
C 

C RELflw  FORMS  CLUSTERS  BY  PUTTING  FACH  CRITICAL  PATH  ACTIVITY  IN 

c SE° A R AT E CLUSTER  AND  THEN  ADDING  EACH  CRITICAL  PATH  ACTIVITY*S 

C ASSOCIATES  TO  ITS  CLUSTER 

C 

NCLUS=KCPR 
DO  3r20  1 = 1 , KC  PR 
NCL  INC  ( I 1=  1 
CL  I NCL I 1,1  1=1 
NlKCl  ( I 1 =NtN'AG(  I 1*1 
I NCL US<  1,1  1 = IC»ITP(  I ) 

IFININAGI I 1.EQ.0  1 GO  To  30  20 
I DUC  K=  N I KcL ( I 1 
DO  3C21  J=  2,1 OUCK 
J J= J-l 

3021  INCLUSI  I , J >=ASSGRP(  I , J J I 
302C  CONTINUE 
C 

C BEl'lW  POOLS  CLUSTERS  FORMED  FROM  ASSOCIATES 

C 

I A = 0 

3031  !A*IAM 

IFI I A.GF.KCPP 1 GO  TO  3030 
I F (NCLUS.EO. 1 1 GO  TO  3030 
I D I A=N I NCL ( 1 A 1 
IFI IOI A.FO.Ol  GO  TO  3C 3 1 
I AA= I A+ 1 

DO  3023  I I = I A A , K C PR 
IDI I =NI NCLI I I 1 
IF  < ID!  I .EO. 01  GO  Tn  .1023 
DO  3f>26  1*1.  I DU 
DO  «<V‘>  Jil.IOII 

ire INCLLSI  II .Jl.ro. INCLUSI  |A, m GO  TO  3027 
3025  CONTINUE 

GC  TO  3021 


THE  CLUSTERS  ARE  POOLED  TOWARD  THE  TERMINAL  NODE 
NCLUS  = THE  NUMBER  OF  NCN-EMPTY  CLUSTFRS 

N I NCL (II  = THF  NUMPFR  OF  ACTIVITIES  IN  THE  T-TH  CLUSTER 
INCLUSI I, J)  = THE  J-TH  ACTIVITY  IN  THF  I-TH  CLUSTER 
NCLINC(I)  = THE  NUMBER  OF  CLUSTERS  COMPRISING  THE  I-TH 
CLUSTER  AFTER  POOLING 

CLINCH  I. J)  = THE  J-TH  CLUSTER  WHICH  HAS  BEEN  POOLED  INTO 
THE  I-TH  CLUSTER 


87 


* 

♦ 


3027  NCLUS=NCI.US- ! 

DO  3029  J= 1 • I D I I 
on  30  2 9 I =1  « I D l A 

IF(  INcLLSI  II. JI.FO. INCLUSI IA.I  ))  GO  TO  3028 
30?9  CONTINUE 

NI NCL ( I A ) = NI NCL ( I A) +1 
INCLUSI  I A, NT NCL ( I A ) >=INCLUS{  I l .J) 

3028  CCNTfNUF 

Ilfs™""  - WSS 

GO  TO  031 
3030  CONTTNUF 


c pflow  describes  clusters  after  pooling  based  on  the  associates 

r 

WRIT F( 6 .30  33  I NCLUS 

3r>  33  FORMAT  ( 1H1  ,10X  , • THFPF  ARE  ».I3.*  NflNF  MPT  Y CLUSTERS  AFTER  POOLING  O 
*N  ThF  basis  of  ASSOCIATES  ONLY.*) 

I T = <- 

00  3034  l^l.KCPB 
IF(NINCL(  I ) .EO.O)  GO  TO  3034 
11=114-1 

1 DUC  J=M  NCL  ( I I 

WPI TF( 6 «30  3s>  l,(  INCLUS  ( I . J>  . J = I . IDUCj  ) 

3035  FORMAT  ( 1H0 .10X,  *THF  Ar.TlVlTfES  IN  THF  *,I3.*-TH  CLUSTER  APE  AS  FOL 
♦LOWS  : I ./,|5X.S0(I3,','»I 
3034  CONTINUE 
C 

C DFSCRIBES  KHFRE  EACH  ACTIVITY  IS  BEFORF  EL  I MI  NANTS  <*RE 

0 CONSIDERED 

r 
r 

C EXAMINE  EACH  ACTIVITY  AND  DETERMINE  WHICH  CLUSTER. IF  ANY.  IT  IS 

C IN. 

C LFFtUI  = C IMPLIES  THAT  THE  I-TH  ACTIVITY  IS  NOT  IN  ANY 

C CLUSTER 

C lEFT(l)  = J IMPLIES  THE  I-TH  ACTIVITY  IS  IN  THE  J-TH 

C CLUSTER 

C 

WPITFI 8.3104) 

3 1 C 4 FORMAT(  lHC.tOX.* The  cluster  to  which  fach  activity  belongs:  •./.  I5X 
* • ( 7 FRO  IMPLIES  THAT  THE  ACTIVITY  IS  NOT  IN  ANY  CLUSTER)*) 

DO  3101  I = 1 , M 

LEFT! I ) ~C 

DO  3 1 C 2 J=l,KCPB 

I FI  NINCLI J ) .FO.C ) GO  TO  3102 

IDUCK=Nl NCL ( J ) 

DO  3110  K=1,IDUCK 

in  I.FQ.INCIUSI  J.K)  ) GO  TO  3107 
3 | 1 0 CONTINUE 
3 1 o?  CONTINUF 

GO  TO  not 
lie  7 IFFTIIIIJ 

3101  WRITFI6.3I03)  I.LEFT(I) 

3103  FORMATIIH  , 1 5X . • THF  »,I3»*-TH  ACTIVITY  IS  IN  fHE  *,t3,*-TH  CLUSTFR 
* • ) 

INDEXL  = t 


480 

481 

482 

483 

484 

485 

486 

487 

488 
4 89 

490 

491 

492 

493 

494 

495 

496 

497 
4 98 

499 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 

511 

512 
51  3 

514 

515 

516 

517 

518 

519 

520 

521 

522 

523 
‘24 
625 

526 

527 

528 

529 

530 

531 
53? 

533 

534 

535 

536 

537 

538 

539 


-1 1 

.X.fl 
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C 

r 

r 

r 

c 

r 

r 

r 


LEFTOVERS  ARF  ACTIVITIES  NOT  IN  CLUSTFPS  ATT FR  ASSOC1ATFS  HAVF 
O^FN  CQN^IOFRFf)  HUT  flFTOPF  FLIMINANTS  HAVF  RfcFN  CONSIDERFO 


OET  e R M l NF  THF  NUMRFp  qE  LEFTOVERS.  NLFFT 

IPPTOtL)  = J IMPLIES  THAT  THF  L-TH  LEFTOVER  IS  THF  J-TH 
ACTIVITY 


1122 

1121 


11?  1 


N1  FFT=  1 

00  IIP?  J=1,M 

1 F (l  FFT(  J ) .NF  ,r>  » r,(J  TO  .1122 
NLFFT^nL.F*  t«-1 

LrcTO(  NI.FF  T ) = J 
CONTI  nut 

WRtTF(6.11?1)  NLFFT 

FOOMftT  ( 1H0 , 1 r X , • ThFPF  AFF  *,11,1  ACTIVITIES  N )T  IN  ANV  CLUSTER  YC 
*T  . • ) 

WP I TF(  t . 132  1 ) 

FORMAT  ( 1H1  , sx .• THF  FLIMINANTS  OF  FACH  NON -CR I T I C AL- P A T H ACTIVITY  A 
*RF  NOW  DE TEPMI NFO: • 1 


r 

C 

c 

r 

c 

c 


FLIMINANTS  FOP  EACH  NON-C R I T I C AL -P AT H ACTIVITY  ARE  NOW 
OF  TE RMl NE  O 

NNNCP  3 THF  NUMiifr  OF  ACTIVITIES  NOT  ON  THE  CRITICAL  PATH 
MINCP(LF)  - THE  LE-TH  ACTIVITY  NOT  UN  THF  CRITICAL  RATH 


NNNICF  = M-K COI 

lt=<- 


or:  f-.'iCZ  U|,M 


SCO) 


S""? 


COPY  AVAILABLE  TO  DDQ  DOES  NOT 
PERMIT  FULLY  LEGIBLE  PRODUCTION 


5005 


Hi 


ACTIVITIES  NOT  ON  THE  CRITICAL  PA 


SOC  Ei 
5907 


J=  I 

IF  (I  ,E  0 . ICRI  TO  ( J ) ) GO  TO  5C0C 
J = JT  1 

IF  ( J .L c .K CRO  ) GO  TO  500  t 
LF=LF*  1 
N JNCPILF) -I 
CONTIMJF 

W& lTC< 5tSr 15 ) FNNCP 
FORMAT!  1H0 .CSX  ,• THFPE  ARF 
*TH.  thfy  ARF  AS  FOLLOWS:*) 

I F (nNNCP. FO.O  ) GO  TO  3124 
DO  FCC6  1=1, LE 
WPTTF ( 6 .5C97  ) I.NUNCP(I) 

FORMAT  ( 1 H .lEiX.I1,'.  * , I 3 ) 

TF< NNNCP. FO.f  ) GO  T0  3124 
LF=r 

LF=LP+i 

IF  ( IStATI  IWWWO)  .n.I  1 CALL  PINV1{6231?7,C0T<  IWWWQ)  .CTIMEI  IWWWO)  . 

* IWWWO) 

RFOCOS  (IWWWO)  = PLOCns(  Iwww0)-C.TIMF  ( IWWWQ)FCpT<  IWWWO) 

21127  CONTINUE 

CTIMEI IWWWO)  = COT( IWWWO) 

C T I Mr  I NONC P(  L E'  ) ) - CUT  ( NONCP(  LI))  «•  T HE  T A *S I GMA  ( NUNCP  < L E ) ) 

IF  I |S T AT (NUNCIMI  F ) ) .EO.  ) ) CALI  M I NV 1 ( 67756.CT I Mt(  NONCPl LF  It* 

* PUT  ( NUNOMl  f ) ) ,MUNCU(l  f ) ) 

prOCPF.  I NFNCp(LL)  )=rFOC'><'(  NONr.PU  f n-COT  ( NOtyCPI  L F ) ) WCT  I ML  C NONCP(  LF  ) 
♦) 

7756  IWWWO  = NCNCP(LE) 

WRI  TF(ft,11  S2  > NCNCP(I.E)  ,CTIMF(N0NCP(1  F)  ) 

US?  FORMAT  I IHO  • ///.  SX,*THf  CUMPI.rTION  T 1 Mf  FOP  TUI  M ACTIVITY 


3126 


SaO 
F 4 1 
S4? 
541 

‘>44 

545 

546 
54  7 
S48 
54V 

550 

551 
55? 
551 

554 

555 

556 

557 

558 

559 

560 

561 

562 
56  3%  S 


564 

565 

566 

567 

568  f 

569 
57C 

571 

572 

573 

574 

575 

576 

577 

578 

579 

580 

581 

582 
581 

584 

585 

586 

587 

588 

589 
59C 

591 

592 

593 

594 
695 
596 
897 

698 

699 


* HAS  BEEN  CHANCED  TD  ‘.F15.fi> 
GD  TD  2?8C 3 
CONTINUE 


89 


3121 

c 

C 

c 

c 

c 

r 

r 

c 


3131 

3 130 
3133 


3135 
31  36 

C 

C 

c 

r 

c 


c 

c 

c 

sole 

c 

c 

c 

c 


c 

c 

r 


501  1 

r 

c 

c 

c 


nr t^pmi ne  the  El 

CP!  r IC AL  PATH 

NE  = 

FGRP(J) 


MlNANTS  OF  the  LE-TH  ACTIVITY  NOT  ON  THE 

THE  NUMRFP  OF  eltminants  for  thf  lf-th 
ACTIVITY  NOT  ON  THF  CRITICAL  PATH 
= thf  J-TH  ElIMINANT  For  THE  LE-TH  ACTIVITY 
NOT  ON  THE  CRITICAL  PATH 


NF=r 

on  3133  K-t.KCPR 

DO  3131  1=1, KK3 

IF( IBB! I ) .FO. ICRITR(K>>  GO  TO  3130 

CONT INUT 

NF=NF+1 

EGRPINE  >=ICRITPIK1 
CONTINiJF 

WPItF! 6,31 33)  NF,NONCP(lE> 

FOpMftT ( 1H0 . 10X  , • ThFRE  ARE  ELIMINANTS  CORRESPONDING  TO  ACTIV 

* 1 T Y •,  M> 

IFlNE.EO.O  ) GO  Tfl  3171 
on  3135  K=1 .NF 

WRITFI  6 . 31  36)  K,  NONCPlLE  ) .FGRP(K) 

FORMAT!  Ih  . 1 4 X • • THF  ‘.I3.‘-TH  ELIMINANT  CORRESPONDING  TO  ACTIVITY 
* ‘ . I 3. • IS  ACTIVITY  • . I 3) 

DFTFBMINE  WHETHFR  NONCP(LF)  is  an  ASSOCIATE 
JA  = 1 IF  NONCP(LE)  IS  AN  ASSOCIATE 
JA  = 2 IF  NONCP(LE)  IS  NOT  AN  ASSOCIATE 


K = NCNC  PIlEI 
JA=  1 

IF(LFFT(K  ) .EQ.O  ) JA=2 
IF(JA.EO.?>  GO  TO  50 1 C 
!T  = LEF  T(K  ) 

thf  IT-TH  cluster  is  FXPANDED  TO  INCLUDE  EL  I M( NANTS 

GO  TC  5011 
CONTINUE 

ITTTT  IS  THE  ACTIVITY  NUMBER  OF  THE  FIRST  FLlMtNANT 

IT  is  THE  CLUSTFR  To  WHICH  THE  FIRST  EL! Ml NANT  CURRENTLY  BELONG 

I TTT  =E  GRP ( I > 

IT=LFFT( ITTT) 

LEFT!  Nfj  NC  P ( L F ) ) = I T 

THE  IT-TH  CLUSTER  IS  EXPANDED  TO  INCLUDE  EL  I Mj  NANTS 

NINCL!  I T > = M NCL  ! IT  > M 
INCLuS!  tT, NINCL!  IT) ) = N0NCP!LE) 

IF  I NF.EC.I  ) GO  TO  3171 
DO  3172  J=JA,NE 


IU  IS  THE  ACTIVITY  NUMBFR  OF  THF  NEXT  ELIMINANT 

if  iu  is  in  cluster  k,  thfn  cluster  k is  pooled  into  cluster  it 


it 


600 

601 

602 

603 

604 

605 

606 

607 

608 

609 

610 
61  1 
612 

613 

614 

615 

616 

617 

618 

619 

620 
621 
622 

623 

624 

625 

626 

627 

628 

629 

630 

631 
532 

633 

634 

635 
6 36 

637 

638 

639 
64  0 


641 

642 

643 


644 

645 

646 

647 

648 

649 

650 

651 

652 

653 

654 

655 

656 

657 

658 
689 


1 

I 

j 


Mil 


I 


90 

IU=EGRP( J> 

K=LFFT ( IU> 

3182  IF(lT.EO.K)  GO  TO  3172 
NCLUS=  NCLUS-1 

I W=  NCL  I Nf  (<) 

DO  318  3 I A = 1 , IW 

LFFT(  ICRITP(CLINCL(K.  I A 1)1= IT 
NCLI NC ( IT  » = NCL INC< I T >♦  1 

3183  CLINCL (IT, NCL  I NCI  IT ) > = CL INCL( K,  I A l 
NCL I NC( K)sn 

IW=NI NCL  <K ) 

N I NCL ( K )= C 
DO  3184  I A=1 , I W 
LEFT < I NCLUSIK  , I A ) ) = I T 
N I NCL  ( I T ) = N I N CL  ( I T ) ♦ 1 
3 184  I NCL US ( IT,NINCL(  IT)  )=INCLUS(K, I A) 

3172  CONTINUE 
3171  CONTINUE 

IF(LE.LT.NNNCP)  GO  TO  3126 
C 

C END  OF  POOLING  PAGED  ON  FLIMINANTS  FXCEPT  FOR  THE  FOLLOWING 

C DESCRIPTION 

C 

WPITFI 6, 31 73  ) NCL  US 

3173  FORMAT ( 1H1  ,05X  , 'THERE  APE  •,I3,*  CLUSTERS.') 

1)0  3178  I = 1 « K C PB 

IF(NINCL(  I ) .EQ.O  » GO  TO  3176 
I DD=NI NCL ( I ) 

WRITE* 6,3174)  NINCL( I ). I .( INCLUSC I, J) ,J*1, IDD) 

3174  FOPMAT(  1H0.10X, •THERE  ARE  *.13.  • ACTIVITIES  IN  THE  •«I3.*-TH  CLUST 
*EP.  THEY  ARE  AS  TOLL OwS 1 • • / • 20 X * 50 < I 3« • • • I ) 

IDUCK=rNCLlNC(  I ) 

WRlTF(S,317S)  NCL I NC  ( I 1 , ( CL I N CL ( I . J > , J=  1 * IDuCK ) 

3175  FORMAT < 1H0 , 1 5X  , I 3 ,»  CLUSTERS  HAVE  BEEN  POOLED  TO  MAKF  THIS  CLUSTER 
*.  they  there  as  follmks:  • ,/,?CX,50(  13  . ’ t • ) ) 

3176  CONTINUE 

WdItE  (6,2776)  SAMSI  l 

2776  FORMAT  (///SX,  • THE  NUMBER  OF  PERCENTILE  COMBINATIONS  EXP 

•LIFITLY  CONSIDERED  IN  DETERMINING  THE  UPPER  BOUNDS  AND  LOWFR  BOUND 
*S*/,6X, • ON  THE  NETWORK  COMPLETION  TIMF  DISTRIBUTION  AND  THE  UPPER 
♦ROUNDS  ON  ITS  MOMENTS  TS  EQUAL  TO  *,I5> 

CALL  CLOCK (XRAN) 

lY'JTS  = KRAN 

write  (6,3238)  IYUTS 

3238  FORMAT  <1H0. SX, 'THE  INITIALIZATION  PARAMETER  FOR  THE  SAMPLING  IS  I 
*Y  = • . TIC) 

C 

C STATEMENT  NUMBER  3124  MAPKS  THE  END  OF  POOLING  CLUSTERS  RASED 

C ON  LEFTOVERS  ANO  FLIMINANTS 

r 

3124  CONTINUE 

c 

c 

. c THE  FINAL  CLUSTERS  HA VF  NOW  BEEN  DETERMINED 

C TmF  ?**N!NCL<I)  RUNS  ARE  NOW  AVERAGED  FOP  ALL  I WITH  NINCLIII>0 

C THE  IZ7(  ) ARE  USFD  TO  REpRFSFNT  ALl  OF  THF  2**NlNCL 

. C POSSIBILITIES. 

C THIS  IS  WHFBF  THE  BINARY  REPRESFNT AT  ION  IS  CONSTRUCTED. 

C 

6DC8  IF(  ICBCP.EO.  1 > ICBCP*? 


664 

(>65 

666 

667 

668 

669 

670 

671 
67  2 I I 

673 

674 

675 

676  i 

677 

678 

679 

680 
681 
682 

683 

684 


I NDEXL  = 2 

I P=C 
| P= I R+ 1 

on  ftr 3 1 I =t , l o 

MOMENT!  IP,  I ) = 0,0 
IF < IR.6T.KCPS > GO  TO  320fl 
lE(NlNCL< IR)  .EO.O)  GO  TO  3200 
IOUCK=NCL I NC  ( IP) 

DO  3410  1 = 1,1  DUCK 

K=f L I NCL<  IR, I ) 

CONTINUE 
L=  ICRI rp(K  ) 

I p = r 

N I R=  C 


NSAVE  = 


NUMBER  OF  PERCEnTILE  COMBINATIONS  IN  ThE  SAMPLE 
PERCENTAGE  HE  THE  TOTAL  NUMBER  OF  PFRCFNTlLE  COMBINATIONS 
Explicitly  consioepfo 

vector  CONTAINING  THE  LOWER  BOUNDS  ON  THE  NETWORK 
COMPLETION  TIME  DISTRIBUTION  TO  BE  AVERAGED  WITH  THF 

upper  rcunds  on  thf  network  to  yield  the  average  network 

COMPLETION  TIME  DISTRIBUTION. 

NIB  ASSOCIATED  WITH  NSAVE. 


!C=NINCL< IR) 

IB  = ?**NINCL t IR  ) 

IDOALL  = 1 MEANS  ALL  PERCENTILE  COMBINATIONS  ARE  EXPLICITLY 
CONS IPFPED. 

inOAII  = 0 MEANS  TO  SAMPLE. 


incALL  = n 

IF  (SAMSIZ ,LE . 0 .OP. SAMS  I Z.GE . IB > IOOALL=l 

DP  ????  1=1, M 

CTI ME<  I ) = COT  I I ) 

GO  TO  2 000 ^ 


STATEMENT  1204  IS  THE  RE-ENTRY  POINT  FROM  ThE  DUAL  SIMPLEX 
ALGORITHM  WHEN  A CLUSTER  AVERAGE  IS  BEING  COMPUTED 


CONTINUE 

IF  (IDLB.EC.l)  GO  TO  R90C 
00  *032  1=1,10 
MOMENTIIR.I)  = MCMENT ( I P , I > 
If  UC0CP.NE.2)  GO  TO  10001 
X=  PI  I N V ( 1 ,NMMP 1 ) 

X=X-I . OD-I 0 

I=C 

1=1  + 1 

IEtX.GT.EDU)  > GO  TO  6420 

NLEFDI  I ) = NLEFOI  I ) ♦ 1 

CONTINUF 

IP=IP+1 

NIB  = NIB+1 


B 1 INVI 1 .NMMP1 )**I 


GFNERATE  next  PERCEnTILE  COMBINATIONS  TO  BE  EXPLICITLY 
CONSIDERED. 


IF  I IDCALL.EQ.C  ) GO  To  20500 


r>  n 
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IF  (TP.GT.ID)  GO  TO  32C  7 78 

RAN?  AM  = ID  7811 

GO  TO  ? C5P  1 78L 

2C5C  0 IF  ( nIB.GT.SaMSIz)  go  TO  3207  783| 

IYUTS  = TVUTS  *539  784 

IF  (IYUTS)  6210.6211.6211  783 

6210  [YUTS  = I YUTS+2147483647'M  786 

621 1 XPAN  = I YUTS  787 

XRAN  = XRAN4. 465661 3F-9  788 

R AN?  AM  = XRAN4DFL0AT ( IB- I ) ♦ 1 789 

20501  CONTINUE  790 

791 

CONVERT  THF  RANDOM  NUMBFP.  RANSAM,  TO  A BINARY  NUMBER  TO  DEFINE  A 792 

C PERCENTILE  COMBINATION.  793 

C 794 

KRAN  = PANSAM  795 

DO  8505  1=1. TC  796 

I HALF  = KHAN/?  797 

12  s KHAN  - I H AL  F*2  798 

L = I NCLUS ( I P , I ) 799 

CtIMF(L)  = I2*FHT(l)-!Z*FL0(L)  ♦ FLO ( L ) 800 

8505  KRAN  = (HALF  801 

CALL  BINVAI&2BO0)  802 

3207  nIB  = NIB-1  803 

IF  (I0LB.FO.1)  GO  TO  9C 1 1 804 

DO  6030  1=1,10  805 

6030  MOMFNT(IR.l)  = MOMENT ( I R . I J/Nl B 806 

IF  (ICBCP.E0.2)  GO  To  6011  807 

GO  TO  3200  808 

3208  WRITE  (6,6362)  809 

6362  FQPMAT  ( I HI  ) 810 

PNTB  = DFLOAT( NIH*1C0 l/OFLOAT! IB)  811 

IF  ( SAMS  I 7 .LF  .0  ) WRITE(6,6045)  812 

IE  ( SA MSI 7.GT.C  ) WR I T F( 6 , 60 5 6 ) SAMSI  Z 813 

6o  56  FORMAT!  5 X , • THE  FOLLOWING  TABLF  WAS  COMPUTED  CONSIDERING  AT  MOST  814 

* * . I 8 . * PERCENTILE  COMBINATION?  .•)  815 

60  a ft  FORMAT ( 5x  • * THE  FOLLOWING  TAHlE  WAS  COMPUTED  CONSIDERING  ONLY  *,  816 

* 18.’  RERCeNTiLF  COMBINATIONS  OR  «.F6.2.*  PFRCENT  OF  ALL  COMB  I NA  817 

♦TICKS. *1  818 

6045  FORMAT  (5X,*  THE  FOLLOWING  TABLE  WAS  COMPUTED  CONSIDERING  ALL  PERC  819 

♦ ENT  ILF  COMBINATION?.*)  820 

DO  10^00  J = 1 , 1 0 821 

t TMAX=9  822 

TMAX=0 • 823 

DO  50? I 1=1 , KCPB  824 

IFININCL! I ) .FQ.O)  GO  TO  5021  825 

IF  ( MOMENT (I.J).LE.TMAX)  GO  TO  5021  826 

TMAX  = MCMFNT(I.J)  827 

1 TM A x= I 828 

5021  CONTINUE  829 

WR I TF  (6.5023)  J.  J. MOMENT! ITMAX .J)  830 

5023  FORMAT  !1H0,5X,*A  LOWER  BOUND,  T- ( • , 1 2 , • t THFT A,L AM8DA > , ON  THE  • 831 

♦ , I?  , * -T H MOMFNT  OF  THE  NETWORK  COMPLETION  TIME  « »,FI5.5>  832 

l 0000  CTNTINUE  833 

C 834| 

C RFC.  IN  THE  PROCEDURE  FOR  DETERMINING  UPPER  BOUNDS  ON  THE  MOMENTS  838 

c of  thf  network  completion  time  and  lowep  bounds  on  the  836 

C DISTRIBUTION  of  thf  COMPLETION  times.  83? 

c «38 

c 639 


IF(NINCL ( I ) .FQ.O ) f.n  TO  6C04  857 

IF<  I .FC  .NCI.  ) GO  TO  60C4  858 

K=MNCL(NCL>  859 

JJ=MNCL<T>  B80 

NI  NCL  ( I ) sf'  861 

DO  6PC  5 J=1 . J J 862 

K = K ♦ 1 863 

6c05  INCLljS(  NCL,K)  = INCLUS(  I . J)  864 

N! NCL< NCL ) = N INCL< NCL ) *J J 865 

6004  CONTINUE  866 

6002  NNCL-NI NCL (NCL ) 867 

C 868 

C 

C 870 

C tHF  UPPFR  BOUNDS,  T ♦ ( P , THET A , L AM0D A > * ARE  NOW  DETERMINED.  87| 

C 872 

C 873 

C FOR  THF  SAkF  OF  NUMERICAL  ACCURACY  THE  PE&T  PROBLEM  874 

c WITH  NEW  ACTIVITY  TIMFS  is  INITIALLY  SOLVED  FROM  SCRATCH  875 

C INSTEAD  OF  UPDATING  AN  OLD  SOLUTION.  AFTER  THIS  876 

c Reinitialization,  thf  pfmaining  critical  path  times  arf  877 

C DFTERMINED  HY  UPDATING  THIS  SOLUTION.  878 

r 879 

DO  6C  C 6 I = I.M  880 

C T I M F ( I > = FHI(  I ) 881 

.6006  COT ( I ) =FHl ( I ) 882 

DO  6007  J=1,MNCL  883 

I = I NCL  US  ( NCL.  , J ) 884 

6007  CtIME( I )=.5*(FLO< I )«-FHI(I  >)  885 

DO  7101  1=1,1 FDF  886 

7101  NL  EFD ( I ) = 0 887 

ICRCP=1  888 

GO  TP  6010  fl8» 

6011  WRITE  (6,6264)  890 

6264  FORMAT!///)  891 

PNI‘1  ■=  DFLnAT(Nl9*100)  / DFLOAT  (IQ)  892 

IF  ( IDOALL.EO.O  ) WR I TF ( 6 , 6C 46 ) NIB.PNIB  893 

IF  ( IOOALL.EO. 1 ) WRITE(6.6C45)  894 

DO  100C9  J=1 , 10  895 

.10009  WRITE  (6.6012)  J , J, MOMENT ( NCL . J » 896 

601?  format  uho.sx.'an  upper  bound.  t+< », 12. • jtheta. lambda),  on  the  • 897 

*.I2.*-TH  MOMENT  of  thf  NETWOPK  COMPLETION  TIME  * *»E15.5)  8)8 

00  69f  0 I = ? , I FDF  899 


U u 


11  = 1-1  900 1 

69  00  NLEFOII)  = NLEFOII)  6 NLEFOII I)  9011 


WRITE  (6.63621  902 

PNJB  = DFLOAT ( N I B*1 00 ) /DFLOAT  ( IB)  903 

IF  ( IDOALL .EG.C  ) WR I T E( 6 ,6046 ) NlB.PNIB  90* 

IF  ( IDOALL. EO. 1 ) WRIT£(6,6045>  905 

WPITF! 6,6423)  906 

6423  FORMAT  { 1H0  ,5x  . *A  LOWER  BOUND  ON  THE  NETWORK  COMPLFTION  T I MF  f>I  STR  907 

♦ibution:  f-< . ; thetajlambda) • ) 908 

K N I R = NIB  9C9 

on  6421  Jsl.IEOF  910 

NSAVE!  I )=NLFFO( I ) 911 

X=NLFFO(I>  912 

X=X/NI3  913 

6*21  WRI Te< 6.6422)  FD ( I ) . THE T A , L AMBO A, X 91* 

6422  FORMAT! 1 7x . «P- ( • ,E15.5, • 1 • ,E15.5,« J • .E1S.5,  •)  = «,£15.S)  915 

C 916 

C Ml  7 

C BEGIN  THE  PROCEDURE  FOR  DETERMINING  UPPER  BOUNDS  ON  THF  918 

C NETWORK  COMPLETION  TIME  0 1 ST  P I HUT  I ON  919 

C 920 

C FOR  THE  SAkF  OF  nUmEPICAL  ACCURACY  THE  PERT  PROBLEM  921 

C WITH  NEW  ACTIVITY  TIMES  IS  INITIALLY  SDLVFO  FROM  SCRATCH  922 

C TNSTFAO  OF  UPDATING  AN  OLD  SOLUTION.  AFTER  THIS  923 

C REINITIALIZATION,  THE  REMAINING  CRITICAL  PATH  TIMES  ARF  92* 

C DETERMINED  by  UPDATING  THIS  SOLUTION.  925 

DO  9C  0 6 1 = 1 ,M  926 

C.TIME!  I ) = FL0I  I ) 927 

9006  COT ( I ) =FLO ( I ) 928 

00  9007  J = 1 , NNCL  929 

1 = I NCLUS ( NCL , J ) 930 

9007  CtIME! I!=.5*(FlO< I )4FHI(I  ))  93l 

DO  91 0 i I=i, IFOF  932 

9101  NLFFDI I ) = 0 933 

ICBCP  = 1 93* 

IOLO  = 1 935 

GO  TO  6010  936 

9011  no  9990  1=2, IFOF  917 

I T = I - 1 938 

9990  Nt  EFD(  I ) = NLEFD(I)  ♦ NLFFDI  II)  939 

WRITE  (6,6264)  9*0 

PNIB  = DFLOAT (NIH*100 ) /DFLOAT  (IB)  9*1 

IF  ( IDDALL ,EO.o  ) WRITE(6,6C*6)  NlB.PNIB  9*2 

IF  ( IDOAIL ,FQ* 1 > WRITE! 6,60*5)  9*3 

WPTTF(6.«423)  9*4 

9*23  FORMAT UHO ,5X ,• AN  UPPER  BOUND  ON  THE  NETWORK  COMPLETION  TIME  DI 5TR  945 

♦IBUTIONI  F*t . ITHFTAJLAMRDA ) » ) 946 

00  942 1 1=1, IFOF  9*7 

X=  NL  FF  D ( I ) 9*8 

X=X/NI3  9*9 

9*21  wPl TE I 6 ,9*22 ) FD ( I ), THETA . LAMBDA . X 950 

9*22  FORMAT  ( 1 7X  . *F«.(  » ,E1  5.5.  • J • .E15.5.*  1 • .F1S.5.  • I = '.E15.5)  951 

952 

THF  APPROXIMATE  NETWORK  COMPLETION  TIME  DISTRIBUTION  953 

WW1TF  (6,6362)  95* 

C 955 

WPl TE I 6 , 9472  ) 956 

9472  FORMAT! IH0.5X.*AN  APPROXIMATE  NETWORK  COMPLETION  TIME  DISTRIBUTION  957 

!<5X,  »F(  . ; THET  * . L AMRDA  ) e .5  * I F*C  . 1 THFT  A .LAM50A  I ♦ F-!*JTh  958 

*ETA. LAMBDA)  >*,//)  9S9 


Tn  n >3  r»  ■»  :3  n Ton  non  nnn 


On  9471  1 = 1 . I EOF 

X = .540FL0AT  <NSAVE(  I I l/KNIB4-,5*DFLOAT<  NLEFD(  I > l/NIB 
9471  *R!  TCI*. 9471J  FO  ( 1 ) . THE  T A . LAMRD  A . X 

9471  FORMAT ( 1 7X  « • F< • .E15.5. • ; • .Fl5.5.» l • .E15.5. • I = ‘.EtS.Sl 
999  WRlTFI6,850) 

STOP 

* END 


«•*«> 
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96* 

966 
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SURPOUTlNF  BINVA<*) 

IMPLICIT  RFAL*8  (A-H.O-Z) 

COMMON  31  I NV,  RFOCOS ,CT IMF.XB1 . INBASE . IHEA D « I T A I L . NMMPI ,NMM,N, I ST AT 
COMMON  M.MP1 

O IMF NS  I ON  I STAT ( 100)  , IHFAOI 6C ) . ITAILI60) .XB1 (4  1) 

DIMFnS ICN  31 INV(41 .4 1 ) . INBASEI40 ) • CT I ME ( 1 00 ) .RFOCOS I 100 ) 

UPOATF  THF  FIRST  ROW  OF  BIINV  AFTER  CHANGING  CT I ME 

00  1 1=?, NMMPI 

91  I N v(  1 . I ) = 0.0 
DO  1 J=?, NMMPI 

1 9 1 I NV Il.I)  = 911 NV ( 1 . I ) 4-  91  INVI J, T)*CTImE( INRASEI  J-l) ) 

UPOATF  VALUE  OF  THE  OBJECTIVE  FUNCTION 
XB1  (1)  = R 1 I NV ( 1 .NMMPI  ) 

returni 

ENO 

SUBROUTINE  9INV1  ( *, TmNFW.TmOLO.ID) 

IMPLICIT  RRAL48  (A-M.O-Z) 

COMMON  R 1 I NV.PEDC.nS  ,CTI  MF.XBl  • InBASE.  IHEAD.  IT*  I L.  NMMPI  .NMM.N.  |STAT 
COMMON  M.MP1 

01 MF NS ICN  ISTAT (lrO ) . IHFAOI 60 ) , I T A I L ( 60 ) . XB I ( 4 1 ) 

OImFNSICN  31INV(41.41  ),  I NBASE ( 40 ) . CT  I ME  I l 00 ) . REOCOS { 1 00  I 

COMPUTE  THF  REOCOS  CORRESPONDING  TO  ONF  CmANGF  IN  CT  IMF 

DO  P 1=1, NMM 

2 !F<  INRASF<  I)  .FO.IO)  11=14-1 
OIFF  = TMNFW-TMOLO 

TMNEW  IS  THF  NFw  T 1 MF  ANO  TMOLD  IS  THE  OLD  TIME  CORRESPONDING 

to  the  single  change  in  ctime 

DO  1 K = 1 , M 

IF< !STAT(K I.EQ.l > GO  TO  1 

REOCns(K)  = RFOCOS  ( K)-0IFF*(BIt  NV  ( I I .IHFADI  K)4-l)  - 
* R 1 I NV  ( II,  I T AI L ( K ) 4-1  )) 

1 CONTINUE 

on  3 K = MP 1 » N 

IF  (ISTAT(K).EO.l I GO  TO  3 

RFOCOS  ( K ) = REOCOS(K)  - 0IFF4B1  INVl  I I.K-M4I  ) 

3 CONTINUF 

UPDATE  THE  FIRST  ROW  OF  BIINV  AFfER  CHANGING  CTIME 
DO  10  1=2, NMMPI 

BIINV  (1  , I )=91  I NV  ( 1 , I ) 4-0  I FF*R  l INV  < 1 1 , I ) 
ir  CONTINUE 


c 


UPOATF  VALUE  OF  THE  OBJECTIVE  FUNCTION 
XR1 ( 1 ) = Bit  NV ( 1 .NMMPI  ) 


967 

968 

969  ‘ 

970 

971 

972 

973 

974 

975 

976 

977 

978 

979 

980 

981 

982 

983 

984 

985 

986 

987  1 

988 

989 

990 

991 

992 

993 

994  ( 

995 


!.i 


■ 

• ' 


996 

997 

998  > 

999  J 
1COO 
1001 
1002 
1303  • 

1004 

1005 

1006 

1007 

1008 
1309 
1010 
1011 
1912 

1013 

1014 
1019 
1016 
101  7 
101 A 
1019 
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Original  Subnetwork  Analysis  Program:  Flowchart 
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APPENDIX  D 

Monte  Carlo  PERT  Simulation  Program 

The  Monte  Carlo  PERT  Si  jlation  Program  will  generate  a random 
sample  of  network  completion  times.  The  required  input  is 

(a)  an  acyclic  network  with  one  sink, 

(b)  parameters  describing  each  activity's  completion  time 
distribution,  and 

(c)  the  size  of  the  random  sample  to  be  generated. 

Currently,  the  program  generates  each  activity's  random  sample  of  completion 
times  from  a chi-square  distribution  with  3 degrees  of  freedom  that  has 
been  linearly  transformed  to  have  specified  15-th  and  85-th  percentilas. 

The  activity  time  distribution  can  be  easily  changed.  The  basic  output 
is  the  ordered  sample  of  random  network  completion  times  and  the  correspond- 
ing empirical  distribution  function.  The  critical  paths  associated  with 
the  sample  of  network  completion  times  are  not  determined. 

The  generation  of  the  random  sample  of  network  completion  times 
Involves  only  one  network  but  varying  values  of  the  individual  activity 
completion  times.  It  is  computationally  faster  to  find  the  network 
completion  time  for  a new  set  of  activity  completion  times  by  "updating" 
the  network  completion  time  for  a previous  set  of  activity  completion 
times  than  it  is  to  start  all  over  each  time.  Since  the  Simplex 
Algorithm  applied  to  the  dual  of  the  PERT  problem  is  ideally  suited 
for  this  type  of  "updating",  the  basic  computational  technique  for 
determining  the  network  completion  times  is  the  Simplex  Algorithm 
(see  e.g.,  G.  Hadley,  Linear  Programming) . 


A listing  of  the  Monte  Carlo  Simulation  Program  is  given  at  the 
end  of  this  appendix. 


Specific  Input  Instructions: 


Card  1.  Col.  1-3  : The  number  of  activities  in  the  network.  Format  (13). 

Col.  4-6  : The  number  of  nodes  in  the  network.  Format  (13). 

Card  2.  Col.  11-15:  The  number  of  random  network  completion  times  to 

be  generated,  Format  (15). 

Col.  21-25:  The  number  of  parameters  needed  for  generating  the 

random  activity  completion  times,  Format  (15) . 

For  each  activity  one  card  with: 


Col.  11-15: 
Col.  21-25: 
Col,  31-40: 

Col.  41-50: 


The  origin  node  of  the  activity.  Format  (15) 

The  terminal  node  of  the  activity.  Format  (15) 
Parameter  1.  The  15-th  percentile  of  the  activity's 
completion  time  distribution.  Format  (F10.4). 
Parameter  2.  The  85-th  percentile  of  the  activity's 
completion  time  distribution.  Format  (F10.4) 


The  nodes  should  be  numbered  1,  2,  ...,  n with  the  sink  being  number  n. 


Example: 

The  program's  input  and  output  are  illustrated  in  terms  of  the  network 


in  Figure  O-l. 
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SAMPLE  INPUT 


2 

2 


12 

336.27 

429.47 

3 

57.47 

89.96 

4 

57.47 

89.96 

5 

57.47 

89.96 

6 

57.47 

89.96 

7 

57.4'' 

89.96 

8 

68.  ‘ii 

107.95 

8 

68.  6 

107.95 

8 

68 . o 

107.95 

8 

68. *6 

107.95 

8 

68.  i6 

107.95 

a 

150.36 

193.49 

10 

333.96 

403.85 

9 

333.96 

403.85 

11 

355.10 

409.85 

18 

141.75 

221.90 

13 

672.36 

783.11 

14 

560.89 

660.00 

15 

560.89 

660.00 

16 

560.89 

660.00 

17 

542.80 

638.71 

18 

111.10 

173.92 

19 

256.03 

346.98 

19 

302.80 

400.67 

20 

311.95 

410.71 

21 

423.58 

530.74 

22 

315.35 

415.71 

20 

7.66 

11.99 

21 

16.77 

22.55 

22 

11.49 

17.99 

23 

39.54 

48.87 

26 

301.91 

400.86 

26 

767.09 

892.74 

27 

350.31 

460.06 

24 

382.32 

464.98 

25 

385.28 

461.54 

24 

11.49 

17.99 

23 

16.28 

23.00 

27 

7.66 

11.99 

25 

20.86 

28.29 

28 

810.17 

976.10 

29 

15.32 

23.99 

30 

15.32 

23.99 

31 

57.47 

89.96 

31 

49.81 

77.96 

31 

53.64 

83.96 

31 

88.12 

137.94 

32 

3.83 

6.00 

13 

49.81 

77.96 

33 

109.01 

152.63 

32 

745.50 

811.32 

32 

714.74 

799.83 

14 
16 

15 
17 
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SAMPLE  OUTPUT 
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N h*  S N K N i0 
* * * P 


OOOOOOOOCTOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


o o o o o o o 
o o o o o o o 

■»  O * 

<A  XVBOA^OOOO 


x o«^Nr'^Moaaot)OAntf\.4NoooN>«4Nrininr>4M9>HSON^>4<oNOO«MAh>»n0n99in<rooeo 

4 r>  «n  »A  t/\  lA  IT  -O  <6  -o  >0  -O  in  m f*>  m «P  I-  O 4 *0  * *Tt  O -♦  rg  ^ N »*  »«  *«  V* -*  (A  <0  4 O 4 »* 

at  in  .4  m n m •*  4 m <a  ia  ia  n in  in  * wi*  nwn  • ni*f* 


Ok 


-I 

2 

Z NNA>IA«h««tt0««OrnaA^A4^«99lO'4NO>4NA44ri>#A^AMA«^O»4>4*4NNAANN«<|Ah 
•C  -*  H-H--^--^HNNNNNNNNNNNNNNNNMNAnAAAA^A«AHHHH 


H* 

3 


O 

«c 

o 


H«NNIMNNn^A4hN«Ah 


O^O'^piNANNHNrOnNNANNNA^^f+A^-f^Ah-NAMOOOOf 
M nmh-^hh-NNNh  h^-nnnN  NNnnHHAAnnMAAflN 


*4444444444<AlAlAiAlA¥MI 


O 

4 


o n 


C 


C 

r 

C 

r 

C 

C 

c 

c 

r. 

r 

c 

c 

c 

r 

r 

c 

c 

c 


c 

1 00 


r 

c 

c 

c 

r. 

c 

c 

r 

• c 
c 
r 
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PROGRAM  LISTING 
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MONTE  CARLO  PERT 
IMPLICIT  RFAL*8  (A-H.O-Z) 

FOP  THF  SAKF  OF  identifying  THE  appropriate  DIMENSIONS.  LET 
M = THE  NUMBER  OF  ACTIVITIES  IN  THE  NETWORK 
NMM  = NUMBER  OF  NODES  IN  THE  NETWORK 
NMM°1  - NMM  1 
N = M 4-  NMM 

NTRIAL  = NUMBER  OF  SIMULATED  SFTS  OF  ACTIVITY 
COMPLETION  TIMES 

np arm  = number  of  parameters  needed  for  generating  random 

ACTIVITY  COMPLETION  TIMES  FROM  A PARTICULAR  DISTRIBUTION 


INTEGER  TAIL!  M),HFAD<  m!.PARM<M,NPaRM),TIMES(NTRjAL) 

ni me ns  r cn  m i nv ( nmmpi .nmmpi i ,ct imeini . xri inmmpi » .y i inmmpi i 
DIME ns  I CN  RPDCOS(N)  . I ST A T < N > . I NB AS E < N MM ) , CDF < 5 > , NSPEC < 20  I 


Currently  thf  dimensions  arf  set  for 

M— 60 

NMM=40 

NTPI AL=?00O 

COMMON  B 1 tNV.RmCOS.CTlMF.xBI.  I NB  ASF  . HE  AD  , T A I L , N MMP  I , NMM  , N • I ST  AT 
COMMON  M , MP 1 

INTEGER  TAIL! 6C ) ,HFAD<  6C  ) .TRIAL .NSPEC (PC) 

D I MF  NS  I CN  INBASE <40 ) . P AR M < 60  . S ) .CDF ( 5 ) , T I MFS < 20 00 ) 

D l MF  NS  I CN  nilNV(4I  ,41  ) .CTIMEt  100  ) 

D I MENS  I ON  XBI  ( 41  ) . Y I < 4]  I ,REDCOS( 109  > . I ST AT  C 10T ) 

DATA  NSPFCI l ) , NSPEC ( 5 I .NSPEC  <9  I , NS PEC  I I 3) .NSPFCI  I 7) /5*  • C NON* / 

DATA  NSPEC< 2)  . NSPEC< 6 I . NSPFCC  10  I .NSPEC! 14) .NSPEC( 18 )/5**F  SP*/ 
data  NS PEC < 1 I .NSPEC  < 7 ) , NSPEC  ( 1 I > .NSpEC  < 15 ) . NSPEC 1 19  I /5* • EC  I F • / 
data  NSPfC ( 4 ) .NSPEC <8 ) . NSPEC < l 2 I ,NSPEC<  1 6> .NSPFC<  20 1/5**  I ED)  •/ 
DATA  BLANKS/*  •/ 

M = THE  NUMRfR  OF  ACTIVITIES  IN  THE  NETWORK 
NMM  = THE  number  OE  NODES  IN  THE  PERT  NETWORK 
OFAD( 5, ICO ) M, NMM 
FORMAT (2TB) 

N=NMM+ M 
MP1 =M+ 1 
NMMpl=NMMFl 
TRIAL  = C 

simulation  variables 

NTRIAL  = NUMBER  OF  SIMULATED  SETS  OF  ACTIVITY 
COMPLFTIDN  TIMES 

NP ARM  = NUMBER  OF  PARAMFTFRS  NEEDED  FOR  GENERATING  RANDOM 
ACTIVITY  COMPLETION  TIMES  from  A PARTICULAR  DISTRIBUTION 
TIMES  = VECTOR  CONTAINING  THE  OPTIMUM  VA(.UF  FOR  FACH  TRIAL 
PAR M ( I • I ) = THE  LOWER  PERCFNTILF  POINT  FOR  THE  I-TH  ACTIVITY 
PAP  M I I , ° ) = THF  UPPER  PERCFNTILF  POINT  FOR  THE  I-TH  ACTIVITY 


8 
9 
1 0 
1 1 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 
3? 
33 
14 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 
5.3 
54 


■i 


■i 

1 


55 


RFAC  (5.2591)  NTRIAL. NPARM 

THF  ACTIVITIES  ARE  DESCRIBED 


IN  TERMS  OF  THFIR  NODES 


56 

57 

58 

59 


c 

c 


C II=THF  TAIL  NODF,  THE  ORIGIN  NODE 

C J J=  THE  HEAD  NODF,  THE  TERMINAL  NODE 


107 


DO  6 10  1=1  ,M 

RF  AD ( 5 , 2 50 1 1 I I , JJ,  (PARMl  l ,J) . J=1  .NPARM) 

2 501  FORMAT  ( 1 OX,  I 5 ,5X,  15  , 5X.5F1C.4) 

T A IL  < I ) = II 
61 0 HEAD ( I )=JJ 

NNPARM  s A*NPAPM 
DO  222  1=1 .NNPARM 

22?  NSprc<I)  = BLANKS 
WPl TT ( 6. ?700  ) 

? 700  FORMAT  ( 1H1  ,15X.  • INITIAL  INPUT*) 

WPI  TE(6,270l  ) 

27q 1 FORMAT ( 1H0, tOX, 'ACTIVITY  ORIGIN  TERMINAL  PARAMETER  I P 

♦ APAMETEF  2 PARAMETER  3 PAR  AMETF  R 4 PARAMETER  5*) 

WRITE  ( 6.250D)  ( NSPFC < II , I = 1 , 20 > 

250  Q TORMAT  ( A 1 X • 4 A A , 4 ( 1 X . 4A 4 ) ) 

DO  2704  1=1, M 

? 704  WRl  TE  ( 6 . 2702  ) I , T A I L<  I ) . HF  AD  I I I , ( P AP  M ( I , J ) , J=  1 , NPAR M ) 

2702  FOPMAT  ( IH  , 1 3 X , I 3.  5x , I 7 , 7X , I 3 , 1 X , 5(  7 X , F 1 0 . 4 > ) 

r 

DO  1 C 4 1 = 1, NMM 

104  l NRA^F  ( I ) = M-f  I 
DO  2001  J=I ,M 
2011  I ST  AT ( J ) = 0 . 

DO  200  2 J=MP1  , N 
2002  ISTAT(J1=1 

DO  in  I 1 = 1 .NMMP1 
Dn  l 2 L =1  . N MMP 1 
12  R 1 I NV  < l ,11)  = 0 . 

1"  B 1 INV(  11,111  = 1 . 

OO  3 P 1=1, NMM 
30  XR 1 ( 1 ) = 0 . 

XBl(NMMPl)  = 1. 

TOLR1  = 1 .CD-I  0 
C 

DO  55610  !=MP1,N 
5561 0 CT IMP ( I ) = 0 . 

350  CONTINUE 
C 

C GENERATE  A SFT  OF  ACTIVITY  COMPLETION  TIMFS 

C 

TRIAL  = TRIAL  + 1 

CAtL  PANTIM  (CTIME.PAPM.M) 

CALL  BINVA  (f,2800) 

C START  THE  SIMPLEX  ALGORITHM 

C SOLVF  T HF  DUAL  PROBLEM 

C THF  NUMBER  OF  VARIABLES  IS  M REAL  ♦ NMM  SLACKS 

C FOR  A TOTAL  OF  N VARIABLES 

r 

2800  DO  23  J = 1 , N 
RATS  = 0. 

IF  { ISTAT ( J)  .60. ! ) GO  TO  52800 
IF  (J.GT.M)  GO  TO  22 

RATS  =-Bl  INV(  1 .HEAD!  J )«•!  )6BI  INVl  1 ,TA  ILC  J»4I  ) ♦ CTlME(J) 

GO  TO  52800 

22  RATS  =-BI INV(  1 . J-M*l  ) 

52800  RE OCOs ( J ) = PATS 
?i  COKTINUF 


! 


f-  i 


27 


32 


37 

33 

C 

C 

c 


38 


*C3 

530 


850 

r 

c 

c 

*0  1 

c 

c 


NUMRFR=P 

OC  27  l -? . NMMP1 

I F I Y 1 (L)  ,Lr.  TOLR1 ) NUMBER =NU MBFR ♦ 1 

IFIN1JM3ER  .EO.  NMM)  GO  TO  40  3 

RMTr>=.qon  ?o 

I RM  j N=0 . 

on  32  I t=2.NMMPl 
IF(YI(tl).LF.  TOl R 1 > GO  TO  32 
RATS  = XR1 I I l ) /VI < t I ) 
rr=rats-pmin 

I F ( PR  ,GF.  0.00)  GO  TO  32 
PMIN^RaTS 

IRMINrM 
CONTINUF 
00  33  J=?,NMMP1 

WW=Pl INV( IRMf N ,J)/Yl!lRMlN  ) 

00  .37  L = 1 iKMMPI 

•3)  INV!  L . J)  =R1  INV(L,  J)  -W#*Y1  (L  ) 

F3 1 INV!  1RMIN  , J ) = WW 

UPOATF  THF  BASIC  VARIAPLFS:  INbaSE  AND  XB1 

1 ST AT C I NO ASF ( IRMTN- 1 ) ) = C 
I ST  AT  ! IRMAX )=  1 

I NR ASF  ( I RM IN-  1 )= IRMAX 
*=XRlllRMIN  )/Yl!lRMlN  ) 

00  38  1 = 1*  NMMP 1 
X R 1 ( I ) = X R 1 I I >-Y) 1 I > *W 
XBH1PMIN  )=W 
GO  TO  2800 
WRITE!  6,530  ) 

FORMAT! 1H0.5X. • NO  FEASIBLE  SOLUTION  EXISTS.  CHECK  YOUR  INPUT  DATA 

*.  • > 

• El TE( 6.850  > 

FORMAT ( 1 hi  1 
GO  TO  08R 

FNO  OF  THE  SIMPLFX  ALGORITHM 

TIMES!  TRIAD  = B 1 I NV!  I . NMMP  I) 

IF  ITS | AL.LT .NTRIAL ) GO  TO  350 

INBASE  IS  A SET  QF  NMM  INTEGER  VARIABLES  WHICH  INDICATE  THE 


109 


c 

COMPOSITION  of  the  current  basis,  for  example  > 

lao 

c 

I NBAS  £ ( K ) = 7 IMPLIES  THAT  THE  K-TH  COLUMN 

IN  THE  BASI S B 

1B1  ! 

c 

corresponds  to  the  7-th  variable 

182 

r 

183 

r 

184 

ISTAT  INDICATES  THE  BASIC  STATUS  OF  FACH  VARIABLE 

185 

c 

1ST  AT ( K 1 = 1 IMPLIES  THAT  THE  K-TH  VARIABLE 

IS 

IN  THF 

186 

r 

DUAL  BASIS 

187 

c 

ISTaT(K)  = C IMPLIES  THAT  THE  K-TH  VARIABLE 

IS 

NOT  IN  THE 

188 

c 

DUAL  BASIS 

189 

c 

190 

c 

191 

c 

192  S 

r 

193  * 

CPDER  THE  RANDOMLY  GEnERATEO  NFTWQRK  COMPLETION  TIMFS 

194 

r 

195 

3 00 C 

LIMIT  = NTRIAL 

196  , 

IPASS  = NTRIAL-1 

197 

ICHNC,  a 1 

198 

DO  400 1 J=t , IPASS 

199 

LIMIT  a LIMIT-1 

200 

IF  ( ICHNG. FO.O)  GO  TO  3060 

201  ’il 

ICHnG  = 0 

202 

DO  4CC  ? 1 = 1 .L I MI T 

203  ,j| 

IF  (TIMES! I > .LF. TIMES!  I 61 ))  Go  TO  400? 

204  $ 

Tf  MP  = TIMES!  141) 

20  5 '| 

TIMFS! I+l ) = TIMES! 1 ) 

206 1 i 

TIMFS! I ) = T F MP 

20  7 | 

ICHNG  = 1 

208  1 -5 

400? 

CONTINUE 

209  j 

4001 

CONT INUF 

210  , 

/- 

211  y 

c 

DFSCPIHE  the  drofrfd  nftkkork  completion  times 

212  N 

c 

213 

3060 

WRITE! 6.9662  ) 

214 

966? 

FORMAT ( 1H1 ) 

215  j 

WRITF16.3012 ) NTRIAL 

21  6 

30i  ? 

FORMAT  (IHo.IPX.'tHF  CRITICAL  PATH  TImES  For  • , 1 5.  • 

TRIALS  ARF  AS 

217 

* FOLLOWS:  ! T!  MF/nOSERVFD  PERCENT  I LE ) • ) 

218 

LINF  = (NTRIAL441/5 

219 

Do  3050  J= l « L I NE 

220 

I ? = ( 1- 11*641 

221  1 

T 3 = J * 5 

22  2 ji 

ICNT  = 0 

2 23  j 

IF  (NTRIAL-I241 .LT.b)  13=  NTRIAL 

224 

DO  3051  K= 1 2 « 13 

225 

ICNT  = ICNT61 

226 

3051 

COF(ICNT)  = DFLOAT! K ) /OFLOAT ( NTR I AL ) 

227 

WRITE  16.3011)  ( T I MES 11)*!=  12.13) 

228 

30  11 

FORM  AT ( l HO  «5X,5! 6X.F15.5)) 

229 

WRITE  (6,3013)  ( CDF  I K ) » K=1 .ICNT) 

230 

3013 

FORMAT  (6X.5I6X.F15. 5) ) 

211 

3 050 

CONT INUF 

2 32 

999 

STOP 

233 

FND 

2 3* 

SUBROUTINF  BINVA!*) 

230 

IMPLICIT  RF  AL  WH  (4-H.II-Z1 

236 

COMMON  4 1 I NV.RLDCOS.CT I Mf  .XBl . InBASE, IHtAD. ITaIL.NMMPI .NMM.N,  I STAT 

237 

COMMON  M . MP  1 

238 

01 Mf  NS  I CN  I ST  AT ( ICO ) . I HE  AD! 6C),JTA1L<60).XB!(4I) 

239  l| 

D I MENS  I ON  J 1 I N V ( 4 1 * 4 1 ) , I N8  ASE  I 4 0 I , C T T ME  ( | 00  » .REDCOSI  100> 

UPDATE  THE  first  row  of  bjinv  after  changing  ctime 

DO  1 I=2,NMMP1 
H 1 INV< 1 . I)  = 0 .0 
DO  1 J=2,NMMP1 

I RllNVU.t)  = RlINV(l.I)  ♦ BI  INV(  J.II  *CTIME(  INBASEIJ-l  )> 

UPDATE  VALUE  OF  THE  OBJECTIVE  FUNCTION 

XBKII  = Bl  INV(  I tNMMPl  ) 

RETURN 1 
END 

SUBROUTINE  RANTIM  I CT I ME , p ARM, M ) 

SUBROUTINE  RANTIM  GENERATES  M RANDOM  TIMES  FROM  A SPECIFIED 
DISTRIBUTION  WITH  PARAMETERS  CONTAINED  IN  PAPMI60.5)  AND  RETURNS 
WITH  THE  RESULTS  IN  CTIMF<99>. 

IMPLICIT  REAL  * 8 <A-H,n-Z) 

DATA  J/0/. IY/1 9**7/ ,TPI /6. 2831 653/ 

DI MENS  I CN  PAR M ( 60 , 5 I , CT I ME (99  > .SAVT I M( 99 1 
IF  (J.NE.O)  C.O  TO  30 

the  FOLLOWING  generates  A CHI  SOUARE  RANDOM  DEVIATE  WITH  3 OF*  S 
TRANSFORMED  TO  MAKF  THE  LOWER  POINT  CORRESPOND  TO  THE  15-TH 
PERCFNTlLF  AND  THE  UPPFR  POINT  CORRESPOND  TO  ThE  85-TH  PERC. 

PARM(3,I>  = THE  PERCFNTlLF  DIFFERENCE 

DO  5 I =1 ,M 

PARM(I,3)  = PARM ( 1 i ? ) — P ARM (1,1) 

IF  (PARMII ,3) ,EO,0)  CTIMF(I)  = PARMII.I) 

J * JM 

IF  (MODI J,?)  ,F0. 1 ) GO  TO  20 
DO  IS  I=|,  M 

IF  (PARMII,  T)«FO*0  ) GO  TO  15 
CT  T ME(  I I = SAVT IM< I ) 

CONTINUE 

RFTURN 


GENERATE  A PAIR  QF  COMPLETION  TIMES  FOR  EACH  ACTIVITY. 

ThE  ROX-MULLFR  METHOD  IS  USED  TO  GFNEPATE  A PAIR  OF 
NORMAL  DEvIATFS 

UI  AND  U2  = UNIFORM  RANDOM  NUMBERS 
W*ns I N(  AN ) = STANDARD  NORMAL  RANDOM  VARIABLE 
W*DCQS(ANI  = STANDARD  NORMAL  RANDOM  VARIAEUE 
CHIT  = A CHI- SOUARE  RANDOM  VARIABLF  — GENERATED  USING  T HF 
METHOD  OF  WILSON  AND  HILFFRTY  - PROC.  NAT.  ACADEMY  OE 
SCIENCE,  1931 

Cl  AND  C2  TRANS=ORM  THE  CHI-SOUAPE  WITH  3 D.F. 

on  ioc?  i = i, m 

IF  ( PARM 1 1, 3). E 0,0)  GO  TO  300  2 
IV  = I Y *65539 
IF  I | Y I 3015.  33  I ft.  .30  I ft 
15  IV  = IY  ♦ 2 1 * 7*8.15*  7 ♦ I 


2*0 
2*1 
2*2 
2*3 
2** 
2*  si 
2*6 
2*7 
2*8 
2*9 

250 

251 

252 

253  | 
25*  ' 

255 

256 

257 

258 

259 

260 
261 
262 
263 

26*  j 

265 

266 

267 

268 

269 

270 

271 

272 

273 
27* 

275 

276 

277 

278 

279 

280 
281 

282  J 

283 
28* 

285 

286 

287 

288 

289 

290 

291 
29? 
293 
29* 

295 

296 

297 

298 

299 


3016  YFt  =IY 


U1  = YFL*.4t56  13F-9 
IY  = ! Y * 6 5 5 3 9 
IF  (IY)  302“).  $*2*,  3026 
3C?5  I Y t-  I Y ♦ 21474P3647  ♦ 1 

30’6  YFL=  IY 

U 2 = YFL*. 46561 3F-9 
W --?.*DLr)G(U?) 

W = r>SQRT(Wl 

AN  = T9I*ul 

C P - 4,51  'J27/°ANM  (1*3) 
ri  = ,79777-PAPM 1 I , 11 *C? 

CTTmF(I)  = ( ( . 3<5?530614*W*DSIN(  AN)  Al  , 335416269)  **3-Cl  1/C2 
SAVTIM1  I ) = ((.  3«?53061  4*W*0C0S(  AN)  ♦ 1 . 1354  1 626.9)  443-C1  1 /C2 
300?  CONTINIJP 
WFTURN 
END 


BASIC  DISTRIBUTION  LIST  FOR 
UNCLASSIFIED  TECHNICAL  REPORTS 


Operations  Research  Copies 

Office  of  Naval  Research 
(Code  434) 

Arlington,  Virginia  22217  3 

Director,  Naval  Research 
Laboratory 

Attn:  Library,  Code  2029  (ONRL)  • 

Washington,  D.C.  20390  2 

Defense  Documentation  Center 
Cameron  Station 

Alexandria,  Virginia  22314  12 

Defense  Logistics  Studies 
Information  Exchange 
Army  Logistics  Management  Center 
Fort  Lee,  Virginia  23801  1 

Technical  Information  Center 
Naval  Research  Laboratory 
Washington,  D.C.  20390  6 


Office  of  Naval  Research 

New  York  Area  Office 

207  West  24th  Street 

New  York,  New  York  10011 

Attn:  Dr.  J.  Laderman  1 

Director,  Office  of  Naval  Research 
Branch  Office 
1030  East  Green  Street 
Pasadena,  California  91101 
Attn:  Dr.  A.R.  Laufer  1 

Lt.  Dennis  Johnston  - ED3 
Manned  Spacecraft  Center 
NASA 

Houston,  Texas  77058  1 

l 

Mr.  William  Dejka 
Code  360 

NaVal  Electronics  Laboratory  Center 
San  Diego,  California  92132  1 

Professor  Lucien  A.  Schroit,  Jr. 

Root*  6731,  Boclter  Hall 

School  of  Engrng.  & Applied  Sci. 

University  of  California 

Los  Angeles,  Californio  90024  1 


Director  Copies 

Office  of  Naval  Research 
Branch  Office 
495  Summer  Street 
Boston,  Massachusetts  02210 
Attn;  Dr.  A.L.  Powell  1 

Director 

Office  of  Naval  Research 
Branch  Office 
219  South  Dearborn  Street 
Chicago,  Illinois  60604 
Attn:  Dr.  A.R.  Dawe  1 

Office  of  Naval  Research 

San  Francisco  Area  Office 

760  Market  St.  - Room  447 

San  Francisco,  California  94103  1 

Technical  Library 

Naval  Ordnance  Station 

Indian  Head,  Maryland  20640  1 

Bureau  of  Naval  Personnel 
Navy  Department 
Technical  Library 
Washington,  D.C.  20370 


STOP  82  2 

Library,  Code  0212 

Naval  Postgraduate  School 

Monterey,  California  93940  1 

Naval  Ordnance  Station 

Louisville,  Kentucky  40214  1 

Library 

Naval  Electronics  Laboratory 
Center 

San  Diego,  California  92152  1 

Naval  Ship  Engineering  Center 
Philadelphia 

Division  Technical  Library 
Philadelphia,  Pennsylvania  19112  1 

Dr.  A.L.  Slafkosky 
Scientific  Advisor 
Commandant  Marine  Corps  (Cods  AX) 
Washington,  D.C.  20380  1 


Purdue  University  Copies 

Graduate  School  of  Ind.  Admin. 
Lafayette,  Indiana  47907 
Attn:  Prof.  Andrew  Whinston  1 

Department  of  Economics 
Tisch  Hall  - 5th  Floor 
Washington  Square,  NY  Univ. 

New  York  City,  NY  10003 

Attn:  Prof.  0.  Morgenstem  1 

Dept,  of  Statistics 
Syntex  Research 
3401  Hillview 

Palo  Alto,  CA  94304  x 

Attn:  Prof.  Stuart  Bessler 

University  of  California 
Department  of  Engineering 
Los  Angeles , California  90024 
Attn:  I'rof.  R.R.  O'Neill  1 

Maritime  Transportation  Research 
Board 

National  Academy  of  Sciences 

2101  Constitution  Avenue 

Washington,  D.C. '20418 

Attn:  KADM  J.B.  Oren  USCG  (RET)  1 

Stanford  University 

Department  of  Operations  Research 

Stanford,  California  94305 

Attn:  Prof.  A.F.  Veinott,  Jr.  1 

Texas  A&M  Foundation 

College  Station,  Texas  77843 

Attn:  Prof.  H.O.  Hartley  1 

Case  Western  Reserve  University 

Cleveland,  Ohio  44106 

Attn:  Piof.  B.V.  Dean  1 

University  of  California 
Center  for  Research  in 
Management  Science 
Berkeley,  California  94720 
Attn:  Prof.  W.I.  Zangwlll  1 

U.S.  Naval  Postgraduate  School 
Department  of  Operations 
Research  and  Economics 
Monterey,  California  93940 
Attn:  Prof.  C.R.  Jones 


Dr.  Del  Cilmer  Copies 

Naval  Weapons  Center  (Code  607) 

China  Lake,  California  93555  1 


Dr.  Paul  Murrill 
Project  Themis 

Department  of  Chemical  Engineering 

Louisiana  State  University 

Baton  Rouge,  Louisiana  70803  1 

Director 

Army  Materials  & Mechanics  Research 
Center 

Attn:  Mr.  J.  Bluhm 

Watertown,  Massachusetts  02172  1 

Commanding  Officer 
U.S.  Army  Ballistic  Research 
Laboratories 

Attn:  Dr.  John  H.  Giese 

Aberdeen  Proving  Ground, 

Maryland  21005  1 

University  of  Iowa 
Department  of  Mechanics 
and  Hydraulics 
Iowu  City,  Iowa  52240 
Attn:  Prof.  W.F.  Ames  1 

Office  of  Naval  Research 
Resident  Representative 
1105  Guadalupe 
Lowlch'  Building 
Austin,  Texas  78701 
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Brooklyn,  New  York  11251 

Naval  Undersea  Warfare  Center 
3202  E.  Foothill  Boulevard 
Pasadena,  California  91107 
Technical  Library 

University  of  Chicago 
Statistical  Research  Center 
Chicago,  Illinois  60637 
Attn:  Prof.  W.  Kruskal 

Stanford  University 
Department  of  Statistics 
Stanford,  California  94305 
Attn:  Prof.  G.J.  Lieberman 

Florida  State  University 
Department  of  Statistics 
• Tallahassee,  Florida  32306 

Atm:  Prof.  I.R.  Savage 

Princeton  University 
Department  of  Mathematics 
Princeton,  New  Jersey  08540 
Attn:  Prof.  J.W.  Tukey 

Stanford  University 
Department  of  Statistics 
Stanford,  California  94305 
' Attn:  Prof.  T.W.  Anderson 

j University  of  California 

Department  of  Statistics 
Berkeley,  California  94720 
Attn:  Prof.  P.J.  Bickel 

University  of  Washington 
Department  of  Mathematics 
Seattle,  Washington  98105 
Attn:  Prof.  Z.W.  Birnbaum 

• 

Cornel  University 
S',,.  Department  of  Mathematics 

* White  IWU 

Ithaca,  N'w  York  14850 
j . Attn:  Proi  j.  Kiefer 


Department  of  Statistics 

Cambridge,  Massachusetts  02139 

Attn:  Prof.  W.G.  Cochran  1 

1 

Columbia  University 
Department  of  Industrial 
Engineering 

New  York,  New  York  10027 
1 Attn:  Prof.  C.  Derman  1 

New  York  University 
Institute  of  Mathematical 
Sciences 

1 New  York,  New  York  10453 

Attn:  Prof.  W.M.  Hirsch  1 

University  of  North  Carolina 
Statistics  Department 
1 Chapel  Hill,  North  Carolina  27515 

Attn:  Prof.  W.L.  Smith  and 

.Prof.  M.R.  Leadbetter  1 

Purdue  University 

1 Division  of  Mathematical  Sciences 

Lafayette,  Indiana  47079 
Attn:  Prof.  H.  Kubin  1 

University  of  California,  San  Diego 
1 Department  of  Mathematics 

P.0.  Box  109 

La  Jolla,  California  92038 

Attn:  Prof.  M.  Rosenblatt  1 

1 Florida  State  University 

Department  of  Statistics 
Tallahassee,  Florida  32306 
Attn:  Prof.  R.A.  Bradley  1 

1 New  York  University 

Deportment  of  Industrial 
Engineering  & Operations 
Research 

Bronx,  New  York  10453 

1 Attn:  Prof.  J.H.K.  Kao  1 

University  of  Wisconsin 
Department  of  Statistics 
Madison,  Wisconsin  53706 
Attn:  Prof.  C.F..P.  Box 


1 


1 


^ it  'iiwnr 


»•*. 


|;s 


3 

I : f 


A 


I 


Logistics  Research  Project  Copies 

The  George  Washington  University 
707  - 22nd  Street,  N.W. 

Washington,  D.  C.  20037 

/‘tin:  Dr.  W.  H.  Marlow  1 

International  Business  Machines 
Corporation 

P.0.  Box  218,  Lamb  Estate 
Yorktown  Heights,  New  York  10398 
Attn:  Dr.  Alan  Hoffman  I 

University  of  California 
Management  Sciences  Research  y 
Project 

Los  Angeles,  California  90024 
Attn:  Dr.  J.R.  Jackson  I 

Harvard  University 
Department  of  Economics 
Cambridge,  Massachusetts  02138 
Attn:  Prof.  K.J.  Arrow  1 

Cowles  Commission  for  Research 
in  Economics 
Yale  University 
New  Haven,  Connecticut  06520 
Attn:  "’rof.  Martin  Shubik  1 

Carnegie-Mellon  University 
Graduate  School  of  Industrial 
Administration 

Pittsburgh,  Pennsylvania  15213 
Attn:  Prof.  G.  Thompson  1 

University  of  California 
Department  of  Economics 
Berkeley,  California  94720 
Attn:  Prof.  R.  Rudner  1 

University  of  California 
Operations  Research  Center 
Institute  of  Engineering  Research 
Berkeley,  California  94720 
Attn:  Prof.  D.  Gale  1 

University  of  California 
Graduate  School  of  Business 
Administration 

Los  Angeles,  California  90024 

Attn:  Prof.  J.  Mnrschak  1 


The  University  of  Michigan 
Department  of  Mathematics,  H.K. 
Ann  Arbor,  Michigan  48104 
Attn:  Prof.  R.M.  Thrall 

Princeton  University 
Department  of  Mathematics 
Princeton,  New  Jersey  08540 
Attn:  ”rof.  A.W.  Tucker 

Case  Western  Reserve  University 
Systems  Research  Center 
Cleveland,  Ohio  44106 
Attn:  Prof.  M.  Mesarovic 

University  of  Texas 
Department  of  Mathematics 
Austin,  Texas  78712 
Attn:  Dr.  A.  Charnes 

Stanford  University 
Department  of  Operations  Research 
Stanford,  California  94305 
Attn:  Dr.  D.L.  Iglehart 

University  of  Delaware 
Department  of  Mathematics 
Newark,  Delaware  19711 
Attn:  °rof.  R.  Kcmage,  Jr. 

Stanford  University 
Department  of  Operations  Research 
Stanford,  California  94505 
Attn:  Prof.  F.S.  Hillier 

Dr,  Claude- Alain  Burdct 
A3st.  Prof.  Industrial  Admin. 
Carnogic-Mcllon  University 
Pittsburgh,  Pennsylvania  15?13 

Stanford  University 
Department  of  Operations  Research 
Stanford,  California  94305 
Attn:  Prof.  G.  B.  Dantzig 

Chief  of  Naval  Research  (Code  436) 
Department  of  the  Navy 
Arlington,  Va.  22217 

Science  Librarian 
Kresgc  Library 
Oakland  University 
Rochester,  Michigan  48043 


Copies 


, y University  of  Connecticut 
Department  of  statistics 
Storrs,  Connecticut  06268 
Attn:  Prof.  H.O.  Posten 

Prof.  Emanuel  Parzen 
Statistics  Department 
State  Univ.  New  York  at  Buffalo 
4230  Ridge  Lea  Road 
Amherst,  New  York  14226 

ARCON  Corporation 
Lakeside  Office  Park 
North  Avenue  at  Route  128 
Wakefield,  Massachusetts  01880 
Attn:  Dr.  A.  Albert 

Stanford  University 
Department  of  Statistics 
• . Stanford,  California  94305 

Attn:  Prof.  H.  Chernoff 

. . Yale  University 

Department  uf  Statistics 
New  Haven,  Connecticut  06520 
Attn:  Prof.  L.A.  Savage 

Rutgers -The  State  University 
Statistics  Center 
New  Brunswick,  New  Jersey  08903 
Attn:  Prof.  H.F.  Dodge 

Yale  University 
. Department  of  Statistics 

New  Haven,  Connecticut  06520 
Attn:  Prof.  F.  J.  Anscombe 

*v> 

Purdue  University 

Division  of  Mathematical  Sciences 

Lafayette,  Indiana  47907 

Attn:  Prof.  S.S.  Gupta 

Cornell  University 
Department  of  Industrial 
Engineering 
Ithaca,  New  York  19850 
Attn:  Prof.  R.E.  Becbhofer 

» 

Mrs.  Barbara  Eaudi 
Univ.  Program  Coordinator,  B.E. 
NASA  Johnson  Space  Center 
Houston,  TX  77058 


Copies  Professor  Geoffrey  S.  Watson, 

Chairman 

Department  of  Statistics 
1 Princeton  University 

Fine  Hall 

Princeton,  N.  J.  08540 

Stanford  University 
Department  of  Mathematics 
1 Stanford,  California  94305 

Attn:  Prof.  S.  Karlin  * 

University  of  Sheffield 
Department  of  Probability 
and  Statistics 

1 Sheffield  10,  ENGLAND 

Attn:  Prof.  J.  Ganl 

University  of  California 
Operations  Research  Center 
1 Institute  of  Engineering  Research 

Berkeley,  California  94720 
Attn:  Prof.  R.E.  Barlow 

Stanford  University 
1 Department  of  Statistics 

Stanford,  California  94305 
Attn:  Prof.  H.  Solomon 

Applied  Mathematics  Laboratory 
1 Naval  Ships  Research  Development 

Center 

Washington,  D.C.  20007 

Systems  Analysis  Division 
1 Room  BE760,  Pentagon 

Washington,  D.C.  20350 
Attn:  Mr.  A.S.  Rhodes,  Op -964 

Department  of  Statistics 

1 University  of  North  Carolina 

Chapel  Hill,  North  Carolina  27515 
Attn:  Prof.  M.R  Leadbetter 

Southern  Methodist  University 
Department  of  Statistics 

l Dallas,  T xas  75222 

Attn:  Prof.  D.B.  Owen 

Isroel  Institute  of  Technology 
Tcchnion 

1 Haifa,  ISRAEL 

Attn:  Prof.  P.  Naor 


Copies 


\ 


1 


1 


1 


1 


1 


1 


1 


1 


Office  of  Gifts  & Exchanges 
Library 

Texas  A&M  University 
College  Station,  Texas  77843 

Archives 

Texas  A&M  University 
College  Station,  TX  77843 


Copies 

2 

1 


